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SUMMARY 


This is the final technical report on the contract 
NAS9-5884 awarded by the National Aeronautics and Space 
Administration to Rice University for the Charged Particle 
Lunar Environment Experiment (CPLEE) . CPLEE is ion-electron 
spectrometer placed on the lunar surface for the purpose of 
measuring charged particle fluxes impacting the moon from a 
variety of regions and to study the interactions between 
space plasmas and the lunar surface. 

Principal accomplishments under this contract were: 

1. Furnishing design specifications to Bendlx 
. Research Laboratories for construction of the 

CPLEE instruments. 

2. Development of an advanced computer-controlled 
facility for automated instrument calibration. 

3. Active participation in the deployment and 
past-deployment operational phases with regard to 
data verification and operational mode selection. 

4. Publication of numerous research papers, including 
such varied facets as a study of lunar photoelectrons, 
a study of plasmas resulting from man-made lunar 
impart events, a study of magnetotail and magneto- 
sheath particle populations, and a study of solar- 
flare interplanetary particles. 

This report describes these past accomplishments in 
detail along with plans for future research programs using 
data from the instrument. 
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Sclentific Motivation for CPLEE 

The scientific justifications for placing a charged 
particle spectrometer such as CPLEE on the lunar surface are 
quite simply to measure the fluxes, energy spectra, and 
charge types of charged particles bombarding the lunar 
surface and to investigate the interaction of these particles 
and other forms of radiation with the lunar surface. 

In one view the moon is a satellite with an orbital 
radius of 60 R„ (earth radii) that carries the CPLEE instrument 

Ili 

through various regions of physical and scientific interest. 

Figure 1 shows a typical lunar orbit in relation to the 
geometry of the earth’s magnetospheric system and interplanetary 
space. The arrows show the look directions of the two 
charged particle analyzers in CPLEE. It is seen that once 
per lunar month the lunar orbit carries the instrument 
through the distant magnetosheath, magnetotail, and regions 
of direct access to the solar wind and finally (during 
periods of new moon) the instrument is viewing into the 
downstream cavity carved into the solar wind by the moon. 

Thus an opportunity is afforded to investigate a wide variety 
of particle phenomena; charged particles in the geomagnetic 
tail (the plasma sheet), magnetosheath fluxes, fluxes at the 
boundary between the geomagnetic tail and the magnetosheath 
(the magnetopause), the solar wind shock, and direct solar 
wind particles. The relative magnitudes and temporal relationshps 
between these particles in the distant lunar regions and 
particles nearer the earth measured by other satellites can 
be measured. 

In another sense, the instrument is capable of measuring 
the interactions of particle and photon radiation with the 
lunar body itself. The size scale of the moon is comparable 
to or larger than typical scale sizes of the particle fluxes, 
and hence one would expect on occasion significant interactions. 



Figure 1 

A schematic plan of the lunar orbit in relation 
to the earth's magnetospheric system and interplanetary 
space. The arrows show the look direcitons of the instrument 


det-ectors . 




- 3 - 


An example would be a possible shock v/ave near the teminator 
regions due to solar-wind lunar limb interactions. A second 
example would be the photoelectron layer generated by solar 
photons striking the surface, and measurement of these can 
only be made by a surface-based instrument. 

Thus there are numerous scientific justifications and 
expectations in placing a charged particle spectrometer on 
the lunar surface. In subsequent sections of this report it 
will be seen how many of these expectations have been fulfilled. 
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The Evolution of CPLEE 

The concept of the CPLEE instrument grew from an active 

program of auroral sounding rocket research at Rice University 

during the years 1963-1968. One of the objects of the 

auroral research program was to make detailed, high time 

resolution measurements of the auroral charged particle 

spectrum over as wide an energy range as possible. Along 

with this requirement went the usual constraints of low 

weight, power, and size made necessary by the constraints of 

the research vehicle. Accordingly, an instrument code-named 

SPECS (Switching Proton-Electron Channeltron® Spectrometer) 

was developed to meet these requirements. SPECS consisted 

of an electrostatic particle deflection system and an array 

0 

of six Channeltron detectors to separate and detect particles 
according to energy and charge sign. In its final evolutionary 
stages, the instrument made a 15-point measurement of the 
spectra of both ions and electrons w r ith energies ranging 
from iJO electron volts to over 20,000 electron volts. The 
time required to complete a measurement cycle could be 
varied, but for most applications was on the order of a few 
seconds . 

SPECS instruments were flown on each of a series of 
three Javelin sounding rockets in February, 1967 and February, 
1968 . In addition, a SPECS instrument was successfully flown 
on a Rice/ONR auroral research satellite code-named Aurora 
I in June of 1967 and operated successfully for over one 
year. Thus the basic detector design was well-proved in 
space applications and environments well before the actual 
CPLEE lunar deployment, and the cognizant scientists were 
able to gain considerable experience in analyzing data from 
the basic instruments years in advance of the CPLEE mission. 


Q 


Registered Trademark, 


Bendix Corporation 
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Description of the CPLEE Instrument 

The CPLEE consists of a box supported by four legs. 

The box contains two similar physical charged-particle 
analyzers, two different programable high-voltage supplies, 
twelve 20-bit accumulators, and appropriate conditioning and 
shifting circuitry. The total weight on Earth is approximately 
2.7 kg (6 lb.), and normal power dissipation is 3*0 W rising 
to approximately 6 W when the lunar-night survival heater is 
on. ' 

Each physical analyzer contains five C-shaped channel 
electron multipliers with a nominal aperture of 1 mm each 
and one helical channel electron multiplier with a nominal 
aperture of 8 mm. These are shown schematically in Figure 2. 

The channel electron multiplier is a hollow glass tube, 
the inside surface of which, when bombarded by charged 
particles, ultraviolet light, etc., is an emitter of secondary 
electrons. In the CPLEE, the aperture of each electron 
multiplier is operated nominally at ground potential (actually 
at 16 V), while a voltage of 2800 or 3200 V (selected by 
ground command) is placed on the other (i.e., anode) end. 

Thus, if an incident particle enters the aperture and secondary 
electrons are produced, these are accelerated and hit the 
walls to generate more secondary electrons, so that a multipli- 
cation to an order of 10^ is achieved by the time the pulse 
arrives at the anode. After conditioning, pulses from each 
electron multiplier are accumulated in a register for later 
readout as described in the following paragraphs. 

As shown in Figure 2, incident particles enter an 
analyzer through a series of slits and then pass between two 
deflection plates across which a voltage can be applied. 

Thus, at a given deflection voltage, the five small-aperture 
electron multipliers make a five-point measurement of the 



Figure 2 


A schematic drawing of the charged-particle analyzer 
in the CPLEE instrument. Particles pass through a series of 
collimators, and are then deflected by a set of electrostatic 
deflection plates onto an array of six channel electron 
multipliers . 
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energy spectrum of charged particles of a given polarity 
(e.g., electrons), while, simultaneously, the large-aperture 
electron multiplier makes a single wideband measurement of 
particles with the opposite polarity (e.g., protons). The 
advantages of simultaneously measuring particles of opposite 
polarity and of simultaneous multiple-spectral samples are 
considerable in studies of rapidly varying particle fluxes. 

In the CPLEE, the deflection-plate voltage, in the 
normal mode, is stepped in the sequence shown in Figure 3* 

As a consequence, the energy passbands shown in Figure 4 
are sampled. Although data acquired by the six sensors are 
not transmitted simultaneously, the six sensors are connected 
to six accumulators for exactly the same time (viz, 1.2 sec) 
and the contents transferred to shift registers for later 
sequential transmission. 

Two analyzers, A and B, point in the directions shown 
in Figure 5* The same deflection voltage is applied to each 
analyzer simultaneously, with counts from 1 . 2-sec accumulation 
time of analyzer A being transmitted while counts from 
analyzer B are accumulating. Thus, each voltage is normally 
on for 2.4 sec with the result that the total cycle time is 
19.3 sec (Fig. 3 ), when allowance is made for two sample 
times when the deflection voltage is zero. On one of those 
two occasions, counts are accumulated as usual to measure 
background or contaminating radiation. On the other occasion, 
a pulse generator of about 375 kHz is connected to the 
accumulators, to verify operation. 

The command link with the ALSEP provides a variety of 
options for CPLEE operation. Aside from the usual power 
commands common to all ALSEP experiments, three commands are 
provided that allow the normal automatic stepping sequence 
to be modified. The sequence can be stopped and then the 
deflection plate supply can be manually stepped to any one 
of the eight possible levels. This is done to study a 



Figure 3 

The deflection voltage stepping and analyzer 
switching sequence of CPLEE. In the automatic mode, a 
complete cycle is completed in 19*3 seconds. In the manual 
mode, the deflection voltage stepper is halted at one of the 
eight possible levels, and the analyzers are read out 
alternately with a cycle time of 2.4 seconds. 
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Figure 4 

A schematic representation of the energy pass- 
bands of CPLEE . The energy passbands are shown for each 
of the three deflection voltage levels of 35 * 350, and 
3500 volts. 
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Figure 5 


The CPLEE system, showing the locations and look 


directions of the two analyzers. ■ 
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particular phenomenon (e.g., low-energy electrons) with 
higher time resolution (2.4 sec). A second set of commands 
allows the electron-multiplier high-voltage supply to be set 
at either 2800 or 3200 V. The higher voltage is used in the 
event the electron-multiplier gains decrease during lunar 
operations. A third pair 'of commands allows the normal 
thermal-control mode to be bypassed in the event of failure 
of the thermostat, thus offering manual control of the 
heaters. A summary of the command structure and of the 
engineering housekeeping measurements is shown in Figure 6. 

The CPLEE apertures are covered with a dust cover to 
avoid contamination during deployment and, particularly, 
during LM ascent. The dust cover, was made doubly useful 
because a Ni^ radioactive source was placed on the underside 
over each aperture. Thus, the sensors were proof calibrated 
on the Moon, and the data compared with measurements made in 
the same way with the same system when the unit was last 
calibrated on Earth. The results of the Ni^ beta source 
tests at times of pre-calibration, post-calibration, and 
post-deployment are shown in Figure 7. 



Figure 6 

The engineering housekeeping measurements, red- 
line limits, and command structure of CPLEE. 
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Figure 7 

A summary of the beta source test performed 

before and after instrument calibration, and after lunar 
surface deployment prior to dust cover removal. 
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Calibration of the CPLEE Instrument 

Although the process of calibration of a scientific 
detector is usually viewed as a prosaic task, some rather 
challenging technological problems had to be solved in this 
endeavor and therefore it is worthwhile to examine in some 
detail the calibration procedure and the system that was 
developed to accomplish this procedure. 

The object of calibration of any charged particle 
detector is to obtain a number which- relates the counting 
rate from the instrument to the particle flux incident upon 
the instrument. For particle fluxes usually encountered in 
space measurement situations, a convenient unit for particle 
flux is particles/square centlmeter-second-steradian-electron 
volt (part. /cm -sec.-ster .-eV) . This unit then measures the 
number of particles crossing a unit square centimeter 
each second, from a cone of arrival directions subtending 
one steradian in an energy interval one electron-volt in 
width. In mathematical terms, if we let C.R. be the counting 
rate of the detector and j (E) be the incident flux, then: 

CO 

C.R. = ; o j (E)A(E)fi(E)e(E)dE 

where A(E), Q(E), and e(E) are quantities determined by the 
detector geometry and are defined in the following manner. 

A(E) is the effective area of the detector (cm ), n(E) is 
the effective solid angle of the detector, which is related 
to the angular response, and e(E) is the efficiency of the 
detector, expressed in counts/particle. Thus if one can 
obtain a plot of the product A(E)ft(E)e(E) as a function of 
energy, then one can in principle unfold the integral relation 
shown above to obtain j (E) if the counting rate is known. 
Parenthetically, the plot of the product A(E)fi(E)e(E) vs. E 
is commonly known as the energy passband of the detector. 
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The task of calibration is then to determine A(E), 

(E), and (E) either singly or in product form. Conceptually 
this would appear simple since one has only to create a 
particle beam of known intensity, impress it upon the detector, 
and observe the counting rate. The problem arises however 
when one considers that one is interested in the response to 
beams distributed in arrival angle whereas beams available 
in the laboratory are most easily generated unidirectionally , 
that is all particles have the same direction. Therefore, 
if a unidirectional beam is applied to the instrument from a 
given arrival direction and , then one actually measure 
the product A(E) (E) at a given energy and arrival direction. 

To obtain (E) it is then necessary to repeat the measurement 
for a number of arrival directions sufficient to include the 
complete angular range of instrument sensitivity. 

Now the magnitude of the calibration problem can be 

appreciated. In CPLEE there are two analyzers that each make 

an 18-point . spectral measurement. Thus, for electrons only 

(say), there are 36 separate energy passbands that must be 

determined. It was desired to take 10 energy points per 

passband, and so this involves 360 separate energy runs. 

Each one of these runs in turn requires measurement at 

numerous angles. In particular, the desired angular grid turned 

out to be 10° x 20° in 0.5° steps each axis, requiring 86l 

q 

separate angle steps. This works out to more than 3 x 10 
separate measurements to calibrate one CPLEE for electrons 
only. Considering that the program called for calibration 
of four units (prototype, two flight units, and flight 
spare) plus calibrations of each for sensititivty to solar 
ultra-violet radiation, it quickly became obvious that 
traditional means involving manual setting of angles, reading 
data, and hand reduction could not possibly fulfill the 
requirements. Accordingly, the decision was made to develop 
a computer-controlled calibration facility that could function semi 
automatically with a minimum of operator attention required. 
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The basic requirements for such a calibration system 
were then: 1. A stable source of electrons capable of 

producing beams with energies ranging from 40 eV to over 
20,000 eV; 2. A mechanical fixture to position the instrument 
at various angles relative to the beam; 3. Fixture drive 
motors and angle position indicators capable of computer 
control and readout; 4. A computer system capable of posi- 
tioning angles, reading data, storing data, and monitoring 
various engineering parameters in an automatic sequence; 5- 
A vacuum chamber of sufficient size to contain the electron 
gun, fixture, and instrument; 6. A number of competent 
personnel to make it all function. 

The various required components, included a SDS-92 
computer, a Yarian vacuum system, and numerous items such as 
power supplies, etc. were purchased and assembled into the 
calibration system. The mechanical accessories (i.e., the 
fixture) were fabricated at Rice and also the logic circuitry 
necessary to interface the computer with the CPLEE, the 
fixture drive motors and angle indicators; and various other 
measuring devices was designed and built at Rice. The 
system became operational in August 1968. It was found that 
approximately one month was required for a complete particle 
and ultra-violet sensitivity calibration of one CPLEE. 

An example of the calibration results for CPLEE is 
shown in Figure 8, showing the analyzer A electron passbands 
for. all eighteen energy channels. These were determined from 
the calibration of CPLEE SN-5, the flight unit for the 
Apollo 14 ALSEP. Three passbands are shown for a given 
channel, corresponding to each of three values of electrostatic 
deflection voltage; 35, 350, and 3500 volts. It is seen 
that the passbands for a given channel are nearly identical 
when the factor of 10 changes in the energy scale are 
considered. In practice, these passbands were numerically 
integrated to obtain a number GFo x aE, with GFo x aE = 

*GF(E)dE and GF(E) = A(E)p (E) e (E) . GFo was set to be numerically 
equal to the peak value of GF(E), and therefore AE- is an 
effective passband width. Operationally then, the fluxes • 

were determined from the relation j(Eo) = C.R./(GFo x AE). 



Figure 8 

The electron energy passbands of Analyzer A, 
CPLEE S/N 5 as derived from the calibration. 
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Deployment, Initial Checkout, and 
Operational History 

The CPLEE was deployed with no dificulty at approxi- 
mately 18:00 G.M.T. on February 5 , 1971. Leveling to within 
2.5° and east-west alinement to within + 2° were to be 
accomplished with a bubble level and a Sun compass, respectively. 
A photograph of CPLEE on the lunar surface is shown in Figure 9« 

It has since been determined by' a careful study of the 
lunar photographs and a comparison of predicted and actual 
solar ultraviolet response profiles that the experiment is 
1.7° off level, tipped to the east, and 1° av. T ay from a 
perfect east-west alinement. This error is well within the 
preflight specifications. 

The CPLEE was first commanded on at 19:00 G.M. T. , 

February 5 > 1971 during the first period of extravehicular 
activity for a brief functional test of 5-min duration. All 
data and housekeeping channels were active, and the instrument 
began operation in the proper initial modes (i.e., automatic 
sequencer, on; electron-multiplier voltage increase, off; 
and automatic thermal control, on). 

A complete instrument checkout procedure was initiated 
at 04:00 and continued until 06:10 G.M.T. , February 6, 1971. 
During this period, data from the dust-cover beta sources 
were accumulated and compared with prelaunch calibrations. 

These tests showed no significant changes in the instrument 
calibration. Also during this period, all command functions 
of the CPLEE were exercised except the forced heater mode 
and dust-cover removal commands. The instrument responded 
perfectly to all commands. After the checkout procedure, 
the CPLEE was commanded to the standby mode to await LM 
ascent . 

Following LM ascent, the CPLEE was commanded on at 
19:10 G.M.T. and the dust cover was successfully removed at 



Figure 9 


A photograph of the CPLEE instrument on the lunar 
surface. The dust cover is in place. The bubble level is 
visible adjacent to the letter "E" and the sun compass > read 
by 'the position of the handling tool shadow, is visible 
at the opposite end of the instrument. 
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19:30 G.M.T., February 6. The CPLEE immediately began 
returning data on charged-particle fluxes In the magnetosphere. 

The instrument temperatures were carefully and continuously 
monitored for 45 days after deployment. It was found that 
the temperature range was nominal, with the internal electronics 
temperature ranging from 58° C at lunar noon to -24° C 
during lunar night. The total lunar eclipse of February 10 
offered an excellent opportunity to determine various thermal 
parameters and to test the capability of the CPLEE to survive 
extreme thermal shocks. A plot of the physical-analyzer 
temperature during the eclipse is shown in Figure 10. The 
maximum thermal shock occurred after umbra exit, with a 
temperature change rate of 25° C/hr. Also from this figure, 
it is possible to derive a thermal time constant of approxi- 
mately 1.9 hr. The CPLEE suffered no ill effects from this 
period of rapid temperature changes. The command capability 
of the CPLEE was used extensively during the 45-day real- 
time support period to optimize scientific return from the 
instrument. Alternate 1-hr periods of manual operation at 
the -35 V step and automatic operation were used to concentrate 
on rapid temporal variations in low-energy electrons. 

Similarly, alternate periods of 350 V manual and automatic 
operation were used to focus on rapid changes in magnetopause 
ions and the solar wind. In fact, the manual operation 
capability and the attendant 2.4-sec sampling interval made 
possible the detection of phenomena that would have been 
impossible to detect otherwise because of sampling problems 
and aliasing. Most of the decisions concerning operational 
modes were based on viewing the real-time data stream. 

On April 8, 1971 at 21:55 the housekeeping monitor for 
the Analyzer B Channeltron high voltage supply showed a 
serious undervoltage condition. Subsequent analysis of the 
data tapes showed that this was a sudden failure, the voltage 
having decreased from its nominal value of 2800 volts to 
approximatley 800 volts within one housekeeping data cycle 



Figure 10 


A plot of the temperature profile of CPLEE during 


the lunar eclipse of February 10, 1971. 
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period. The instrument was placed into standby for the 
remainder of the lunar day period and operation with Analyzer 
A only was resumed on April 16, 1971- A meeting was held 
with representatives of JSC, Bendix, and Rice University 
where it was determined that the most probable cause of the 
failure was either a failure of a component in the high 
voltage rectifier assembly or an arc track developing in the 
high voltage supply encapsulating material. 

The instrument continued to operate returning Analyzer 
A data only until June 6, 1971. At that time a low voltage 
condition (1800 volts) appeared in the Analyzer A Channeltron 
high voltage supply, rendering the Channeltrons inoperative. 

The decision was made to place the instrument into standby 
and continue testing during the subsequent lunar night 
period . 

It was discovered that the instrument would operate for 
short periods during the lunar night (1/2 hour) before the 
high voltage decreased to the point where Analyzer A was 
inoperative. This pointed to a temperature-dependent effect 
wherein normal voltage could be developed when the instrument 
was cold but local heating effects would soon cause degradation. 

It was also learned that, as time went on, longer and longer 

periods of lunar night operation were possible and beginning July of 

1972 operation during the entire lunar night was possible. 

Furthermore, operation was possible for portions of the 

lunar day. In December, 1972 at the commencement of the 

real-time support period for the Apollo 16 mission the 

instrument was turned on at lunar morning and continued to 

operate continuously throughout the lunar day. However, in 

April, 1973 a second partial failure of the high voltage 

supply occurred which again restricted operation to lunar 

night periods only. This condition has continued to the 

present day. 
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Scientific Results of the CPLEE Program 

In this section we discuss the scientific results of 
the CPLEE program. As this report includes reprints of most 
of the CPLEE scientific publications, only summaries will be 
presented here. 

A prominent feature of data obtained while in the 
geomagnetic tail was a stable, low-energy electron popula- 
tion with energies ranging from 40 to 200 electron volts. 

Although possibly due to a low energy magnetospheric population, 

a plot of their flux during the lunar eclipse of February 

10, 1971 revealed their origin to be of photoelectric 

origin (Figure 11). Their correlation with sunlight upon 

the lunar surface shows that these fluxes were due to photoelectrons 

generated by solar photons striking the lunar surface. A 

plot of their energy spectrum is shown in Figure 12, showing 

that a) the photoelectrons were essentially isotropic over 

the upper hemisphere and b) the highest energy was about 200 

electron volts. Recalling that the instrument detects those 

photoelectrons which are returning to the surface, this 

means that the sunlit lunar surface potential when the moon 

is in the magnetotail away from the plasma sheet is on the 

order of 200 volts with respect to the surrounding medium. 

The impact of the Apollo 14 Lunar Module offered an. 
opportunity to study any charged particle fluxes resulting 
from the impact. The impact point was 66 km west of the 
CPLEE location. Figure 13 shows a time history of particle 
fluxes resulting from the impact. This plot shows fluxes of 
70 ev ions and 65 ev negative particles. It is seen that 
approximately one minute after impact, two plasma '’clouds" 
passed the instrument. These ions could not have been generated 
at the moment of impact, for had that been the case the 
particles would have appeared less than one second after 
impact. Rather, it was postulated that the impact produced 
expanding neutral gas clouds, and constituents of this cloud 



Figure 11 


The counting rate of Channel 6, +35 volts sensitive 
to electrons with 50 eV< E< 150 e.V during the lunar eclipse 
of February 10, 1971* The disappearance of the fluxes 
during the eclipse shows that the electrons are of photo- 
electric origin. 
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Figure 12 


The lunar surface photoelectron spectrum measured 
by CPLEE. The geometries relative to the surface and to the 
incoming photon flux are shown in the inset. 
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Figure 13 

Ion and electron fluxes resulting from the impact 
of the Apollo lM Lunar Module upon the lunar surface. 
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were then Ionized and energized by interaction with the 
solar wind. The implications of the results are that if 
neutral gas is released upon the lunar surface by any means 
then it could subsequently be carried away relatively rapidly 
by strong interactions between the gas and the solar wind. 

A systematic study of particle fluxes in the geomagnetic 
tail, the plasma sheet, was undertaken in order to determine 
the average characteristics of these particles and their 
response to geomagnetic storm activity. One question to be 
considered is that of shadowing, that is "does the moon 
effectively sweep out magnetic flux tubes of their particles 
and thus prevent CP LEE from observing these particles?" 

This question was answered by comparing CPLEE particle 
observations with magnetic field data from the Explorer 
35/Ames Research Center lunar-orbiting magnetometer. The 
magnetometer indirectly indicates the presence of the plasma 
sheet by a decrease in the magnetic field strength owing to 
the , diamagnetism of the plasma sheet particles. It was 
discovered that in 90 % of the cases when magnetic field 
depressions were observed, particle fluxes were observed by 
the CPLEE instrument. Thus for the majority of the time 
plasma sheet fluxes are not shadowed from the lunar surface. 

The average plasma sheet electron and ion spectra are 

shown in Figure 14. The lowest energy portion of the spectrum 

is due to the lunar photoelectrons. The plasma sheet spectra 

represent electron and ion populations with n = 0.1, T = 

0 0 

200 eV, n^. = 0.1, T^ = 2.5 keV. These particle temperatures 
were found to be equivalent to temperatures measured in the 
plasma sheet at 20 Rg by the Vela satellites, but the number 
density was on the average a factor of 5 less than that 

measured at 20 R . 

0 

A statistical study of plasma sheet encounters was done 
in order to determine the vertical extent, or thickness of 
the plasma sheet at the lunar distance. The thickness was 
determined to be 5 + 2 Rg, somewhat thinner than the value 



Figure 14 


The average plasma sheet electron and ion spectra 
observed at the lunar distance. The lowest energy portion 
of the spectrum was due to photoelectrons, and had to be 
subtracted from the data in order to obtain the fitted spectra 


shown 
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of 6-12 Rg measured in the range 20-30 Rg from the earth by 
other satellites. 

The response of the plasma sheet to geomagnetic substorms 
has also been investigated. It was found that at the onset 
of a substorm a sudden, dramatic enhancement in the plasma 
sheet flux appeared at the lunar orbit followed by a disappearance 
of plasma until substorm recovery. The thinning, or disappearance 
of the plasma sheet at substorm onset is a phenomenon observed 
at distances closer to the earth, but the precursor spike 
seems to be unique to the lunar distance and is quite possibly 
a result of a pileup of plasma along the leading edge of the 
thinning front. 

Strong flows were observed in the netural sheet on two 
occasions during substorms. At the time of substorm onset 
the magnetic field normal to the neutral sheet was observed 
to change direction and simultaneously a strong anti-solar 
ion flow was observed. This is strong evidence that at the 
time of substorms a neutral line, or a line of magnetic 
field connection across the neutral sheet, is formed between 
the earth and the moon. Plasma is accelerated away from the 
neutral line, travelling both toward the earth and backwards 
toward the moon. 

On April 9, 1971 a large world-wide magnetic storm 
occurred. Data from CPLEE, the Solar Wind Spectrometer 
Experiment (SWSE), the Explorer 35/ARC Magnetometer, and 
solar wind and magnetic field detectors on-board the Explorer 
33 satellite were used for a detailed study of the response 
of the distant geomagnetic tail to this storm. Figure 15 
shows the geometry of the magnetosphere, the location of the 
solar wind discontinuity prior to the storm, and the locations 
of the detectors. 

The storm was first signaled by ground-based magnetometer 
disturbances at 0428 G.M.T. At 0543 magnetosheath particle 
fluxes appeared at the moon and remained until after 1000. 
Referring to the pre-storm geometry of Figure 15, it Is seen 


Vi 



Figure 15 


The geometry of the earth's magnetosphere relative 
to the solar v/ind discontinuity and the locations of the 
moon and the Explorer 33 satellite immediately prior to 
the geomagnetic storm of April 9, 1971. 
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that the entire geomagnetic tail was compressed until the 
magnetopause boundary crossed the position of the rnoon. The 
radius of the geomagnetic tail decreased from 26 to 16 R E . 

This was accompanied by an increase in the tail magnetic 
field strength from 10 gammas to 32 gammas. Multiple crossings 
of the magnetopause boundary after 1000 G.M.T. coincided 
with changes in the angular direction of the solar wind 
flow, showing the influence of solar wind flow direction 
upon the orientation of the geomagnetic tail. Finally, when 
the moon was in the magnetosheath between 0545 and 1000 
G.M.T. an electron population characteristic of the plasma 
sheet was found to be superimposed upon the normal low- 
energy magnetosheath electron population. This shows an 
enhanced loss of plasma sheet particles into the magnetosheath 
during periods of high magnetic activity. 

The large solar flare event of August, 1972 was an 
event of unusual magnitude and fortunately at the time the 
moon. was outside of the geomagnetic tail regions in interplanetary 
space. Therefore there existed the opportunity of studying 
particle fluxes resulting from the flares unmodified by 
interaction with the earth's magnetic field. Plots of the 
counting rates of analyzer A, channel 1 at the deflection 
voltages +0 (background), -35 volts and -350 volts for the 
period is shown in Figure 16. The count rate scales are 
displaced by a factor of 10 vertically for the sake of 
clarity. It is seen that after day 217 at 1200 hours the 
counting rates were independent of deflection voltage, 
indicating that these counts were due to high energy solar 
flare particles that were capable of penetrating the instrument 
case. This interpretation was confirmed by cosmic ray 
detectors on other satellites. However, between 0200 and 
about 1200 on day 217 there are counting rates in the -35 
and -350 channels that are not matched by corresponding 
rates in the background channel, indicating a population, of 
electrons present. A temporal history of the energy spectra 



Figure 16 

Counting rates in CPLEE Analyzer A channel 1 at 
three deflection voltages during the solar flare events of 
August 1972. The plots have been displaced by a factor of 
10 vertically for clarity. The identical counting rates at all 
deflection voltages after 1200 GMT on day 217 shows that 
these counts were due to high energy protons penetrating the 


instrument case. 
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of these electrons is shown in Figure 17 . Here the horizontal 
scale is energy and the vertical scale is flux. The spectra 
represent averages over sixteen minutes. The most likely 
explanation for these electrons is that they were generated 
at the earth's bow shock by the interaction between the 
shock and the strongly disturbed solar wind, and that these 
electrons subsequently propagated upstream to the location 
of the moon. 

In summary, the principal scientific results of the 
CPLEE program were: 

1) Observation of a layer of lunar photoelectrons above 
the sunlit lunar surface. From these data the lunar 
surface photo yield and- potential in the high- 
latitude geomagnetic tail were calculated. 

2) Observation of a strong interaction between the 
solar wind and neutral gas clouds produced by 

the Apollo 14 Lunar Module Impact. This interaction 
produced charged particles with energies ranging 
up to 100 eV. 

3) Determination of the average spectral and spatial 
characteristics of the plasma sheet at the 
lunar distance. 

4) Discovery of neutral line formation between the 
earth and the moon and of strong anti-sunward plasma 
flow during magnetic substorms. 

5) Determination of the response of the distant 
geomagnetic tail and magnetosheath to the geomagnetic 
storm of April 9, 1971* 



Figure 17 


Electron spectra observed in August 1972 resulting 


from the solar flare events. 
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6) Observation of a hot electron gas in interplanetary 
space resulting from the August, 1972 solar flare 
event . 

7) Discovery of an energetic electron population in the 
dawn-side magnetosheath that is strongly influenced 
by geomagnetic activity. These particles were 
shown to be capable of being a source of plasma 
sheet particles. 

Only summaries of the scientific results have been 
given here. For further details the reader should consult 
the publication reprints in Appendix B. 
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Conclusion 

The CPLEE program has been seen to have been an outstanding 
technical and scientific achievement. Of no small consequence 
has been the educational benefits of the program. One 
student (Frederick J. Rich) received his Ph.D. degree in 
1973 as a result of CPLEE data analysis and another student 
(Patricia R, Moore) will receive a M.S. degree in the spring 
of 1974 and will continue research into CPLEE data leading 
to the Ph.D. degree in 1975- Furthermore, the technical 
capability of our laboratory was significantly upgraded by 
the demands of the program. 

Analysis of CPLEE data is by no means complete. Small 
but significant fluxes of particles bombard the lunar surface 
throughout the lunar night when the instrument is viewing 
into the solar wind cavity downstream of the moon. It is, 
of course, the natural proclivity of an investigator to look 
at data with the largest counting rates first, and so these 
lunar night data have been largely neglected. However, 
studies of these data in conjunction with ion data from the 
Suprathermal Ion Detector (SIDE) and from the Explorer 
35/ARC magnetometer hold promise for a deeper understanding 
of interaction between the moon and the solar wind and of 
particle energization processes at the earth’s bow shock. 

These investigations have been proposed and accepted for the 
Post-Apollo Program of Data Analysis and Synthesis. 
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ABSENCE OF THE PLASMA SHEET AT LUNAR DISTANCE 
DURING GEOMAGNETICALLY QUIET TIMES 

WILLIAM J. BURKE and DAVID L. REASONER 

Department of Space Science, Rice University, Houston, Texas 77001, U.S.A. 

{Received 18 October 1971) 

Abstract— Plasma data from the Apollo XIV Charged Particle Lunar Environment Experiment 
(CPLEE) are presented to show that contrary to previously published analyses the plasma sheet 
does not extend to the lunar orbit with a thickness of 8 R B . Two electron spectral types are 
observed: (1) low energy photoelectrons with no statistically significant medium and high 
energy fluxes, and (2) double peaked medium and high energy electrons. The second type is 
observed either coincident with auroral substorms or at the center of the tail during quiet times. 

These spectra are one to several orders of magnitude less intense than plasma sheet spectra 
measured near 20 Re. 

Detailed analyses of the spatial, temporal and energy characteristics of the near-Earth 
plasma sheet (<30 R E ) have been carried on by many investigators, (cf. review articles by 
Gringauz (1969), Vasyliunas (1970) and references therein, as well as Hones et aL (1971)). 
Preliminary reports have been presented by Nishida et aL (1969), concerning magnetically 
quiet periods, and by Prakash and Binsack (1971), concerning a disturbed period, extending 
these observations to 60 R E using plasma measurements made by the M.I.T. Faraday cup 
aboard Explorer 35. Both conclude that the plasma sheet does in fact extend to lunar 
distances in approximately the same form and characteristics as at closer distances. 

In this paper we report data from four months of observation by the Moon based 
detector CPLEE (Charged Particle Lunar Environment Experiment). We conclude that 
during geomagnetically quiet times the plasma sheet does not extend to the Moon. 
Although our conclusion is consistent with the disturbed time report of Prakash and 
Binsack (1971), we find no consistancy with that of Nishida et aL (1969). 

DESCRIPTION OF THE INSTRUMENT 

CPLEE was deployed on the lunar surface at lunar coordinates 3° 40' S latitude, 
17° 27' W longitude by the Apollo 14 astronauts on February 5 1971. The instrument has 
been described in detail by O’Brien and Reasoner (1971), but the salient features relevant to 
this study will be repeated here. 

The instrument contains two identical charged particle analyzers that are similar to the 
particle analyzer SPECS described by O’Brien et aL (1967). The look direction of one 
analyzer (Analyzer A) is toward the local lunar vertical, and the look direction of the other 
(Analyzer B) is 60° from vertical, toward lunar west. As the deployment site is 2° south 
lunar latitude, the look directions remain within 3-5° of the plane of the ecliptic. 

Each analyzer contains a set of collimators and electrostatic deflection plates to separate 
particles according to energy and sign of charge. Particles are detected by an array of five 
1 mm diameter channel electron multipliers placed to one side of the central axis of the 
deflection system, and a single large aperture (8 mm diameter) multiplier placed on the 
opposite side of the central axis. Thus for a given voltage on the deflection plates, the five- 
multiplier array measures particles of a given charge sign (e.g. electrons) in five differential 
energy ranges and the single large aperture multiplier measures particles of the opposite 
charge sign (e.g. ions) in a single wide-energy range. 
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A stepping high-voltage supply furnishes ±3500, ±350, ±35, and 0 V in sequence to 
the deflection plates. There are two steps at zero volts. At one a background measurement 
is made while a test oscillator signal is introduced at the other for data link verification. The 
normal 8-step sequence requires 19*2 seconds to complete. However, the automatic 
sequence can be interrupted by ground command, and any of the 8 deflection voltage levels 
can be selected. In this way, a particular range of particle energy can be observed with a 
time resolution of 2*4 seconds. 

The energy passbands of the various channels at the three deflection voltage levels (3500, 
350, and 35 V) are shown in Fig. 1. Note that the total energy range is 40 eV to 50 keV. 



PARTICLE ENERGY ( keV) 


Fig. 1 . Rectangular equivalent energy passbands of CPLEE at the three deflection 

VOLTAGES. 


OBSERVATIONS 

Although in this section we show only the first month of CPLEE observations during 
the magnetotail passage of February 1971 in some detail, they are typical of data from 
subsequent observations. The YZ projection of the Moon’s trajectory has been plotted in 
both solar ecliptic and solar magnetospheric coordinates from the beginning of February 7 
to the end of February 13 in Fig. 2. At approximately 0300 hours UT on February 8, 
1971 CPLEE passed from the dusk magnetosheath into the magnetosphere. After several 
encounters with the dawn-side magnetopause it left the magnetosphere at approximately 
0830 U.T. on February 12. The four-day period of tail transit was marked by few magnetic 
fluctuations on Earth, where ^K. p for 24 hour periods ranged between 11 and 19. We are 
thus able to examine the plasma conditions at 60 R E while the magnetosphere was in both 
quiescent and moderately active states. 

Crossings from the magnetopause into the magnetosphere are marked by sharp decreases 
in the counting rate of all CPLEE channels, particularly in those sensitive to low energy 
electrons (40-200 eV). The low energy electron fluxes fall to relatively low and stable levels. 
During the entire passage of the magnetotail the low-energy electron fluxes remained at 
these levels except in short-lived (<1 hr) bursts when the count rate increased in these 
channels by a factor of about 2. Low energy fluxes in the magnetotail are found to be iso- 
tropic i.e. the fluxes measured in detector A are substantially the same as in B. We assert 
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Fig. 2. The YZ projection of the Moon’s trajectory from the beginning of February 8 

TO THE END OF FEBRUARY 12. 

The beginning of each day is indicated by an arrow. Symbols □ represent solar ecliptic 
coordinates and are spaced 6 hours apart. Symbols + represent solar magnetospheric 
coordinates at 2 hour intervals. The solar ecliptic boundaries between the dawn and dusk 
magnetosheaths and the magnetotail are also indicated. 

that the stable low-energy electrons seen during quiet times (at the lunar surface) in the 
deep tail are photo-electrons and not part of an ambient plasma. In Fig. 3 we have plotted 
the counting rate of Analyzer A, Channel 1 at a deflection voltage of “35 V, sensitive to 
electrons with E w 40 e V for all of February 10, a day encompassing a lunar eclipse (500- 
1000 UT). Note that as the lunar surface in the vicinity of CPLEE is shadowed from the 
Sun during lunar eclipse, the electron counts fall to near zero levels and return as soon as the 
Sun reappears. If the counts were predominantly due to ambient plasma they would be 
observed despite the occultation of the Sun (O’Brien and Reasoner, 1971). 

Figure 4 shows a record of the February magnetospheric crossing. The counting rates 
in the Analyzer A, Channel 3 at “350 V (~750 eV electrons) are averaged over 5 minutes 
then recorded at the closest integral value, A similar record of the same channel at —3500V 
(~7-5 keV electrons) shows the same features. A cursory glance reveals that the uniformity 
of low energy electron counts stands in marked contrast to medium and high energy count 
rates. Medium and high energy electron fluxes reaching the surface of the Moon seem to be 
of two basic kinds. The first kind marks quiet times in which the detectors see fluxes within 
one standard deviation of the background level of the instrument. We refer to times such as 
February 8, hours 16-21; February 9, hours 8-13; February 10, hours 3-13; and February 
3 1 , hours 2-9. For periods of up to ten hours in length, the counting rates remained at or 
near their background levels. Departures from background to the second spectral type are 
of relatively short duration ( < 1 hr). The rise times from background to peak count rates 
are generally between 2 and 10 minutes. Although these fluxes generally appear first in 
detector A, isotropy is established within a few minutes. 

Phenomenologically, these enhanced spectra seem to appear under two quite different 
conditions. We have examined ground magnetograms from College, Fort Churchill, and 
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HOURS 

Fig. 3. CPLEE Analyzer A, Channel 1 at —35 V deflection for February 10, 1971 1 

ILLUSTRATING THE STABLE PHOTO-ELECTRON BACKGROUND WHICH DISAPPEARS AS THE MOON 
GOES INTO ECLIPSE FROM ABOUT 500 TO 1000 UT. 

Kiruna, and have found that the relatively short bursts (<1 hr) occur within 15 minutes of 
the onset of a substorm observed at one or more of these stations. Although they were 
measuring positive ions, we find our observation of enhanced plasma fluxes at the lunar 
surface associated with substorms on Earth consistent with that of Garrett et aL (1971). 
The longer period of enhanced flux, February 9, hours 17-24, occurred during a period of 
geomagnetic quiet, but the Moon was close to the center of the tail. It is tempting to assume 
that we are observing neutral sheet plasma. However, in the absence of local magnetic 
field data, this remains only an hypothesis. It is important to note that the Moon moved 
from Z s ^— 0*95 or Z HM = 3-46 at 1700 hours to Z SE = 0*65 or Z SM = 2*59 at 2400. 
The subscripts SE and SM indicate solar ecliptic and solar magnetospheric coordinates 
respectively. In neither system of coordinates did CPLEE detect plasma even 1 R K thick, 
much less 8 R E . 

In Fig. 5a we have plotted a typical spectrum of quiet time electron observations. The 
heavy black line indicates a background flux corresponding to one count per second 
Quiet time low energy electrons (<200 eV) are almost exclusively photoelectrons; medium 
and high energy fluxes hover near the background of the detector. 

The second spectral type. Fig. 5b ; shows a slight increase in low electron fluxes, with 
the major enhancement in the higher energy channels. Peak differential fluxes near 500 eV 
have been observed to rise as high as 10 4 electrons/cm 2 -sec-ster-eV. Generally too, we find 
a second peak near 5 keV of ~10 3 electrons/cm 2 -sec-ster-eV. 

The shape of our ‘background spectrum’ (Fig. 5a) suggest that the second peak may well 
be due more to our detector than to nature. For two reasons, we argue that this is not the 
case: 

(1) We have observed spectra where the second peak is absent or attenuated. 

(2) As part of the calibration of CPLEE, Ni 63 beta sources of known strength were 
placed at the apertures of detectors A and B, (O’Brien and Reasoner, 1971). The Ni 63 
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Fig. 4. Four days of data from the February 1971 magnetotail passage by CPLEE, 
Analyzer A, Channel 3, —350 V. 

There is a two hour data gap, 1000-1200 UT on February 8. From ~10;00 of February 11 
through the end of the day the detector was looking at the Sun, putting the contamination 
count rate well above electron counting rates. 


spectrum calculated from the counting rates and the independently-measured geometric 
factors matched fluxes observed by other investigators (Archuleta and DeForest, 1971). If 
the secondary peak were due to an instrumental defect we would have found a similar 
deviation from the Ni* 3 spectrum. 

DISCUSSION AND CONCLUSIONS 

The question thus remains; how do our observations compare with known plasma sheet 
measurements? 

Hones et at. (1971) have presented a synoptic view of plasma sheet and auroral zone 
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Fig. 5(a). Typical CPLEE photo-electron spectrum taken from February 10 ~ 400 LIT. 
Fluxes in A and B were approximately the same. The background flux is that which would be 
observed from one background count in each channel. During quiet time* the medium and 
high energy counts were at or dose to this level. 



Fig. 5(b) Typical enhancement event observed by CPLEE, Analyzer A. 
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spectra in Fig. 1 1 of their article. For the purpose of comparison we have reproduced some 
of their results along with our own. Plasma-sheet electron spectra at 1 7 R E are of two kinds : 
(1) cool electrons, as in spectrum A of Fig. 6. This spectrum peaks at >300 eV at flux 



ENERGY {eV) 

Fig. 6. Enhanced Spectrum at the Moon compared with plasma sheet and auroral 

zone spectra. 

Spectra A and B are ‘cold’ and ‘hot’ electrons observed by VELA in the quiet and post-storm 
inflation of the plasma sheet. Spectrum C is the enhanced CPLEE measurement. Spectrum D 
is a double peaked electron flux observed in an auroral arc. The background flux of VELA and 

CPLEE are also indicated. 

levels ~10 5 electrons/cnr-sec-ster-eV. The spectrum is Maxwellian down to the background 
level of the VELA detectors and shows no high energy tail. (2) Hot electrons, spectrum B, 
appear in the post- storm inflation of the plasma sheet. These fluxes peak near 1 keV at 
levels of several times ID 4 electrons/cnr-sec-ster-eV. No trace of a secondary peak was 
reported. The peak fluxes of this ‘hot’ spectrum is several times the height of the highest 
flux we have observed at the Moon. 

We have displaced the spectrum observed by Westerlund (1970a) because of the character- 
istic double peak found in auroral arcs. There however, the spectral similarity ends, as we 
are below Westerlund’s flux level by three or four orders of magnitude. Since the auroral 
oval is the locus of points in the ionosphere mapping into the plasma sheet (Vasyliunas, 
1969), if the plasma sheet extended to the Moon we might expect to find spectral similarities. 

Our observations indicate that during quiet times the plasma sheet does not extend out to 
the Moon. Observed plasma has spectral peaks near 0*5 and 5 keV. However the fluxes are 
one to several orders of magnitude lower than near Earth plasma sheet measurements. The 
upper limit to the quiet time plasma flux at the Moon is 5 x 10 2 electron$/cm 3 -sec-ster-eV 
at E — 0*5 keV. These enhancements observed to be coincident with auroral substorms 
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could be interpreted as a temporary expansion or outward flow of the plasma sheet to 60 R E . 
The single observation of enhanced fluxes during a quiet time (February 9, hours 17-24) 
could have been due to the plasma sheet. If this is the case however, its thickness would be 
less than 1 R E . 

Note added in proof 

Since the submission of this article a number of people have argued that even if the plasma sheet extended 
to 60 R E ' CPLEE might not be able to detect it because of shadowing. A preliminary study of ion fluxes at 
the surface of the Moon, as well as a calculation allowing for motion of the plasma sheet due to wagging of 
the tail, lead us to conclude that lunar shadowing is not a dominant phenomenon in this experiment. 
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Measurement on the Lunar Surface of Impact-Produced Plasma Clouds 
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Simultaneous enhancements of low-energy ions and negative-particle fluxes due to the im- 
pact of the Apollo 14 lunar module were observed by the lunar-based charged-particle lunar- 
environment experiment (CPLEE) * The impact occurred 66 km away from CPLEE, and the 
time delay between impact and flux onset was approximately 1 min. It is argued that the ob- 
served charged particles could not have energized at the instant of impact but rather that the 
impact produced expanding gas clouds and that constituents of these clouds were ionized and 
accelerated by some continuously active acceleration mechanism. It is further shown that the 
acceleration mechanism could not have been a static electric field but rather is possibly a con- 
sequence of interaction between the solar wind and the gas cloud. 


The ascent stage of the Apollo 14 lunar 
module Antares impacted on the lunar surface 
on February 7, 1971, at OOh 4.5m 24s GMT. 
Shortly after the impact, a lunar-based charged- 
particle detector based 66 km away detected 
fluxes of low-energy positive ions and negative 
particles with intensities a factor of 10 greater 
than the ambient fluxes. The ion and electron 
enhancements exhibited near-perfect temporal 
simultaneity, and we report here preliminary 
studies of these impact-produced plasma clouds. 

The measurements were made with the 
charged-particle lunar-environment experiment 
(CPLEE) deployed as part of the Apollo 14 
Alsep instrument array at Fra Mauro. The 
CPLEE instrument is conceptually similar to 
the switching proton-electron channeltron spec- 
trometer (Specs) described in detail by O’Brien 
et al . [1967]. Two identical particle analyzers 
are housed in the unit. One analyzer, labeled 
A, is pointed toward the local vertical, and the 
other analyzer, labeled B, is pointed 60° from 
vertical toward lunar west. 

We refer the reader to O’Brien et al. [1967] 
for a detailed description of the particle anal- 
yzers and report here a few salient features of 
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the instrument relevant to this report. Charged 
particles are deflected by a set of electrostatic 
deflection plates according to energy and charge 
sign into the apertures of an array of six chan- 
nel electron multipliers, and, at a given deflec- 
tion-plate voltage, an analyzer makes meas- 
urements of fluxes of particles of one charge 
sign (e.g., electrons) in five energy ranges and 
particles of the opposite charge sign (e.g., ions) 
in a single energy range. Normally the instru- 
ment steps through a series of six deflection 
voltages plus two background steps every 19,2 
sec. However, the automatic sequence can be 
halted by ground command and the deflection 
voltage stepped to any one of the eight levels, 
with a consequent reduction of the sampling 
interval to 2.4 sec. Before the impact, the de- 
cision was made to operate the instrument in 
the manual mode at a deflection voltage where 
the instrument was sensitive to negative par- 
ticles in five energy ranges centered at 40, 50, 
65, 95, and 200 ev, respectively, and sensitive 
to positive ions in a single energy range with 
peak response at 70 ev and half -intensity points 
at 50 and 150 ev. As will be seen, this decision 
proved extremely fortuitous. 

The Antares impact occurred at lunar co- 
ordinates 3.42°S latitude and 19.67 °W longi- 
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tude, a point 66 km west of CPLEE, at OOh 
45m 24s GMT on February 7, 1971. The geom- 
etry of the impact event is shown in Figure 1. 
The terminal mass and velocity were 2303 kg 
and 1.68 km/sec, respectively, resulting in an 
impact energy of 3.25 x 10* joules (G. V. 
Latham, private communication, 1971). The 
lunar module contained approximately ISO kg 
of volatile propellants, primarily dimethyl 
hydrazine fuel (CH 8 NHNHCH 8 ) and nitrogen 
tetroxide oxidizer (N 2 0,). There was an addi- 
tional source of energy if one considers the 
possibility of these hypergolic propellants com- 
bining. The heat of oxidation of dimethyl hy- 
drazine is 3.3 X 10 7 joules/kg [ Goodger , 1970], 
and hence, if all of the propellants were able to 
combine, the energy released would be ^3 X 
10® joules maximum or comparable to the im- 
pact energy. 

In Figure 2 are shown the counting rates of 
channel 6 of analyzer A, measuring positive 
ions with energies of 50 to 150 ev per unit 
charge, and of channel 3 of the same analyzer, 
measuring negative particles with energies of 
61 to 68 ev, for the period OOh 44m 53s to OOh 
48m 55s GMT on February 7, 1971. 

As can be seen from Figure 2, the counting 
rates before and during Antares impact were 
reasonably constant. The preimpact counting 
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rates in the low-energy negative-particle chan- 
nels were due to an ambient population of 
photoelectrons that was present whenever the 
lunar surface in the vicinity of CPLEE was 
illuminated. We note as proof of this assertion 
that these ambient fluxes disappeared entirely 
during the total lunar eclipse that occurred a 
few days later on February 10, 1971. The back- 
ground rate of the ion channel (channel 6) 
was due to various sources of contamination, 
including a small (~$%) contribution from 
electron scattering within the analyzer. (This 
effect was well documented in preflight calibra- 
tions.) The ion-channel background level was 
essentially constant throughout the lunar orbit 
and hence was not due to a low-energy tail of 
the solar-wind ion flux. 

Beginning at T + 4S sec, a series of pro- 
nounced enhancements above background levels 
in counting rates of both the ion and the nega- 
tive-particle channels was observed, with the 
data dominated by two major enhancements 
centered at T + 58 sec and T + 74 sec, respec- 
tively. Because the enhancements were observed 
simultaneously in particles of both charge types, 
we refer to these events as plasma clouds. 

The same data for analyzer B oriented 60° 
from vertical toward lunar west (i.e., toward 
the impact point) are shown in Figure 3. The 



Fig. 1. A sketch of the geometry of the impact event, showing the location of the impact 
point relative to the location of CPLEE and the Apollo 12 Side instrument [ Freeman et ai. } 
1971], Also shown are the incident solar -wind and assumed interplanetary magnetic-field 
directions. 
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UNf VE RSAL TIME ( DAYS/ HOU RS/MI N. /SEC ) 

Fig, 2, The counting rates of channel 3 and channel 6 of analyzer A at — 35 volts, meas- 
uring 65-ev negative particles and ions with 50 < E < 150 ev with peak response at 70 ev, 
respectively, showing the particle fluxes resulting from the lunar module impact. 


background level of the ion channel in analyzer 
B is considerably higher than that for analyzer 
A for reasons that remain unknown to us. The 
effect however was well documented in a series 
of postdeployment tests with N™ beta sources 
that were attached to the underside of a re- 
movable dust cover protecting the analyzer 
apertures during deployment and lunar-module 
ascent. As for the analyzer A data, however, 
we measured the flux enhancement as the dif- 
ference between the instantaneous counting 
rate and the background level. From compari- 
son of Figures 2 and 3, one can note that the 
flux enhancements were essentially simultaneous 
in two directions, but the ion flux measured by 
analyzer A was 5 times higher than the ion flux 
measured by analyzer B. The geometric factors 
of the corresponding sensors in analyzers A and 
B are essentially identical, and hence the rela- 
tive flux magnitudes can be directly compared 


by comparing the relative counting-rate en- 
hancements above the background levels of the 
channels. On the other hand, the negative-par- 
ticle flux measured by analyzer A was only 
one-third as great as the negative -particle flux 
measured by analyzer B, Examinations of Fig- 
ures 2 and 3 also show sporadic enhancements in 
particle fluxes occurring after the two initial 
plasma clouds. We consider that these enhance- 
ments were also the result of the impact but 
were due to various 'aftereffects/ perhaps sec- 
ondary impacts of ejecta material, (See dis- 
cussion below.) In this paper, therefore, we 
concentrate on the first- two well-defined en- 
hancements. 

The detailed characteristics of the two domi- 
nant plasma clouds are shown in Figure 4, a 
plot on an expanded time scale of the negative- 
particle fluxes in five energy ranges and ion 
flux in a single energy range measured by ana- 
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lyzer A. The plot shows clearly that the nega- 
tive-particle enhancement was confined to en- 
ergies less than 100 ev, since the 200-ev flux 
remained essentially constant throughout the 
event. The plot also shows that the enhance- 
ments of all the particles measured were simul- 
taneous to within the temporal resolution of 
the instrument (2.4 sec). 

The negative-particle spectrum is seen (Fig- 
ure 4) to vary throughout the event in both 
the magnitude of the fluxes and the shape of 
the spectrum, A comparison of the preimpact 
negative-particle spectrum and the spectrum 
during the enhancements is shown in Figure 
5. The first spectrum was measured at 00k 
42m 38s, or during the period of stable, ambient 
fluxes some 3 min prior to impact. The second 
spectrum was measured at 00k 46m 21s, or 
during the first enhancement. The differing 
spectral shapes are clearly seen in this figure, 

It might well be questioned whether the 
flux enhancements at T 7 -f- 58 and T -f- 74 sec 
were actually initiated by the Antares impact. 
Indeed, in the time period of approximately 2 


days following the impact event, when CPLEE 
was in the magnetosheath, several rapid en- 
hancements in the low-energy electron fluxes 
by up to a factor of 50 were observed. How- 
ever, these other enhancements were not cor- 
related with positive-ion flux increases, and, in 
fact, the event referred to here is the only such 
example of such perfectly correlated low-energy 
ion and negative-particle enhancements seen in 
this time period. In addition, careful monitoring 
prior to the impact revealed that the fluxes 
were relatively stable, constant to within a 
factor of 2 over time periods of a few minutes. 
These facts lend credence to the belief that 
this situation is a valid cause and effect one. 

Further confidence in our interpretation that 
the flux enhancements were artifically impact- 
produced rather than of natural origin is gained 
by noting that although no such plasma clouds 
were previously detected resulting from impact 
events, Freeman et ol . [1971] reported detec- 
tion of positive-ion clouds with the Apollo 12 
suprathermal ion detection experiment (Side), 
which they concluded resulted from the Apollo 
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Fig. 3. Same as Figure 2, except showing data from analyzer B. 
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Fig, 4. Ad expanded view of the data of Figure 2, showing details of the two prominent 
peaks. Fluxes computed from five negative -particle energy ranges and one ion energy range 
are shown. 



ENERGY (ev J ENERGY {ev ) 

Fig. 5. Electron spectra measured by analyzer A for two periods. The first is a few minutes 
prior to impact, and the second is the time at the height of the first large peak in Figure 2, 
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13 and 14 Saturn IV-B stage impacts. Further- 
more, the positive-ion components of the plasma 
clouds reported here were also detected by the 
Side at the Apollo 12 site, located 116 km west 
of the impact point (J. W. Freeman, Jr., pri- 
vate communication, 1971). 

It is concluded, therefore, that the impact 
of the ascent stage of the Apollo 14 lunar- 
module was responsible for the positive- and 
negative-particle fluxes observed by CPLEE, 
and these fluxes are referred to as plasma 
clouds. The salient features of the event are 
the time delay between the impact and the flux 
enhancements (~60 sec) and the simultaneous 
appearance of positive and negative particles. 

There are two possible interpretations 
of these data in a gross sense. It can be as- 
sumed that the particles were created and 
energized at the instant of impact or that the 
impact created an expanding neutral gas cloud 
and the components of the neutral cloud were 
ionized and accelerated by mechanisms that 
were more or less continuously active and in- 
dependent of the impact itself. 

If it is assumed that the particles were en- 
ergized at the instant and point of impact by 
some unknown mechanism, it is necessary to 
explain the subsequent behavior of the plasma 
clouds. 

According to this hypothesis, the plasma 
clouds had an average travel velocity of ^1 
km/sec and horizontal dimensions of 14 and 7 
km, respectively, for the first and second clouds. 
Noting that the positive and negative particles 
appeared simultaneously, a mechanism must be 
found to explain both the cloud containment 
and the relatively slow propagation velocity. 
It can be postulated that the positive-ion-di- 
rected velocity was of the order of the inferred 
plasma-cloud propagation velocity (^1 km/ 
sec), and then one can appeal to ambipolar 


diffusion to contain the negative particle com- 
ponent, if it is assumed that the negative par- 
ticles observed were electrons. Several calcu- 
lated parameters of 50-ev charged particles of 
various masses are listed in Table 1, and it is 
seen from this table that, to fit the foregoing 
hypothesis, the ion mass would have to be on 
the order of 1000 amu. Since the gas released 
at impact probably consists mainly of vaporized 
lunar-module propellants and lunar-surface ma- 
terials, we would estimate ion masses in the 
range 25-100 amu, but it is difficult to see how 
mass 1000 ions could have been created. In- 
deed, this assumption is substantiated by the 
observation of the Apollo 13 Saturn IV-B im- 
pact ion cloud by Freeman et ol. [1971] with 
the Apollo 12 Side instrument. The mass-ana- 
lyzer portion of the instrument showed peak 
ion fluxes in the range 66-90 amu/unit charge. 

Rejecting the hypothesis that the particles 
'traveled in straight-line paths between the im- 
pact point and CPLEE, there is still the pos- 
sibility that the particles could have been 
energized at the instant of impact and the 
trajectories influenced by a local magnetic field 
or that the plasma cloud could be magnetically 
confined. The measurements of the lunar-sur- 
face magnetic field by the Apollo 12 lunar 
surface magnetometer [ Dyal et al 1970] 
showed a steady field of 38 zt 3 y, with vector 
components B x = 26 y, B y — 13 y, and B, = 
—24 y. (In this coordinate system, x is south, 
y is east, and z is the direction of the local 
vertical.) Measurements of magnetic fields at 
two locations near the Apollo 14 site were made 
with the lunar portable magnetometer [ Dyal 
et al., 1971]. These sites (1 and 2) were located 
approximately 500 and 1700 meters from 
CPLEE, respectively. At site 1 the field was 
103 y with components B x =r 24 y, B y = 38 y, 
and Be = — 93 y, and at site 2 the field was 


TABLE 1. Plasma Cloud Parameters 


Cyclotron Radius, km 

Particle Velocity, Energy Density, : 

Energy, ev' Charge Sign Mass, amu km /sec ergs/cm® 36-y Field 100-y Field. 


50 

+ 

1 

100.0 

5.6 X J0"“ 

30 

10 

50 

+ 

25 

20.0 

28.0 X 10-“ 

150 

50 

50 

+ 

100 

10.0 

56 X 10-“ 

300 

100 

50 


1000 

1.0 

560 X 10-“ 

3000 

1000 

50 


m e 

4300 


0.7 

0.23 


1298 


Reasonek and O'Brien 


43 y with components £* = 19 y , B y = —36 y, 
and B c — — 15 y. The reader is referred to the 
above references for a complete description and 
discussion of the lunar-surface magnetic fields, 
but the data do suggest [Dyol et d 1971] 
that surface magnetic fields of several tens of 
Y*s exist over wide regions of the lunar surface. 
By contrast, magnetic-field measurements by 
the lunar-orbiting Explorer 35 spacecraft showed 
values of 10-12 y 800 km above the lunar sur- 
face IN ess et <d., 1967]. From these data we 
might postulate that the plasma clouds were 
magnetically confined in the enhanced mag- 
netic field close to the surface. However, when 
it is recalled that, according to the hypothesis, 
the dimensions of the two clouds were 14 and 
7 km and when it is argued that the cyclotron 
radii of the particles can be no larger than the 
cloud dimensions, it is seen from Table 1 that, 
the ions would have to be predominately of 
small masses (he., protons). We have argued 
above, however, that the ions probably have 
masses in the range 25-100 amu, and these ions 
would have cyclotron radii (see Tabic 1) too 
large by a factor of at least 5 to fit the ob- 
served data. 

Therefore it appears that it is impossible to 
reconcile the observed data with the hypothesis 
that the charged particles were energized at 
the instant of impact and then propagated in 
some manner to the location of CPLEE. The 
time delay between impact and observation by 
CPLEE and the relatively short duration of 
the enhancements were seen to require, depend- 
ing on the mode of propagation chosen, either 
extremely large (^1000 amu) or extremely 
small ('-'1 amu) ionic masses, and it was 
argued that such extreme values are highly 
unlikely. 

An alternate hypothesis is that the lunar- 
module impact produced expanding annular gas 
clouds and the components of the gas clouds 
were then ionized by solar photons or other 
mechanisms and subsequently energized by a 
continuously or erratically active acceleration 
mechanism. These fluxes were then observed by 
CPLEE only when the expanding annular gas 
cloud was in the vicinity of the instrument. 
Thus, according to this hypothesis, the velocity 
of 1 km/sec deduced from the impact-CPLEE 
distance and the delay time is a characteristic 
velocity of the gas-cloud expansion. The fact 


that there were two large enhancements, and 
by inference two gas clouds, can be explained 
by noting that the lunar-module impact tra- 
jectory was at a low (~10° ) elevation angle, 
and this could, of course, lead to secondary 
impacts following the primary impact. 

We can only speculate as to the mechanism 
responsible for energization of the charged par- 
ticles. We note that the solar magnet ospheric 
coordinates of CPLEE at the time of impact 
were Tsar ~ 34 Be and Zsv = 21 Bb and the solar 
elevation angle was 30°. Examination of the 
complete CPLEE data records before and after 
the impact show that the impact occurred just 
before the instrument crossing from the inter- 
planetary medium into the magnetosheath. 
Therefore the solar wind had direct access to 
the lunar surface at the time of the impact 
event. 

Manka and Michel [ 1970 ] calculated the trajec- 
tories of ions created near the lunar surface 
and accelerated by the V xB electric field of the 
solar wind. Although their calculated electric- 
field values (2-4 v/km) are certainly of sufficient 
magnitude to produce the observed particle 
energies, there are two observational features 
of these impact data that cause the hypothesis 
of acceleration in a static electric field to be 
rejected immediately. The first feature is that 
energetic particles of both charge signs appeared 
simultaneously, and the second feature is that 
positive ions resulting from the impact were 
detected both by CPLEE located east of the 
impact site and by the Apollo 12 Side located 
west of the impact site (J. W. Freeman, Jr., 
private communication, 1971 ). 

The solar-wind energy density is ~80 X 10~ w 
erg/cm 3 , and, by comparing this value with the 
range of plasma-cloud energy densities calcu- 
lated from the measured flux (see Table 1), it is 
seen that the solar wind is energetically capable 
of being the energy source. Whether or not 
interaction between the solar wind and a gas 
cloud can actually accelerate particles to the 
observed energies and fluxes is unknown, al- 
though Alfven [1954] and Lehnert [1970] 
pointed out that strong interactions may occur 
between magnetized plasmas and neutral gases. 

In summary, these lunar-module impact data 
indicate a situation of interaction among a neu- 
tral gas cloud, the solar wind, and possibly 
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local lunar magnetic fields, and thus a unique 
problem in plasma physics is presented. 
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The charged particle lunar environment experiment (CPLEE), a part of the Apollo 14 
lunar surface package, is an ion-electron spectrometer capable of measuring ions and electrons 
with energies between 40 ev and 50 kev. The instrument, with apertures 26 cm above the 
surface, lias detected a photoelectron gas layer above the sunlit lunar surface. No detectable 
flux above 200 ev has been observed. Experimental data for periods while the moon was in 
the earth's magnetotail for electrons with energies 40 ev < E < 200 ev follow a power-law 
spectrum j(E) zz j n {E/En)~ ti with 3.5 < & < 4. In the absence of photoelectrons with E > 200 
wc assume that the surface potential is at least 200 volts. The modulation of this potential 
in the presence of intense plasma-sheet fluxes has been observed. Also, a detailed history of 
the February 10, 1971, total lunar eclipse, to determine the source distribution of high-energy 
solar photons, is presented. A classical pcnumbral-umbral behavior indicates (hat at the time 
of the eclipse the emission of higher-energy photons was uniform over tile solar disc. Numerical 
solutions- for the variation of electron density and potenlial as functions of height above the 
lunar surface were obtained. The solar photon spectrum T(hp) f obtained from various experi- 
mental sources, and (he photoelectron yield function of [he surface materials, Y{hv), are two 
parameters’ of the solution. Energy spectra at the height of the measurements for various 
values of Y{hv) were computed until a fit to experimental data was obtained. Using a func- 
tional form Y{hv) = [Y 0 (hv - W)/(W/ 2)] for 6 ev g hv < 9 ev and Y(hv) - Y« for 
hv >9 ev, where the lunar-surface work function W was set at 6 ev, we calculated a value of 
Yo = 0.1 electrons/photon. The solution also showed that the photoelectron density falls by 
5 orders of magnitude within 10 meters of the surface, but the layer actually terminates several 
hundred meters above this height. A hydrostatic model of the photoelectron layer has also 
been developed. It is shown that the numerically calculated pressure, density, and potential 
can be approximated by solving the hydrostatic equations with an equation of slate P/n v * zz 
constant out to 200 cm from the surface. Beyond this height, the equation of stale shifts 
toward the isothermal case, P/n — constant. 

The general problems of photoelectron emis- 
sion by an isolated body in a vacuum and in a 
plasma have been the objects of several investi- 
gations. For example, Medved r 1965] has 
treated electron sheath formations about bodies 
of typical satellite dimensions. Guernsey and 
Fu [1970] have considered the properties of 
1^1 infinite, photocnjitting plate immersed in a 
dilute plasma. Grobman and Blank [1969J ob- 
tained expressions for the lunar surface poten- 
tial due to photoelectron emission while the 
moon is in the solar wind. E. Walbridge (un- 
published manuscript, 1970) developed a set of 
equations for obtaining the density of photo- 
electrons as well as the electrostatic potential 
as functions of height above the surface of the 
moon while the moon is in the solar wind. By 
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assuming a simplified form of the solar photon 
emission spectrum, he could provide analytic 
expressions for these quantities. 

In this paper, wc report on observations 
of stable photoelectron fluxes with energies 
between 40 and 200 ev by the Apollo 14 
charged-particle lunar environment experiment 
(CPLEE). These observations, made in the 
magnetotail under near-vacuum conditions, are 
compared with numerically calculated photo- 
emission spectra to determine the approxi- 
mate potential difference between ground and 
CPLEF/s apertures (26 cm). Numerically cal- 
culated density and' potential distributions, 
when compared with our measured values, help 
us estimate the photoelectron yield function of 
the dust, layer covering the moon. 

We have also developed a hydrostatic model 
for a photoelectron gas in equilibrium above 
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Fig. 1. Five-minute averaged counting rates for CPLEE, analyzer A, channel 1, at —35 
volts, measuring 40-ev electrons on February 8, 1971. After 0300 UT, counting rates fell from 
high magnetosheath to stable photoelectron levels. 


the surface of the moon. An equation of state 
P = constant n y is postulated, where y is a 
free parameter to be determined from the nu- 
merical analysis. 

The Instrument 

Complete descriptions of the CPLEE instru- 
ment have been given by O’Brien and Reasoner 
[1971] and Burke and Reasoner [1972]. The 
instrument contains two identical charged- 
particle analyzers, here referred to as analyzers 
A and B. Analyzer A looks toward the local 
lunar vertical, and analyzer B looks 60° from 
vertical toward lunar west. 

The particle analyzers contain a set of electro- 
static deflection plates to separate particles ac- 
cording to energy and charge type, and an array 
of 6-channel electron multipliers for particle 


detection. For a fixed voltage on the deflection 
plates, a five-band measurement of the spec- 
trum of particles of one charge sign and a 
single-band measurement of particles of the 
opposite charge sign are made. The deflection 
plate voltage is stepped through a sequence of 
3 voltages at both polarities, plus background 
and calibration levels with zero voltage on t^e 
plates. A complete measurement of the spec-* 
trum of ions and electrons with energies be- 
tween 40 ev and 50 kev is made every 19.2 sec. 
Of particular relevance to this study are the 
lowest electron energy passbands. With a de- 
flection voltage of —35 volts, the instrument 
measures electrons in five ranges centered at 40, 
,50, 65, 90, and 200 ev. With H-35 volts on the 
deflection plates, electrons in a single energy 
range between 50 and 150 ev are measured. At 
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the next higher deflection voltages, ±350 volts, 
the energy passbands given above are scaled 
upward by approximately a factor of 10. 

Observations 

In this section we present data from the Feb- 
ruary 1971 passage of the moon through the 
magnetotail. Because these are typical, display 
of the data from subsequent months would be 
redundant. At approximately 0300 UT on Feb- 
ruary 8, CPLEE passed from the dusk-side 
magnetosheat-h into the tail. The 5-min aver- 
aged counting rates for analyzer A, channel 1, 
at —35 volts measuring 40-ev electrons, are 
plotted for this day in Figure 1 . Almost identi- 
cal count rates are observed in analyzers A and 
B during this period of observation. As CPLEE 
moves across the magnetopause, the counting 


rate drops from ~ 200/cycle to the magnetotail 
photoelectron background of 35/cycle (1 cycle 
= 1.2 sec) . Enhancements at ~0530 hours and 
at ~0930 hours correspond to plasma events 
associated with substorms on earth \ Burke and 
Reasoner, 1972]. There is a data gap from 1000 
to 1200 hours. With the exception of the short- 
lived (<1 hour) enhancements, the detector 
shows a stable counting rate when the moon is 
in the magnetotail. 

Our contention is that these stable fluxes ob- 
served in the magnetotail during periods of low 
magnetic activity are photo electrons generated 
by ultraviolet radiation from the sun striking 
the surface of the moon. In support of this 
thesis, we have reproduced the counting rates 
observed in the same detector on February 10 
when the moon was near the center of the tail 



HOURS 

Fl |- 2 - Five-minute averaged deep-tail counting rates of 40-ev electrons on February 10. 1971. 
The lunar eclipse (0500 to 0900 UT) is marked by vanishing photoeleetron counting rates. 
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(Figure 2). First, we note that the stable count 
level is the same at the center as it- was when 
CPLEE first entered the tail. Secondly, from 
about 0500 to 1000 hours the moon was in 
eclipse. During this time, the counting rates go 
to zero. As the moon emerges from the earth's 
shadow, the counting rates return to their pre- 
eclipse levels. If the stable low-energy electrons 
were part of an ambient plasma, rather than 
photoelectrons, the counting rates would not be 
so radically altered as the moon moved across 
the earth's shadow. 

A detailed plot of the photoelectron flux on 
an expanded time scale for the period of lunar 
entry into the earth's penumbra and umbra is 


shown in Figure 3. The times of penumbral and 
umbral entry are indicated, and these times 
were computed from ephemeris data appropriate 
to the lunar coordinates of CPLEE. In addition 
to confirming our earlier arguments, this plot 
also shows that the high-energy (>40 ev) pho- 
tons are radiated essentially uniformly over the 
solar disc. If these photons were emitted from 
a few isolated regions, then one would observe 
sharp transitions in the flux as the regions were 
progressively shadowed by the limb of the 
earth. As can be seen, however, the curve in 
the penumbral portion is smooth, with no obvi- 
ous discontinuities. Furthermore, one notices 
that the counting rate falls by more than 2 



UNIVERSAL TIME 

Fig. 3. Expanded time-scale plot of photoelectron fluxes for the period of penumbral and 
umbral entry of the lunar surface region around CPLEE (analyzer A) during the lunar 
eclipse of February 10. The data are the counting rate of channel 6 at +35 volts measuring 
electrons of 50 ev < E < 150 ev. 
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Fig. 4. An example of ultraviolet contamination of analyzer A from —1200 of February 11 

to —0300 of February 12, 1971. 


orders of magnitude as the moon traverses the 
umbral and penumbra-1 regions (Figure 3). The 
minimum counting rate of ~ 15/sec for this 
channel (channel 6) is very near the instru- 
ment background level, but even if these were 
all due to an ambient plasma our argument that 
the pre-ccl ipse counting rates were due entirely 
to photoelectrons is not significantly altered. 

It could be argued that the observed count- 
ing rates were due to photons scattering within 
the detectors themselves and not due to ex- 
ternal photoelectrons. This, however, is not the 
case. Preflight calibrations with a laboratory 
ultraviolet source showed enhanced counting 
rates only when the angle between the look 
direction of the detector and the source was 
less than 10°. Given the 60° separation between 
the look directions of analyzers A and B } it 
would be impossible for the sun, essentially a 


point source, to produce identical counting rates 
in both analyzers simultaneously. There are 
Limes when we do observe ultraviolet contami- 
nation in one or the other channel. An example 
of such contamination is shown in Figure 4 
from 0600 hours,' February 11, to 0900 hours, 
February 12. As the sun moves across the 
aperture of analyzer A, the counting rates in- 
crease a full order of magnitude. During this 
period, analyzer B continued to produce typical 
deep tail counting rates. Note that as the de- 
tector came out of ultraviolet contamination it 
encountered typical magnetosheath plasma. At 
—0345, it passed back into the magneto tail, 
then at — 0800 returned to the magnetosheat.h. 

A typical spectrum of photoelectrons shown 
in Figure 5 was observed by analyzer A at 
—0400 hours on February 10, shortly before the 
moon entered penurnbral eclipse. The dark line 
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marks the differential flux equivalent to a back- 
ground count of one per cycle in each channel 
(channels 1-5). For all five channels, with the 
deflection plates at —35 volts, the differential 
flux is well above this background level. During 
geomagneticnlly quiet times, no statistically sig- 
nificant counts arc observed when the deflection 
plates are at —350 or —3500 volts, correspond- 
ing to electrons with E > 500 ev [Burke and 
Reasoner, 1972]. 

With the exception of periods of ultraviolet 
contamination in analyzer A, we always ob- 
serve nearly the same counting rale due to 
photo electrons in analyzers A and B. For all 
purposes, we can say that the spectrum dis- 
played in Figure 5 is just as typical for ana- 
lyzer B. We have found no case of anisotropy 
in the photoelectron fluxes. In all cases, too, we 
found that the photoeleetron spectra observed 
in both analyzers had dose to a power-law de- 
pendence on energy. If we write the differential 
flux in the form j(E) = jx is between 

3.5 and 4. In the next, section, the details of 
this spectrum are more carefully studied. 

Also in Figure 5 we display a schematic cross 
section of our instrument as it is deployed on 
the surface of the moon. The apertures of both 



Fig, 5. Typical photoeleetron spectrum ob- 
served by CPLEE at the lunar surface in the high- 
latitude magnetotail. 


analyzers are elevated 26 cm from ground. 
Their geometry is such that they observe only 
electrons with a component of velocity in the 
downward direction. Since wc continually ob- 
serve photoelcctrons with energies up to ^200 
ov, wc must assume that, the lunar surface po- 
tential is at least 200 volts during these times. 
This measurement will seem high to those 
familiar with the work of E. Walbridge (un- 
published manuscript, 1970) and Grobman and 
Blank [1969], who calculate a surface potential 
that is at least an order of magnitude lower. 
The difference is that their models deal with 
photocmissions from the surface of the moon 
in the presence of the solar wind. Our measure- 
ments in the magnetotail are made under near- 
vacuum conditions. After further analysis of 
the problem, we return to considerations of the 
surface potential. 

To summarize: During geornagnetically quiet 
times, wheu the moon is in the magneto tail and 
not in eclipse, stable photoeleetron fluxes with 
energies between 40 and 200 ev arc observed. 
These fluxes are isotropic and obey a power 
law, E~ M } where /a is between 3.5 and 4. From 
the fact that CPLEE is observing downward- 
moving electrons, we conclude that in the mag- 
netotail the lunar surface pot.cn lial is of the 
order of 200 volts. 

A Hydrostatic Model 

Our observation of steady photoeleetron 
fluxes that are isotropic over the lower half- 
plane for much of the moon's passage through 
the magnetotail suggests that wc can make the 
following assumptions about the physical situa- 
tion: 

1. The solar radiation flux at the lunar 
surface and the photoeleetron flux produced by 
the lunar surface are constant in time. 

2. There are equal probabilities for emitting 
photoelectrons into equal solid-angle elements 
in the upper half-plane. 

3. In the equilibrium situation, the net cur- 
rent out of the surface is zero. 

4. The lunar surface may be approximated 
by a flat infinite plane. Physical quantities vary 
only with height above the surface. (The co- 
ordinate system is such that X , Y y Z increase 
toward the local vertical, west, and south, 
respectively.) 
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The first assumption assures a constant flux 
of photoelcctrons at the surface of the moon. 
In the Appendix, we show that our postulate 
of an isotropic production of photoelectrons is 
sufficient to explain our observation of isotropic 
fluxes at. 26 cm. Assumption 3 demands that in 
a vacuum all photoelectrons be trapped by the 
electrostatic field above the surface. The surface 
potential under vacuum conditions is deter- 
mined by the most energetic photoelectrons. 

In the equilibrium situation, the governing 
equations are Poisson’s equation, 

= ~ iw <2 n ( x ) (i) 

and the equation of conservation of momentum, 


d<f>(x) = , _ j_ ap(») 
dx n(x ) dx 


( 2 ) 


Here the density and pressure are defined by 


n(x) = J f(v, x) d'v 

i‘(x) = | J mv 2 f(w, x) d’v 

The density can be eliminated from equations 
1 and 2 to show that 



= -8*7* 

J x dx 

= ±1 SnmX*)] 1 " 


dx f 


(5) 


In order to insure a potential that, decreases 
as x increases, we must choose the minus sign 
in equation 5. Equating the right-hand sides of 
equations 4 and 5, 


ya t ^ 2 dn 
q U dx 

= (87T «) ,/ V /1 

_ /„ \ i/2 


* = 2 (f) 

/ (y /2) —2 J 

1 n dn 

q \8t/ 

J n 0 


For y / 2, 

o 12/(7 -2) 

1 ~~ — y — X J (b ft ) 

and for y ^ 2 

n(x) — n n e _x/ * 
with the characteristic length 

X - (iV2 irqW) W2 (66) 

Equations 6a and 6& can be substituted into 
4 to get: 



4 > o 


However, to solve for all three quantities wc 
must postulate an equation of state 


(j>(x) = 


2P 0 

Xgn 0 



7-2 X y /<T ~ 2) 
7 X/ 


dx 


(7a) 


P(x) = P_o 
n y (x) 7i u y 


(3) 


where a is a constant,, and y is related to the 
poly trope index v by the relationship y = ( v 
+ l)/v. The values of physical quantities at the 
surface arc denoted by zero subscripts. In solv- 
ing equations 1-3, we demand that, as x — » oo, 
all physical quantities go to zero. 

Pressure can be eliminated by differentiating 
equation 3 and substituting into 2: 


7^2 


— 4>(x) - (1 — e 

qn 0 


= < 7 »> 

7=2 

For y 2 (equation 7a), there are two 
formal solutions, depending on whether y = 
1 or not: 


0<f> ya y- 2 dn 

- = _j — 

dx q dx 


(4) 


If we multiply (1) by d<j>/dx, using (4) on the 
right-hand side, and integrate in from infinity, 
we get 


<t>o ~ 4>(x) = 


Po 7 


qn 0 (y — I) 

(7 , 


y ^ i, 2 
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<t>o ~ <t>(x) = — — In (l + (7a") 

qno \ A/ 

7=1 

Th e solutions for the electrostatic potential, 
equations 7 7a" and 76, satisfy the boundary 
conditions = t£ 0 ~ 200 volts and <£(c o) = 
0. Here 0* = (A/qWn)y/(y - 1). If 1 < y 
< 2, <f>(x) goes smoothly to zero as x — *■ co. 
If y > 2, <£(:r) goes to zero at x ~ Ay/(y — 2). 
However, of 0 < y < 1, both boundary con- 
ditions cannot be satisfied for finite <£„. It would 
be possible to have an equation of state with 
finite <£ 0 and y < 1 near the surface, but the 
value of y must shift to a value greater than 1 
beyond some height. The numerical analysis of 
the following section shows that this is the case. 

Numerical Analysts 

General theory . The variations of photoelec- 
tron density and electrostatic potential above 
the surface of the moon can be calculated nu- 
merically. Again we approximate the lunar sur- 
face by an infinite plane, with the x direction 
normal to the surface, and assume spatial varia- 
tions of physical quantities only with the height. 

At a height x above tlic surface, the electron 
density is J /(v, x) d*v, where /(v, x) is the elec- 
tron distribution function. Tf we assume an 
isotropic flux at the surface, the Liouvillc 
theorem can be used to show that the distribu- 
tion function is independent of angles at all 
heights. Writing 

ctv = (2 E/my /2 dE dU 

and integrating over solid angles, we obtain for 
the density 

n(x) = 4?r [ (2E/rn) W2 i(E 7 x) dE (8) 

*f n 

Since the distribution function is a constant 
along particle trajectories, f(E, x) = f(E 0 , x ~ 
0), where E = E 0 — q[<j> 0 — By chang- 

ing the variable of integration from E to 
equation 8 can be expressed 

n(x) = 4tt 

■f {2m(E 0 - q \4> 0 -*(x)])} 1/2 

J 0[<£o-<Kx) I 


To calculate the distribution function of 
photoelect rons at the surface, wc consider the 
quantity 


j(E 0 ) dE(\ 

-[/; 


J{hv) Y(hv) p(E f)} hv) dhv 


] dE u 


( 10 ) 


the upward-moving flux of photoelectrons 
emitted from the surface with energies between 
E n and E 0 -f dE where I(hv)d(kv) is the 
flux of photons reaching the lunar surface with 
energies between hv and hv -j- d(hv). The 
quantum yield function Y (hv) gives the num- 
ber of electrons emitted by the surface per in- 
cident photon with energy hv, and p(E 0 , hv)dE 0 
is the probability that an electron emitted from 
the surface, due to a photon with energy hv, 
will have a kinetic energy between E { > and E 0 -f 
dE». The quantity p(E ih hv) is normalized so 
that 



1 


W is the work function of the lunar surface 
material. 


The total upward -moving flux at the surface 
is S« p (x — 0) = f<Tj(Eo)dE«. But 


S uv (x = 0) 




6, <t>, 0)v 0 2 dv 0 


• i sin 0 dO d<f> 

Since Vo — t? D [i cos 0 + j sin 0 cos <j> + k sin 0 
siii</>] and / is independent of angle, 

<5.„(x = 0) 

= »/ f(E 0 , X = 0) dE 0 (11) 

j o m 

Thus 


f(E 0> x = 0) = 


to 2 j(£ 0 ) 

2ttE 0 


( 12 ) 


7?,(x) = 2 


J Ql 


1 


{2 m(E a — q\4>„ — <^(x)]) | 1/2 




0 


• /C&'o, x — 0) dE, 


(9) 


( 13 ) 
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The potential as a function of height is 
evaluated by multiplying Poisson's equation, 
= — ^qn{x) y by d$/dx and integrating 
in from x — oo to get 

(d<t>/dx) 2 ~ — 8ir q f - n(<h') <J4>' (14) 


where we have written 


1. Adopted a work function of lunar ma- 
terial of 6 ev. 

2. Assumed a photoelectron yield function 
of the form 


Y(h>) - 0 
Y{hv) = Y 0 


hv — 6 


hv < 6 ev 
6 < hv < 9 ev 


r»o o * 0 

I n(x f )(dip/dx') dx f — j n(<f>') cW 


A further integration out from the surface 
gives us the potential at a point x . 

Computational methods and results. To de- 
termine the upward-moving differential flux at 
the surface, upon which the distribution func- 
tion, number density, and potential depend, wc 
must first calculate the integral in equation 
10. The solar photon differential flux at 1 AU, 
I(hv), is taken from Friedman [1963] for the 
range 2000 to 1800 A and from Hinteregger 
et al. [1965] for the range 1775 to 1 A ixnd, is 
plotted in Figure 6. Following the suggestion 
of Wat, bridge [1970], we have: 


Y(hv) — Y 0 hv > 9 ev 

(15) 

where Y 0 is a free parameter of our calculation. 

3. Chosen a probability function 

p(E, hv) = 6 E(E t - E)/E, 0 < E < E y 

( 16 ) 

p{E, hv) = 0 E > E, 

where 

E, = hv - W hv > W 

Ei = 0 hv < W 

In general, the probability function is a 
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complicated function depending on the nature 
of the photoemission material. However, Grob- 
man and Blank [1969] have shown that for the 
purpose of calculating equation 10 any broad 
function with zeros at 7? = 0 and E — E x and 
a width AE ~ hv will suffice. A plot of p(E : 
hv) is shown in Figure 7 for various values of 
E 

The upward-directed differential flux in elec- 
trons/cm 2 sec ster ev for the values Y 0 = 1, 
0.1, 0.01 were numerically computed and have 
been plotted in Figure 8. We have also inserted 
the photoelectron differential flux observed by 
CPLEE at 26 cm. The Liouville theorem allows 
us to set a lower bound on F 0 of 0.1. That is, 
if there were no potential difference between 
the ground and 26 cm, the yield function 
would be 0.1 eiect rons/photon. After estimating 
the potential difference between 26 cm and 
ground, we can also determine an upper bound 
on F ft . 

Solving the integro-differentinl equation 14 
for <f)(x) involves an integration from the sur- 
face outward, with an assumed value of ^> 0< 
However, the expression for d<f>/dx involves an 
integral from infinity in to x } or equivalently 
from </> — 0 to <j>(x). By the expedient of divid- 
ing the integral into pieces in E 0 space and 
using an analytic approximation to the func- 
tion in cadi of these intervals, a solution 


was effected. In this way, it was only necessary 
to know the values of <j> — <j>(x) and <j> = 0 at 
the end points of the interval, and the solution 
would proceed. In Figure 9 we show families of 
solutions for $(x ) with several values of the 
parameter F 0 . 

The value of F ft calculated by assuming no 
potential difference between the surface and x 

— 26 cm was 0.1. Figure 9 shows that, for F 0 

— 0.1, the potential difference <j>(x = 0) — 
4>(x — 26 cm) is only 3 volts. Obviously, we 
could now use an iterative procedure, modifying 
our spectral measurement at 26 cm to obtain 
the surface spectrum according to the equation 
j{E } s) j(E 0) 0) and hence obtain a new 
estimate of F 0 . However, the procedure is 
hardly justified considering the small potential 
difference (~3 volts) and the energy range of 
the measured photoelectrons (40 to 200 cv). 
Hence from our measured photoelectron fluxes 
and numerical analysis we obtain a lunar sur- 
face potential of at least 200 volts and a value 
of the average photoelectron yield of Ft, — 0.1 
electrons/photon . 

Lunar Surface Potential 0 o 

The presence of photo electrons at 200 ev 
without significant fluxes in the next highest 
energy channel at 500 ev leads us to conclude 
that the lunar surface potential is at least 200 



Fig. Probability function that a photon of energy hv will cause the lunar surface material 
to emit a photoelectron of energy E , with different values of 2?i z= hv — W. 
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Fig, 8- Numerically computed phofoeleetron 
spectra emitted for the yield functions Yo =; 1, 
0.1, and 0.01 electrons per photon. The photo- 
electron spectrum measured by CFLEE (shown 
by plus signs) is found to fall close to the To = 
0.1 line. 

volts. The data of Hinteregger et ai. [1965] 
show significant solar photon fluxes up to 400 
ev, and presumably the lunar surface potential 
under vacuum conditions could be 400 volts. 
However, we detected no photoelect rons with 
E > 400 ev, although it should be noted that 
the extrapolation of our measured spectrum 
(Figure 5) to 400 ev is below the instrument 
background. We have therefore adopted a con- 
servative value of 200 volts as the lunar surface 
potential for the calculations in the preceding 
and following sections. 

The lunar surface potential can be decreased, 
however, by the presence of a hot. ambient 
plasma that furnishes an electron return cur- 
rent that partially balances the emitted photo- 
electron current. In effect, the highest-energy 
photoelectrons can escape from the potential 
well, since electrons from the ambient plasma 
furnish the return current to balance these 


escaping photo electrons. Quantitatively, if F t 
is the net negative flux to the lunar surface 
from the ambient plasma, and is the 

emitted photoeleetron energy spectrum in units 
of electrons /cm 2 sec ev, then 



and this equation ran be solved for ^ (l , the 
lunar surface potential. 

Our measurements of phofoelectrons were 
taken during periods in the magnet of ail when 
all the channels of the instrument except the 
lo west-energy electron channels were at back- 
ground levels. Thus we can establish an upper 
limit to the electron flux from the ambient 
plasma for electrons with 40 ev < E < 50 kev. 
Figure 5 shows the 'background spectrum/ 
calculated by converting the background count- 
ing rate of count/sec to equivalent, flux in 
each of the energy channels. Integrating over 
this spectrum and converting to flux over the 
hemisphere gives F * < 3.4 X 10 E elect r on s/ertr 
sec. We feel that this is a valid upper limit, 
as the range of measurement in energy includes 
the peak energy of the plasma-sheet spectrum 
(^1 kev) . 

Wo note that Vasyliunas [1908] obtained an 
upper limit to the electron concentration for 



HEIGHT (CM) 

Fig. 9. Numerically computed potential dis- 
tribution above the lunar surface for several val- 
ues of the yield function Fo. For Y n = 0.1, the 
potential difference between ground and 26 cm is 
about 3 volts. 
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locations outside the plasma sheet based on 
Ogo 3 data. The relation expressed was NE^ rz < 
10~ 2 cm -3 kev ,/a , where N is the electron density, 
and Ed is the energy at the peak of the spectrum. 
For an isotropic plasma in which the hulk 
motion can be neglected relative to the thermal 
motion, the electron flux to a probe is given by 
F 9 = Nv/2(tt) u 2 . Applying the appropriate con- 
version factors, the expression of Vasyliunas 
results in an upper limit to the electron flux of 
E* < 5.6 X TO 7 clcctrons/cm 2 sec. 

The emitted photo elec Iron energy spectrum 
j(E u ) is shown in Figure 8. The procedure 
involved in calculating I he surface potential <j> a 
is to integrate the function j(E) from the maxi- 
mum energy of 200 ev backward until the 
total flux is equal to the upper limit of the 
return flux. The computation was done for V Q 
— 1 , 0.1, and 0.01, and for the two values of 
the upper limit, of the return flux derived above. 
The results are shown in Table 1. 

The lower half -height of the channel 5 energy 
passband is 160 ev. Hence the surface potential 
could be as low ns 160 volts and still result in 
particle fluxes in channel 5. This estimate of 
the potential is seen to be not inconsisicnt with 
values of F n — 0.1 and F 8 < 3A x 10°, result- 
ing in a surface potential (Table 1) of 114 volts. 

One rather obvious prediction of our argu- 
ments about the surface potential is that, when 
the electron flux reaching the surface of the 
moon is high enough, surface-generated photo- 
electrons with energies in the range of our de- 
tector should vanish. This happens in the solar 
wind and magnetosheath, but it is impossible 
for CPLEE to provide conclusive observational 
evidence in the presence of contamination by 
solar wind and magnetosheath electrons that 
the > 40-ev photoelectrons are not returning to 
the surface and hence entering the detectors. 
Electron densities >1 cm" 3 and temperatures 


TABLE 1. Lunar Surface Potential 


Electron Flux 

To 

volts 

3.4 X 10 6 

1.0 

181 

3.4 X 10 6 

0.1 

114 

3.4 X 10« 

0.01 

44 

5.6 X 10 7 

1.0 

96 

5.6 X 10 7 

0.1 

36 

5.6 X 10 7 

0.01 

8 


^15 ev, such as those commonly encountered in 
the solar wind and magnet os heath, provide 
much higher fluxes in the 40- to 200- ev range 
than photoelectron fluxes observed in the high- 
latitude magnet otail. 

On April 9, 1971, a worldwide magnetic storm 
was observed by CPLEE \ Burke et at., 1972] 
in which the magnetosheath moved in to F,™ — 
15 R e . As the magnetopause moved out past the 
moon, intense plasma-sheet fluxes were ob- 
served. At this time, the fluxes observed in the 
40-, 50-, and 70-cv channels fell below photo- 
electron levels. In Figure 10(a) we show the 
count rates for the 40-ev and 500-ev electron 
channels from 1100 to 1200 UT on April 9. A 
heavy line has been drawn at the 40-ev photo- 
electron level. As the 500-ev count rate rises, 
the 40-ev count rale falls below this line. From 
1128 to 1135 UT, when the 500-ev count rate 
dropped, the 40-ev channel returned to the 
photoelectron level. In Figure 10(5), we have 
plotted the electron spectra at 1110 (plasma 
sheet) and at 1135 (photoelectron). The inte- 
gral flux at 1 1 10 was calculated to be 2.5 X 
10 s cm" 2 sec" 1 . Such a flux reduces the surface 
potential to less than 10 volts. Thus we are 
observing a modulation of the surface potential 
by plasma -sheet fluxes into the lunar surface. 

Discussion 

For the sake of comparison with the predic- 
tions of the hydrostatic model, we have plotted 
the numerically calculated density, pressure, 
and potential difference from x — 0 out to a 
height of 200 meters in Figure 11. In Figure 
12 the pressure is plotted as a function of num- 
ber density to obtain the equation of state. We 
find that, from 0 < x < 30 cm, the value of y 
is about 0.5. From 30 to 100 cm, y drops to a 
value of 0.2, then recovers to about 0.5 out to 
x — 2000 cm. Beyond this point, y shifts toward 
a value greater than 1. 

The dashed lines in Figure 11 represent the 
density and potential difference as computed 
from equations 6a and 7fl' using y = 0.5. A 
surface pressure of ~2 X 10“* erg/cm" and 
density ^6 x 10 4 cm" 3 give a value of \ = 
2.1 cm. Out to x = 200 cm, the numerically 
calculated potential agrees quite well with the 
hydrostatic prediction. Beyond this height, the 
potential difference rises less steeply than the 
y — 0.5 prediction. However, this can be under- 
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Fig. 10. (a) Counting rates for the 40-ev and 500-ev channels between 1100 and 1200 UT, 
April 9, 1971. <f>) Electron spectra observed at 1110 and 1135 UT, showing plasma-shcct and 
photoele citron fluxes. 


stood in terms of the shift to larger than unity 
values of 7 required by equation 7a' if the 
boundary condition <£(oc) = 0 is to be met. 

The density curve is much more sensitive to 
fluctuations in the value of 7. The variations in 
the region 30 < x < 300 cm correspond to 
potential differences of 3 to 15 volts from 
ground. We note that in Figure 8 the photo- 
electron flux generated at the surface has sharp 
breaks in this region. Evidently the photoclec- 
trons can be broken up into three groups of 
low (0 < E < 1 ev), medium (1 < E < 10 
ev), and high (10 < E < 200 ev) energy. 
Where one distribution dominates over the 
others, a value of 7 is established. Fluctuations 
in y are found in the transition regions between 
populations. 

Feuerbacher et ai . [1972] have measured the 
photoclcctron yield of a lunar fine sample. In 
the photon energy range 5 < E < 20 ev, they 
found a yield function that reaches a maximum 
value of -—'0.08 at 15 ev, then drops to 0.01 



Fig- 11. Numerically computed values of elec- 
tron density, potential difference, and pressure as 
a function of distance from the lunar surface. The 
dashed line represents a 0.5 power-law curve. 


p( X J (ergs/cm 3 ) 
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Fig. 12. Plot of photoclectron pressure against density lo determine the local equation of 

state. 


at 20 ev. B. Fitton (private communication, 
1071) has suggested that our value of 0.1 is 
more a measurement of the CP LEE instrument 
case yield function than that of the lunar sur- 
face. We find it difficult to understand how 
(his could be the case. 

First, were OP LEE an electrically isolated 
package, the requirement, that the net current 
to the instrument be zero would result in a 
measured yield function that is representative 
of the case material. Relative to the lunar sur- 
face, tiie case would bear a positive potential in 
order to maintain an enhanced photoclcctron 
density in its immediate vicinity. 

The fact is, however, that CPLEE is not 
electrically isolated. CPLEE is connected to the 
central Alsep station. Further, at any given 
time, only the top surfaces and one or two sides 
of Alsep packages are illuminated by the sun, 
while the remaining area is shadowed. These 
unilluminated surfaces provide receptor areas 
for return current from the photoelectron gas. 
Thus no large potential difference can develop 
between CPLEE and the lunar surface, and 
phot oelectrons emitted at the lunar surface 
and at the CPLEE case are indistinguishable. 
Geometrical considerations of electron trajec- 


tories would lead us to expect that the bulk of 
the photoelectrons measured by CPLEE were 
emitted at the lunar surface at least several 
meters from CPLEE. if the Alsep instrument 
cases had a photoelectron yield much larger or 
much smaller than the yield of the lunar sur- 
face, then one would expect a perturbation of 
the photoelectron flux in the vicinity of Alsep. 
This perturbation would depend not only on the 
ratio of photoolectron yield, but, more impor- 
tant, on the ratio of the area of the Alsep in- 
struments to the area encompassed by the tra- 
jectories of E > 40 ev phot oelectrons from a 
point source. This last, area is of the order of 
the square of the scale height of photoelectrons 
with E > 40 ev (~10 meters. Figure 11). Since 
the Alsep area is nr, the ratio of areas is 
of the order of 2%. Thus, even if the yield of 
the instrument cases were a factor of 10 greater 
or smaller than the surface yield, the flux per- 
turbation would only be 20%. 

Second, were CPLEE measuring its own 
photoelectrons, one would expect to observe 
changes in the relative fluxes observed in the 
two analyzers with solar zenith angle as the 
electron cloud surrounding the instrument ad- 
justs to changing ilium inn lion. Specifically, the 
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ratio of the flux observed by analyzer B (look- 
ing 60° west of vertical) to the flux observed in 
analyzer A (looking to the vertical) should be 
larger after lunar noon than before. Our data 
show isotropic photoelectron fluxes across the 
entire magnetotaii. 

Summary and Conclusions 

In this paper, we have reported the obser- 
vation of stable, isotropic photoelectron fluxes 
26 cm above the lunar surface. In the energy 
range 40 < E < 200, the flux obeys a power 
law of the form j{E) = jo(#/En) _/i , where \l 
is between 3.5 and 4. Because these fluxes were 
moving down, we conclude that, in the near- 
vacuum conditions of the high-latitude inag- 
nctotail, the lunar surface potential is at least 
200 volts. The modulation of the surface po- 
tential in the presence of intense plasma-sheet 
fluxes has also been observed. It was shown 
that these electrons can be explained in terms 
of the measured solar photon spectrum produc- 
ing an isotropic flux of phot del eetrons at the 
surface. A photo electron yield function of 1 d 
~ 0.1 elec tron/pho ton was calculated. Finally, 
wc have shown that the numerically calculated 
pressure, density, and potential distributions 
can be approximated by the solutions to a set of 
hydrostatic equations that employ an equation 
of state P/n' n — constant out to 200 cm from 
the surface. Beyond this height, the equation 
of state shifts toward the isothermal case, P/n 
— constant. 

ArrENDix 


P(x) — j mvv/(v 7 x) cfv 

= J vvP(E, ft, X) dE d£l 

Here we have used d'v = ^dv sin t) dO -- 
{2E/m'Y'’dE da and have defined the direc- 
tional differential density, flux, and pressure 

n(e, a, x) = (5) , 2 ^ e - x) 

S(E, a, x) ^ -k ]{E, a, x) 042) 

P(E, ft, TC) ' (“)"’»’, <!. « 

The angular, dependence of these quantities is 
contained only in the distribution function. To 
compare one of the directional differential 
quantities at a point X L with its value at the 
ground, X = 0, we use the Liouville theorem: 

h(Ei, Qi, *i) = /o(«o, *o) 

The distribution function is a constant along 
particle trajectories. Subscripts 0 and 1 indi- 
cate the value of the quantity at x = 0 and an 
arbitrary height x = X*. 

In our model, we assumed an isotropic photo- 
electron production (at X = 0) over the upper 
half-plane in velocity space. Since all electrons 
arc trapped in a conservative field, isotropy is 
maintained over the whole of velocity space at 
X ~ 0. Thus U is independent of G 0 . From the 
conservation of energy 


Here we present a justification for using a 
scalar pressure in the equation of state (equa- 
tion 3). 

By definition, the number density, flux, and 
pressure are 


n(x) = J }(v, x) d'v 

= j N(E, ft, X) dE i 
S(x) = j v/(v, x) d*v 

--- J 0 S(K, ft, X) dE 


da 


(A 1) 


£?„ = £1 + q[4> f> ~ <*>(*,)] 
or 


h{Ei, «i, V.) 

= /o(Ei + <?[<£„ - <K*i)] U 3 ) 

Thus if photoelectrons are isotropic at the 
ground, they are isotropic at X lf This explains 
our observation of isotropic fluxes, measured at 
X — 26 cm. In this case, the directional dif- 
ferential pressure is also isotropic, and on inte- 
gration reduces to the scalar form used in the 
text. 

Being independent of angle, the distribution 
function Is an even function about v x = 0, v v 
— 0, and v r — 0. The Vlasov equation 
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§1 i 6/ _ e 60 6/ _ n 

cU x dx tn dx ft;* 

iias nontrivial moment* solutions only when 
multiplied by v M v 9 * t or v j; v*. These mo- 

ment equations take the form 

i~n(*)(V\x)) + ^f x n(x) = 0 (A4) 
and 

n(v*(x)v„ 2 (x)) 


and (44') becomes 

Sn(x)w 2 (x) 2e d<j> 


dx 


+ n(x) — 0 
m dx 


and (A5'), 


dn(x)w\x) 2 e 30 2 

JZ n(x)w = 0 


m vx 


Expanding, we get 


. 2/^1 dn(x)w\ x) 2ed<f> 


w 2 {x) 


^ -f — ^ n ( x \ 
dx ^ m dx { \ 




+ y(x)w 2 (x) 


dw 2 (x) 

dx 


n(v*(x)v*(x)) ( 

+ ~ 7 ^ «(*){fx 2 (*)) — 0 
well ' w/ 

t; n(oZ(x)v\x)) 


+ -~n(x)(, , x \x)) = 0 
m dx v x w/ 

where 

n(x){t)i(x)) - J x) 

Since the distribution function depends only on 
velocity, equations ( A4 ) and (.45) can be 
written 


at , ‘ w + »to ,,( * ) "° (/H,) 
#.«*) + $ !!«*>-» mm 


where 

h.(*) = f «*>, *) do 
Jo 

If, for example, we had a locally Maxwellian 
gas 


then 


f(v, x) 


— W ( J ) -»*/»’ ( 1 ) 

3/2 3 / \ # 

7T lt> (x) 


T / rj \ _ (2i> — l )H ra(s)tt> ai ' 2 (x) 

2K '‘ ^ eyv i l 


Since the bracketed term is zero and r?. (x)tc> 2 (a:) 
> 0, to*(x) is a constant. This is the isothermal 
case, whose solution n(x) — m 0 exp {—gl^o 
— <j>(x)]/mvf} is well known. 

Tn general, however, (A 4') and (45') cannot 
be solved without assuming a distribution func- 
tion, from which the equation of state can be 
determined from (A5') . 
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A geomagnetic storm beginning with a sc occurred on April 9, 1971. During the storm 
the charged particle lunar environment experiment at the Apollo 14 site, the solar wind 
spectrometer experiment at the Apollo 12 site, anti the Ames magnetometers on Explorer 
35 took data in the magnetosheath, at the magnetopause, in the plasma sheet, and in the 
high-latitude geomagnetic tail. The MIT Faraday cup and Ames magnetometers on-board 
Explorer 33 monitored the solar wind. The data show that the storm was caused by a 
corotating tangential discontinuity in the solar wind, the magnetopause position is strongly 
dependent on the attack angle of the solar wind, and the tail field strength was indirectly 
measured to increase from 10 to 14 7 after the sc. During the main phase the field strength 
in the tail was observed to increase to between 2S and 34 7. This increase is consistent with 
a thermal and magnetic compression of the tad radius from ~26 to ~16 Re. The data also 
show that plasma sheet electrons were observed to drift into the magnetosheath, the thick- 
ness of the magnetopause was measured by spatially separated instruments to be >3000' km, 
and the plasma and magnetic field data are consistent with a model for currents flowing 
along magnetic field lines in the boundary. 

In this paper we investigate the effects of 
the disturbed solar wind on' the topology of the 
magnetotail by correlating solar wind data with 
data taken near the magnetopause at lunar 
distance during the geomagnetic storm of April 
9, 1971. The topology of the magnetotail during 
geomagnetic storms has been the object of 
several studies. Sugiura et gl. [1968] have 
traced the propagation of a sc to Ogo 3 at a 
radial distance of ^10 R E behind the earth. 

They argue that the increase of field strength in 
the tail near the time of the sc is due to magnetic 
flux being transferred from the dayside instead 
of being due to a thermal compression of the 
tail by the solar wind. In the closed field line 
region of the near-earth tail the magnetic field 
strength has been observed to undergo diamag- 

Copyright © 19?3 by the American Geophysical Uflion. 


netic decreases due to the presence of the ring 
current during the main phase of storms [ An- 
derson and Ness, 1966] . However, in the deep 
tail outside the plasma sheet the field strength 
increases [Behannon, 1968]. As much as half 
of this high-latitude field increase has been 
attributed to the transport of dayside flux into 
the tail [Behannon and Ness, 1966]. 

The positions of the magnetopause during 
quiet and disturbed times have been studied 
by Behannon [1968, 1970], Fairfield [1971], 
Mihalov et al. [1970], and Ho'we and Binsack 
[1972] . Behannon^ data show that the magneto- 
pause at lunar distance has been observed to 
move in from quiet time positions of Y se — 
—21 and +25 R B to Y SK = — 15 and 21 R E on 
the dawn and dusk sides, respectively, during 
disturbed periods. 

In this paper we report on simultaneous ob- 
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servations made by several instruments of the 
magnetic storm of April 9, 1971. These obser- 
vations include (1) proton and electron 
flux measurements by the Bice University 
charged particle lunar environment experiment 
(CPLEE) at the Apollo 14 site, (2) proton 
fluxes by the Jet Propulsion Laboratory solar 
wind spectrometer experiment (SWSE) at the 
Apollo 12 site, (3) solar wind proton fluxes by 
the MIT Faraday cup experiment aboard Ex- 
plorer 33, and (4) magnetic measurements by 
the Ames Research Center magnetometers 
aboard Explorer 33 and Explorer 35. 

Data from these instruments are used to 
reconstruct a detailed history of the storm. 
Ground-based magnetometer records are pre- 
sented first for purposes of identifying the 
general temporal behavior of the storm. Solar 
wind particle and magnetic field characteristics 
are next displayed. The response of the geo- 
magnetic tail at the lunar distance is displayed 
by using data from CPLEE, SWSE, and the 
Explorer 35 magnetometer. The properties of 
the boundary layer region and an intense 
plasma sheet near that boundary are discussed. 

INSTRUMENTS 

The CPLEE, deployed at the Apollo 14 land- 
ing site, has boon described by O'Brien and 
Reasoner [1971]. The instrument contains two 
charged particle analyzers similar to an 
ion-electron spectrometer node-named Specs 
[O'Brien et al . , 1967]. The look direction of one 
analyzer is toward the lunar vertical, whereas 
the other is 60° from the lunar vertical toward 
lunar west. The normal voltage stepping se- 
quence requires 19.2 sec to complete a 15-point, 
spectrum over the energy range 40 ev to 20 kev 
for electrons and positive ions. 

The SWSE was placed on the surface of the 
moon as part of the Apollo 12 scientific package 
and has been described by Snyder et al. j 1970 ] * 
This instrument consists of seven Faraday cups, 
one of which faces the lunar vertical. The other 
six are arranged symmetrically about it but 
looking 00° away from the vertical. By mea- 
suring the relative and absolute currents in each 
detector the intensity and directions of solar 
wind and mngncloslienth fluxes reaching the 
surface of the moon are calculated every 28.1 
sec. Because of their high sensitivity thresholds 


the Faraday cups do not respond to magneto- 
spheric plasma. Thus SWSE data used in this 
paper are from magnetosheath and boundary 
layer protons. 

The MIT Faraday cup aboard Explorer 33 
has been described by Lyon et al. [1968]. This 
instrument, sensitive to electrons and ions with 
energies between 100 ev and 4 kev, determines 
the density, bulk, and thermal speed and the 
direction of plasma flow every 328 sec. 

The Ames magnetometers aboard Explorer 
33 and Explorer 35 have been described by 
Mikola v et al. [1968] and consist of three 
orthogonal flux-gate magnetometers. One vector 
measurement is made every 6.14 sec; in this 
paper, 81. 8-sec averages are presented. During 
periods when Explorer 35 was in the shadow 
of the moon and/or occulted by the moon with 
respect to the earth, no data were received. 

The locations of the various instruments in 
relation to the quiet time magnetosphere are 
shown in Figure 1. Explorer 33, monitoring 
the solar wind, is near (Z^, Y t0 , Z x „) — (5, 
—34, 33) Rt>. On April 9 at 1200 UT the moon 
was at {Y» c , Z^) — (15.9, — 3.8) R « and was 
well within the quiet time magneto! ail On the 
lunar surface, SWSE and CPLEE are separated 
by 182 km. The maximum separation between 
the lunar-based instruments find Explorer 35 is 
6 7? ril (l R m = 1738 km). 

Observations 

A selection of magnetograms for the first 16 
hours of April 9, 1971, is shown in Figure 2. 
For at least 3 hours before l ho sc at 0428 the 
horizontal components of magnetic field were 
at quiet time levels for all stations. The sc was 
observed most prominently at low-latitude sta- 
tions located in the afternoon sector. The initial 
phase of the storm lasted between 1.5 hours 
at Guam and 2.5 hours at Fredericksburg. By 
0630 the ring current had established the main 
phase depression at the low-latitude stations. 
The recovery phase began at ^1000 UT. Re- 
maining features of the magnetograms, particu- 
larly those of high-latitude stations, are pre- 
sumably from local effects such as subs form 
electrojets. We need not, dwell on these effects, 
but. we have presented the traces for the sake 
of completeness, 
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Solar Wind 

A cursory glance at the Kp indices for the 
first 6 months of 1971 (Solar Geophysical Data, 
1971.) indicates that the storm of April 9, 1971, 
was caused by a fast-moving stream co rotating 
in the solar wind. The February and April 
passages of this structure have been observed 
by the Explorer 33 plasma detector [ Howe et 
al.y 1971]. 

The solar wind proton and magnetic field 
data observed by Explorer 33 are displayed in 
Figure 3. Here N (per cubic centimeter) is the 
proton density, V (kilometers per second) is the 
bulk speed, <t> is the solar wind longitude, W is 
the most probable thermal speed, B is the inter- 
planetary field strength, and <p and 0 are the 
longitudinal and latitudinal field directions in 
the solar equatorial coordinate system. The 
angle is positive if the flow is from the east 
of the sun. Aberration due to the earth's orbital 
motion about the sun has been removed from 
the data. 

For our analysis of the effects of the storm 
time solar wind on the magnetotail we note: 

1 . For 4 hoiirs before the sc, N, V, W, and 
B were near steady values of 14 cm - ", 350 1cm/ 
sec, 35 krn/sec, and 4 y. 

2. At ~0435 the density, flow speed, and' 
magnetic field suddenly increased to ~30 cm -3 , 
400 km/sec, and 67, respectively, and the ther- 
mal speed decreased to 25 km/sec. This event; 
with a field and density increase and temperature 
decrease, is similar to a (+ , +, — ) tangential 
discontinuity reported by Burlaga [1968], If 
indeed the 0435 event was a tangential dis- 
continuity, the normal to the surface of discon- 
tinuity is given by Bi xB 2 /|B 1 xB 2 |. In R, T, 
and N coordinates [Turner and Siscoe i 1971] the 
preevent magnetic field was B j = (—1-5, 1.6, 
— 3.9) 7 and changed to B 3 == (—3, 3-3, —5.3) y. 
Thus a surface normal n = (0.76, 0.65, —0.3), 
which is illustrated in Figure 1, is given. Such 
a normal is consistent with the fact that the sc 
was observed only in the afternoon sector and 
that the leading edge of the fast-moving stream 
arrived at the Explorer 33 satellite 7 min after 
the sc. 

3. At ^0542 the solar wind direction shifted 
by 5° further from east of the sun and the 
interplanetary field increased to ~18 7. 

4. Between 0945 and 1000 the density de- 
creased to ^15 cm -s , the thermal speed ih- 



Fig. 1. Position of instruments with respect 
to the earth, the quiet time magnetopause, and 
the bow shock just prior to the sc event of April 
9, 1971, at 0426 GMT. 

creased to ^75 km/sec, and the flow shifted to 
a direction from east of the sun. The change in 
flow direction from east to west as the flow 
velocity increased from 375 to 600 krn suggests 
that the event at 1000 UT is a second tangential 
discontinuity within the solar wind structure. If 
it is so, the i?, AT, and T normal to the discon- 
tinuity is ft ~ (0.82, 0, 0.55). 

5. Between 1230 and 1240 the flow direction 
was again from east of the sun. After this the 
flow was from west of the sun, and the speed 
increased to over 600 km/sec. 

Magnetotail 

The responses of the magneto tail to the 
change in solar wind conditions as measured 
by CPLEE, SWSE, and the Explorer 35 mag' 
netometer are shown in Figures 4, 5, and 6 for 
hours 0300 to 1300. In Figure 4 the differential 
fluxes of electrons with energies near 50 and 600 
ev and protons of energy 1 kev are plotted. 
The density, flow speed, and thermal speed of 
protons observed by SWSE are plotted in Fig- 
ure 5. The magnetic field is plotted in Figure 6 
as the field magnitude, solar equatorial longi- 
tude <p, and solar equatorial latitude 0. 

For the sake of clarity, we have divided the 
lunar and near-lunar observations into three 
periods: from 0300 to 0543, from 0543 to 1030, 
and from 1030 to 1300. 
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Fie:. 2. tf components for ground-based mag- 
netometers for the April 9, 1971, geomagnetic 
storm. The Z component for Cambridge Bay 
is also given. 


From, 0300 to 051(3, Prior to the sc at 0428 
the high-latitude field strength was ~T0 y . 
The moon entered the plasma sheet, at 0400 as 
evidenced by the field strength decrease in the 
Explorer 35 magnetometer data (Figure 6). 
Behannon [1968] observed a delay of 12 dr 2 
min between the sc and a high -latitude tail 


field increase near lunar distance for the geo- 
magnetic storm of September 14, J 066, Thus 
we might expect to see an effect at the. lunar 
distance at ^0440. However, between 0441 and 
0447 the 180° shift in the magnetic field direc- 
tion indicates that the moon crossed the neutral 
sheet. A high-latilude tail field increase would 
not be observed at the instant of onset because 
of the masking effect of the neutral sheet. After 
0447, however, the magnetic field strength was 
8-11 y, but CP LEE data show that the moon 
was still in the plasma sheet. From pressure 
balance considerations we calculate a high- 
latitude tail field strength of 14 y, or an increase 
of 4 y Over the prestorm value, which was due 
to the sc propagating to the lunar distance. 

From 051(3 to 1030, After several encounters 
with the magnetopause the moon passed into 
the magnetosheath at 0543 and remained there 
until after 1000. This period coincides with that 



Fig. 3. Solar wind plasma and magnetic field 
data for April 9, 1971, by Explorer 33. 
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Fig. 4. Differential fluxes measured by the 
50- and 600-ev electron and 1-kev proton channels 
of CPLEE during the hours 0300 to 1300 of April 
9, 1971. 

when the solar wind was flowing from cast of 
the sun. Within the period, two regions of ob- 
servations are evident: magnetosheath and 
boundary layer. Characteristics of the magneto- 
sheath are discussed below. During times of 
magnetosheath observations we note that the 
solar wind flow direction was >8° from east of 
the sun. Encounters with the boundary layer, 
characterized by rapid decreases in density and 
flow speed, were observed at 0607, 0643, 0714, 
and 0800. Allowing 15 to 20 min for solar wind 
information to propagate from the position of 
Explorer 33 to the moon, we find that, these 
encounters coincide with periods when the solar 
wind flow was ~5° from east of the sun. Since 
none of the other solar wind characteristics 
show a consistent variation, we conclude that 
at lunar distance the position of the magneto- 
pause is strongly dependent on the flow direc- 
tion of the solar wind. That is, the dynamic 
pressure of the solar wind so completely domi- 
nates magnetic field tension that at lunar dis- 
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tance the tail responds to the solar wind by 
aligning itself with the flow. 

During the hours that the moon was in the 
magnetosheath the SWSE measured particle 
densities between 30 and 70 cm" 3 until the 
rapid decrease at ~0945. The magnetosheath 
flow speed remained near 300 km/sec. The 
thermal speed of the protons stayed between 
50 and 70 km/sec until 0945, when it increased 
to over 100 km/sec. 

The fluxes observed in the low-energy elec- 
tron channels of CPLEE have been fit to 
Maxwellian spectra and show temperatures 
between 15 and 20 ev. Temperatures remained 
within this interval for the entire period and 
are comparable to the 10- to 12-ev temperatures 
that we observe in the quiet time sheath. This 
behavior is consistent with observations by 
Montgomery et al. [1968] that, the solar wind 
electron temperature remains fairly constant 
despite wide variations in proton temperatures. 
We would stress that the temperature given 
above is an upper limit. The effect of a positive 
surface potential due to photoeleetrons would 
lead us to overestimate the temperature by a 
few electron volts. On the other hand, the 
medium energy electrons 0.5 < E < 2 kev 
show spectra typical of what we normally 
observe in the plasma sheet rather than the 
magnetosheath. We suggest that the origin of 



Fig. 5. Proton density, flow speed, and thermal 
speeds measured by SWSE during the hours 0300 
to 1300 of April 9, 1971. 
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Fig, 6. Magnetic field strength as well as the 
solar equatorial azimuthal and latitudinal direc- 
tions as measured by Explorer 35 for the hours 
0300 to 1300 of April 9, 1971. 

this population is the magnetosphere. The gra- 
dient drift of plasma sheet electrons tends to 
carry them to and possibly through the dusk 
magnetopause. 

The SWSE and CPLEE entered the magneto- 
sheath during a magnetometer data gap. After 
0545, when data resumed, the field magnitude 
was slightly larger than it was at 0537, and its 
orientation was significantly different: south- 
ward {0 between —40° and —75°) with a 
distinct dusk to dawn component (<p ~ 250°). 
This configuration characterized the magneto- 
sheath field data up to the data gap beginning 
at 0714. 

Between 0547 and 0714 a remarkable coin- 
cidence is found between the interplanetary 
and magnetosheath fields. In both, the field 
strength was ~ 15-18 y, with longitudinal direc- 
tions near 270° and strong southward compo- 
nents. The stronger southward component in 
the magnetosheat.h field was probably due to 
the effects of field lines draping over the mag- 
net otail. The similarity of interplanetary and 
dusk side magnetosheath fields at lunar distance 
is consistent with observations by B shannon 
and Fairfield [I960]. 

From 1030 to 1300. During the hours after 
the moon reentered the magnetotail a multiplic- 
ity of phenomena was observed. The high- 
latitude magnetic field strength was in excess 
of 30 y. There were at least eight crossings of 
the neutral sheet. In the plasma sheet, CPLEE 


observed particle densities of the order 1 cm _a , 
a factor of 10 higher than our observed quiet 
time density. The electron and proton tempera- 
tures were ^300 ev and kev. 

Because of the information afforded concern- 
ing the boundary layer between the tail and 
sheath, special attention is directed to phenom- 
ena observed between 1205 and 1300. First, we 
compare the magnetic fields observed by Ex- 
plorers 33 and 35. Their magnitudes are com- 
parable. However, although the interplanetary 
field varies in longitudinal and latitudinal angle 
ranges of 135° to 250° and 0° to —40°, the 
field observed at Explorer 35 was in the ranges 
0 C to 45° and 0° to —15°. Since the field ob- 
served by Explorer 35 shows features of neither 
the interplanetary medium nor the magnetotail, 
we postulate that Explorer 35 spent this hour in 
a transition region between the tail and sheath. 
A model that reproduces the field characteris- 
tics observed by Explorer 35 during this hour 
is developed in the appendix. 

From 1132 to 1137 and from 1205 to 1228 
the 40-ev electron channel of CPLEE was at 
the photoelectron counting level, and the 500-ev 
channel was at background. Generally, this 
situation would indicate that the moon was in 
the high -latitude tail. However, the flux ob- 
served in t-he 200-ev channel of cm" 5 sec" 1 
ster -1 ev" 1 is a factor of 6-10 higher than the 
photoelectron level [Reasoner and Burke , 
1972]. In both instances the magnetometer 
showed field orientations of the boundary re- 
gion. This result suggests that there is a mech- 
anism in the boundary region responsible for 
the separation of ~200-ev electrons from the 
plasma sheet population. We return to this 
point, below in the discussion. 

Magnetosheath ions were observed by SWSE 
and CPLEE for the final time between 1230 
and 1240, a time when the solar wind flow 
was ^5° from east of the sun. However, the 
increased ion flux was not accompanied by a 
proportional increase in magnetosheath elec- 
trons. Observed electron fluxes are more typical 
of the plasma sheet population. If, as was 
argued above, both Explorer 35 and the moon 
spent this period in the boundary layer, the 
boundary at this time was at least 3000 km, 
the Y* nt separation of Explorer 35 and the 
instruments. 
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Discussion 

The position of the magnetosphere is deter- 
mined by the combined effects of the dynamic, 
thermal, and electromagnetic pressures of the 
solar wind. Our observation of the magneto- 
pause changing position in response to changes 
in the solar wind flow direction enables us to 
decouple dynamic from thermal and electro- 
magnetic effects. 

Bekannon and Ness [1966] have argued that 
the storm time increase in- the high-latitude 
field strength at 30 R E is due more to the 
carrying of flux from the front side into the 
tail than to thermal pressure from the magneto- 
sheath plasma. We find that at lunar distance 
the converse is true. 

By assuming constant pressure across the 
tail we can use the high -latitude field strength 
and ignore the magnetospheric plasma in the 
pressure balance equation 


Bl 

8tt 


Bl 

8t r 


+ Pi, + Pes 


The subscripts t and s represent tail and sheath 
quantities. Although the moon was in the mag- 
netosheath between 0600 and 0713, we observed 
plasma densities of ^65 cm' 3 , ion thermal 
speeds of ^65 km/sec, and an electron tem- 
perature of ~15 ov. These give ion thermal 
pressures of = 1.3 X 10" 

dyne/cm 2 and an electron thermal pressure 
— N?kT,, — 1.5 x 10" dyne, /cm 3 . During 
this time the magnetic field in the sheath was 
about 18 y, mostly in the ~Z direction, and 
thus a magnetic pressure of about 1.3 X 10" 
dyne/cm* is given. From this we calculate a 
total high-latitude field strength of about 32 y 
necessary t.o balance the combined magnetic 
field and particle pressure in the magnetosheath. 

Assume that we can represent the magneto- 
tail at lunar distance as a cylinder with the 
field pointing toward the earth in the northern 
half and away from the earth in the southern 
half. If such a field were compressed owing to 
an increase in external pressure, the conserva- 
tion of flux equation shows that B X R* = B*R?. 
Subscripts 1 and 2 refer to precompression and 
postcompression values. 

During the period 0600-0700 the solar wind 
had a velocity of ~450 km/sec with a direc- 
tional flow of about 5° from the east of the 


sun. Thus an observer in the earth's frame of 
reference would see the solar wind as flowing 
almost, radially from the sun. The radius of 
the cylinder is then approximately the square 
root of the sum of the squares of the Y tnt and 
Z Km positions of the boundary. 

On the basis of five lunar passages through 
the magneto tail we have determined the average 
dawn and dusk magnetopause positions. On 
the dawn side we find (Y sm 2 -f Z s j) l/Sl = 23.8 
± 3.7 U B} and (Y,J + Z*m*) t/S = 28.6 =fc 2.8 
R E on the dusk side. Thus au average tail radius 
of about 26.2 d= 2.3 R E is obtained. The error 
limits represent the standard deviation of the 
ensemble of data points and not the error 
associated with a single measurement. This 
value is within the limits found by Mihalov et 
ol . [1970]. Wc assume the value of ^26.2 R E 
as the quiet time radius of the tail and use 
Apollo 14 ephemeris data that show storm 
time positions of the magnetopause at 1200 to 
be [Y,l + Z a ,*y f * — 16 R e . Prior to the storm 
the high-latitude field strength was ^10 y. Flux 
conservation gives a value of 

B = 10(26. 2 2 /16) = 27 y 

Thus, of the 32-y field required for pressure 
balance, 27 y can be accounted for by thermal 
compression of the tail. Presumably the other 
5 y must come from flux that is carried from 
the front side to the tail as a result of the en- 
hanced solar wind. Using the Mead and Beard 
[1964] model of the magnetosphere, Sugiura et 
al. [1968] have calculated the magnetic field 
increase in the tail. The change of magnetic 
field strength AB is related to the prestorm 
values of the flux in the tail as well as pre- 
storrn and intrastorm solar wind density and 
velocity by 

= B { [{N,/N x ) h/ \V 2 /V,y n - 1 ] 

Again, B L — 10 y, A 7 , = 15 cm -8 , iV 3 ~ 60 cm' 3 , 
K — 350 km/sec, and V * ~ 450 km/sec. Thus 
the total change in the tail field due to flux 
transfer from the front is A B ~ 4.4 y. 

Since the transference of flux from the day- 
side to the night side takes place before the tail 
is thermally compressed [Sugiura et al ., 1968], 
it mi glit seem more proper to add 4.4 y to the 
10-y prestorm high-latitude tail field and then 
compress the tail. The 32-y field strength re- 
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quired by pressure balance implies that the tail 
radius should be R z — 20.2 (14/32) 1/2 =; 17 R E 
We see no way of distinguishing between a tail 
radius of 10 and 17 22*; a solar wind angle 
error of 1° would account for this difference. 
In either case it is obvious that most of the 
high-latitude tail field strength enhancement 
during the storm's main phase is due to thermal 
and magnetic compression of the tail. 

Two rather obvious objections can be raised 
against the analysis given above: it compares 
magnetosheath observations made at ~0600 
with rnagnetotail observations made at ^1200, 
and it ignores the effects of tangential magnetic 
stresses due to the strong southward component 
in the interplanetary field. 

The first difficulty disappears when one con- 
siders the fact that, although the solar wind 
particle density is loss at 1200 than it was at 
0000, the temperature has increased, and the 
magnetic field is about the same 15-18 y. Thus 
the sum of thermal and magnetic pressures is 
about the same at 0600 and 1200. 

It is true that flux has been carried from the 
front side to the tail because of tangential 
stresses. However, it is also true that this effect 
drives the magnet.ospheric convection system in 
a complicated way so that flux is carried back 
to the front side. On a long time scale, as 
much flux ius is carried to the tail by tangential 
stresses is also carried back again to the front. 
How much is carried to the. back on a short 
time scale is unknown. The fact that our calcu- 
lation using the Mead-Beard model [Mead and 
Beard f 1964] and then compressing the tail 
gives results consistent with our calculation of 
the tail radius and observed high-latitude field 
strength suggests that it would introduce a 
correction of a few gammas at most. 

Model for boundary thickening. We now 
consider the effects of compressing a cylindrical 
magnetotail of radius R L , high-latitude field 
strength and plasma sheet thickness B, to 
values R 2f B 2f and I/». If R 3 = R t /a, flux con- 
servation requires that B 2 — For simplic- 
ity, we assume that L 2 — LJa and that the 
field gradient across the plasma sheet is of order 
B/L in the northern and southern lobes. 

The gradient drift of plasma sheet electrons 
that do not cross the neutral sheet is of order 
V 0 = fx/gL } where p is the magnetic moment 


of the electron. The drift is directed toward 
the dusk side magnetopause. If the magnetic 
moment is conserved during compression, the 
ratio of the postcompression to precompression 
gradient drift for a particle of given energy 
and pi tell angle, V ff5 /V^, is of order a. If the 
compression is isothermal, V ff JV 0l is of order a. 

For the April 9 storm, a zz 1.65, and the 
ratio V,.yv gl is bounded: 

1.65 < VjV 9l < 4.5 

The flux of magnetospheric plasma into the 
boundary is of order nVg with n 2 ~ lOn, ; the 
flux into the boundary increases by a factor of 
between 16 and 45, depending on the nature 
of the compression. 

The effects of a magnetosheath flow of plasma 
on the boundary layer between the sheath and 
the tail have been xSludied by Parker [1967a, 6], 
Lercke [1967], Su and Sonnerup [1971], and 
Willis [1970]. Because of the different gyro- 
radii of electrons and protons a charge separa- 
tion region will form at the magnetopause. 
However, unless electrons from the ionosphere 
or magnetosphere also move into the region to 
provide quasi- charge neutrality, the boundary 
will tend to thin in response to the strong 
electric fields in the boundary layer. 

If there is a current along magnetic field 
lines in the boundary layer, a magnetic field 
perpendicular to the main field is generated. 
Parker [19676] has shown that, unless the 
solar wind electron pressure exceeds a certain 
critical level, the magnetic pressure will carry 
the newly generated field out into the magneto- 
sheath where it is oonvected away. A small-scale 
nonequilibrium situation exists that gives the 
boundary a thickness between 1 and 10 s km. 
On the other hand, using a resistance-capaci- 
tance circuit analogy, WiUis [1970] has argued 
that, unless the solar wind is constant for a 
full day, the unstable boundary layer cannot 
form by way of ionospheric electrons moving 
to the magnetopause. During moderately quiet 
times we expect that the tow density of the 
plasma sheet found at lunar distance would 
not be able to provide enough particles to estab- 
lish charge neutrality in the boundary. Flectric 
fields would then insure a narrow boundary 
with little or no charge separation. In this 
case there could be little parallel current in the 
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boundary. This thin boundary is consistent with 
observations of sharp transitions in the mag- 
netic field signature [Mihalov et al. y 1970] and 
in low-energy electron fluxes \Bwke and Rea- 
sorter , 1972] passing between the sheath and 
tail. Recently, thick boundary layers, marked 
by spatial gradients iri the proton flux, have 
been reported by Intriligator and Wolfe [1972] 
and Howe and Siscoe [1972]. Since in the 
Parker model the thickness of the boundary 
is defined by the separation between magneto- 
sheath ions and electrons, without concomitant 
magnetosheath electron data, it is not possible 
to say whether the observed 'thick’ boundaries 
are 'thick' in Parker’s sense of the word. 

During the magnetic storm the enhanced 
flux of plasma sheet electrons into the boundary 
provides charge neutrality and allows magneto- 
sheath protons to penetrate further than mag- 
netosheath electrons. A current along the field 
lines in the boundary should develop, and the 
boundary thickness increase [Parker, 19676]. 
A calculation in the appendix shows that the 
observed plasma characteristics would produce 
perturbations in the magnetic field similar to 
those observed by Explorer 35. 

Further, the requirement of charge neutrality 
in the boundary layer is not absolute. As long as 
the electric field is weak enough that it does not 
drag magnetosheath electrons into the boundary, 
the requirements for a thick boundary layer are 
satisfied. Thus, as plasma sheet electrons 
gradient-drift into the boundary layer, they 
undergo a E x B drift that carries them in the 
ztZ directions in the northern and southern 
lobes. The effect is to increase the thickness of 
the plasma sheet near the boundary. Because 
the gradient drift is energy dependent, the time 
for crossing the boundary layer is longer for 
lower-energy electrons. Figure 7 shows the 
trajectories of the guiding centers of electrons 
with energy E { and E 2} where E { > E 2} in the 
boundary layer. 

We submit that the above-mentioned mech- 
anism for the transport of electrons from the 
tail into the sheath is responsible for the electron 
spectra observed by CPLEE in the magneto- 
sheath. In Figure 8 (a, 6), typical electron 
spectra observed in the sheath and in the plasma 
sheet are shown. With the exception of the 
flux depletion in the 200-ev channel the spectrum 
| in Figure 8a resembles a superposition of 


plasma sheet and magnetosheath populations. 
Evidently, >500-ev electrons can get through 
the boundary layer where 200-ev electrons 
cannot. 

From 1132 to 1137 and 1205 to 1228 while 
the moon was in or near the boundary layer, 
CPLEE observed 200-ev electrons with fluxes 
of ' 10* cnT* sec -1 ster" 1 ev" 1 , whereas the 40-ev 
channel was at x>hotoelectron level and the 
500-ev channel at background. If this flux 
value is inserted into Figure 8a, as it is in 
Figure 8c, the total spectrum does indeed come 
to bo a superposition of plasma sheet and 
magnetosheath populations. 

Alternative explanations of the thick bound- 
ary suggest themselves. If, for example, there 
is a connection between the interplanetary and 
tail magnetic fields, magnetosheath plasma 
could flow into the boundary layer. Another 
possibility is that there exists a scattering 
mechanism that causes magnetosheath particles 
to diffuse a finite distance into the tail. 

The first counterproposed model implies the 
existence of stresses at the boundary as opposed 
to our model where interplanetary field lines 
slide freely along the magnetopause. Then too, if 
such a connection exists, both protons and elec- 
trons should have access to the boundary layer. 
This view is in contradiction to our observation 
of magnetosheath protons without corresponding 
increase in the flux of magnetosheath electrons. 

A difficulty exists for the second model in 
devising a scattering mechanism that would 
permit protons to diffuse across lines of flux 
without carrying magnetosheath electrons with 


2 



Fig. 7. Intersection of the dusk side boundary 
layer and the plasma sheet. The different tra- 
jectories of electrons are due to the slower 
gradient drift of the lower-energy electron. 
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Fig. S. (a, b) The CPLEE electron differential flux spectra (a) observed in the magnet.oshea.tli 

and (6) observed in the plasma sheet. 


them, (our observation). If such a scattering 
mechanism does exist, it can oniy be viewed 
as a triggering device for the Parker mecha- 
nism. 

To summarize, during the storm time com- 
pression of the tail the flux of plasma sheet 
electrons gradient-drifting into the boundary 
is enhanced. This enhancement allows a 
separation between magnetoshea t.h electrons 
and protons to develop. The deeper penetrating 
magnetosheath proton flux establishes a current 
parallel to the magnetic field lines, and the 
resulting unbalanced electromagnetic forces 
cause the boundary layer to expand. The resid- 
ual charge separation electric field that exists 
in the boundary iayer produces a drift of 
plasma sheet electrons away from the neutral 
sheet. Because of the energy-dependent differ- 
ence in time needed for plasma sheet electrons 
to cross the boundary, <200-ev electrons drift 
in such a way as not to get through the bound- 
ary. Electrons with energies >500 ev do got 
through. The spectrum of electrons observed 
in the magnetosheath is a superposition of 
typical magnetosheath electrons and the >500- 


ev plasma sheet electrons that penetrate the 
thickened boundary layer. 

Summary and Conclusions 

During the geomagnetic storm of April 9, 

1971, Explorer 33 was in the solar wind on the 
dawn side of the earth, and the moon was in 
a position normally welt inside the dusk side of 
the magnetot.aiL A rise in solar wind particle 
density and flow speed as well as magnetic 
field strength occurred at 0435, 7 min after 
the observation of a sc by afternoon sector 
geomagnetic stations. At 0543, as the inter- 
planetary field increased to 18 y and the flow 
changed to >5° from cast of the sun, the 
magnetosheath appeared at the position of 
the moon. The moon remained in the sheath 
until -—'1000 when the flow direction shifted to 
the west. From our analysis of the history of 
the storm we have shown : 

1. The sc and geomagnetic storm were caused 
by a fast-moving stream cor ota ting in the solar 
wind. Our identification of the leading edge 
as a tangential discontinuity propagating near 
the garden hose angle in the solar wind based € 
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on computation of the discontinuity is con- 
sistent with the time delay between the ob- 
servations of the sc in the afternoon sector and 
the arrival of the disturbance at Explorer 33. 

2. The observations of the magnetopause 
in response to changes in solar wind flow 
direction, although no other solar wind param- 
eters varied in a consistent way, suggest that 
the position of the magnetopause is strongly 
dependent on the attack angle of the solar wind. 

3. The tail field strength shortly after the 
sc was indirectly measured to increase from 
^10 to — 14 y and is consistent with the pre- 
dictions of the Mead-Beard model. During the 
main phase the tail field strength was directly 
measured to be between 28 and 34 y. This 
finding is consistent with a thermal and mag- 
netic compression of the tail from — '26 to ^16 

n E . 

4. Plasma sheet electrons were lost th rough 
the magnetopause and became a component of 
the sheath flow. 

5. By simultaneous observations of the 
magnetopause with spatially separated detec- 
tors, its thickness was observed to be >3000 km 
on at least one occasion during the storm. 

6. Plasma and magnetic field observations 
are consistent with a model first proposed by 
Parker for currents flowing paralled to the 
boundary. 

Appendix 

In this appendix we present a calculation to 
check the consistency of plasma and magnetic 
field observations with the boundary layer 
model. 

At -—1200 the magnetic field indicates that 
Explorer 35 left the northern high -latitude lobe 
of the magnetotail. The main component of 
the field was toward the earth and had a 
strength of —TO- 15 y. The Y component of 
the field was 5-10 y and directed in the ,-j-K 
direction; The Z component was positive with 
values of 1-3 y. Between 1230 and 1240 the 
particle flow velocity was — '250 km/sec, and 
the density of plasma sheet electrons was 
—>0.1 cm 4 . At — 1245, Explorer 35 encountered 
the neutral sheet and then returned to the 
northern lobe of the tail. 

For the sake of simplicity, wo approximate 
the intersection of the magnetosheath and 
plasma sheet (Figure 7) by a cylinder of radius 



ENERGY (eV) 

Fig. Sc. Composite spectrum made by in- 
serting the 200-cv differential flux observed in 
the boundary layer. 

U (Figure 9a). The magnetic field at a point r 
within the cylinder is 

B(r) = — f j(r')r' dr' 
rc J u 

If we assume a current distribution 
j(r') = jo(l - r'/R) 

then 

bw _%m 

The term (r/R) [1 — (2r/3ft)] is plotted as 
a function of r/R in Figure 96 and has a 
maximum of 0.37 at r/R = Here jo .= 
UaC/V, where n a is the density of unbalanced 
magnetosheath protons; this number has to be 
of the order of the density of electrons supplied 
by the plasma sheet. Furthermore, since B v 
and B* are both positive oiily in the upper 
right-hand quadrant of Figure 9a, w r e can put 
a lower limit of 3000 km on R . Putting a 
value of {r/R)[l — (2r/3R)] = y 3f we get 
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X 

R 

Fig. 9. («) Schematic representation of inter- 
section of the magnetosheath and plasma sheet 
as a current-carrying cylinder, (b) Plot of the 
function (r/R) [1 — {2r/3I2)j as a function of 
r/R . 

B > 3.5 y. Since we have probably underesti- 
mated VR (by as much as a factor of 2), the 
calculated perturbation magnetic field is con- 
sistent with its observed strength. 

This calculation does not show that the 
boundary layer mode! is cor reel, only that it 
is not inconsistent with the data. 
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The plasma sheet at lunar distance is investigated with the use of data from the charged 
particle lunar environment experiment (CPLEE), complemented with data from the Explorer 
35/ARC magnetometer. It is shown that the presence of the lunar surface docs not appreciably 
affect measurements of the plasma sheet characteristics by the lunar-based CPLEE instru- 
ment. In particular, the lunar surface generally does not shadow plasma sheet particles. This 
may be due to rapid random, passage (>40 km/sec) of magnetotail field lines with respect 
to the lunar surface or to diffusion of plasma sheet electrons into the flux tubes in contact 
with the lunar surface. The plasma sheet is generally observed as a rapid increase in ob- 
served particle fluxes and a simultaneous decrease in field strength. Analysis of the CPLEE 
data shows that, for the plasma sheet at lunar distance, typical quiet time parameters are 
n = 0.10 ± 0.05 enr 3 , JcT e — 200 ± 50 ev, and kTi = 2.5 ± 0.75 kev. The typical ranges for 
the parameters are 0.05-0.20 chi ' 8 for n, 175-325 ev for kT e , and 1-5 kev for kTi. The typical 
total (plasma plus magnetic) energy density is in the range 150-350 ev/cm 3 . A statistical 
analysis of the CPLEE data shows that the plasma sheet in the midnight sector has a 
thickness of 5 R n ± 2 Re. Geomagnetic activity reduces the probability of encounters between 
the moon and the plasma sheet. This finding is consistent with the concept that geomagnetic 
activity increases the variance of the plasma sheet from its average location and/or decreases 
its thickness. 

The plasma sheet and the magneto tail have 
been explored extensively and reported upon 
out to a geocentric distance of ~30 R B [Bame 
et d. 7 1967; Vasyliunas, 1968; Gringauz, 

1969; Fairfield and Ness, 1970 ; Hones et al. y 
1971a] but have not been as well explored 
for large geocentric distances. The most com- 
plete reports of the plasma sheet and the 
magnetotail at lunar distance (~6Q R B ) 
are the report by Meng [1971] of high-energy 
(E > 20 kev) electron data, the statistical sur- 
veys by Meng and Mihalov [1972a, 6] of flux- 
gate magnetometer data, and the report by 
Nishida and Lyon [1972] of low-cnergy (0. 1-3.0 
kev) electrons measured with a Faraday cup. 

In a preliminary report, Burke and Reasoner 
[1972] stressed that papers by Nishida et al. 

[1969] and Prakash and Binsach [1971] to the 
effect that the plasma sheet at lunar distance 
is similar in geometry and spectral characteris- 
tics to what had been observed at 18 R E were 
in disagreement with charged particle lunar en- 
vironment experiment (CPLEE) observations. 

The absence of plasma observations for long 
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periods of time seemed to indicate that, if the 
plasma sheet extended to the moon, it had a 
thickness of less than 1 R E . This more complete 
analysis of CPLEE data indicates substantial 
agreement with Nishida and Lyon [1972] on 
the question of geometry. Disagreement contin- 
ues on spectral characteristics. 

The plasma sheet at lunar distance is studied 
with the use of data obtained with the CPLEE. 
The CPLEE was deployed on the lunar surface 
on February 5, 1971, during the Apollo 14 
mission and has been described by 0 7 Brien and 
Reasoner [1971] and Burke and Reasoner 
[1972]. It contains two electrostatic analyzers, 
each having an aperture of 4° by 20°. The 
aperture of analyzer A is directed toward lo- 
cal lunar vertical, and the aperture of ana- 
lyzer B is directed 60° from lunar vertical 
toward lunar west. The look directions in solar 
ecliptic coordinates during the lunar cycle are 
shown in Figure 1. The instrument can obtain 
a 15-point electron and ion spectrum in the 
Tange 40 ev to 20 kev and a background mea- 
surement in 19.2 sec or a five-point electron or 
ion spectrum in 2.4 sec. 
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Fig. 1. Lunar orbit projected into the ecliptic plane, showing the position of the CPLEE 
with respect to the geomagnetic environment and the look directions of the physical analyzers. 



Y 

Fig. 2. Projection of the lunar orbit onto the 
Y-7, plane in SM and SE coordinates, showing 
the passage of CPLEE through the magneto tail 
for February and March 1971. 


General Survey op the Data 

For illustrative purposes the CPLEE and 
the Explorer 35/ARC magnetometer data for 
the passage of the moon through the magneto- 
tail during February 1971 are presented here in 
Figures 4-7. The solar ecliptic (SE) and solar 
magxietospheric (SM) Y-Z coordinates for the 
February and March 1971 magneto tail passages 
are displayed in Figure 2. The Kp indices for 
the February 1971 tail passage are shown in 
Figure 3. The CPLEE data .are presented as 
raw count rates (counts/ 1.2 sec) versus time. 
Background count rates have not been sub- 
tracted from the data presented. Generally, 
the background is less than 10 counts/ 1.2 sec 
and is often less than 3 counts/1.2 sec. The 
Explorer 35/ARC magnetometer data are pre- 
sented in solar equatorial polar coordinates. 
The many gaps in the display are all due to 
gaps in the available data. The large displace- 
ment of <p from the expected values of — 0° and 
—'ISO 0 on February 10 and 11 (Figures 6b and 
7b) is due to a temporary bias added to the 
data by the on-board computer because of a 
total lunar eclipse that occurred on February 
10 between —0600 and — '0900 hours. At other 
times the field direction angles have an accuracy 
of ±:10° for B > 8 y. The data shown in Fig- 
ures 4-7 are typical of data obtained during 
subsequent tail passages. 

There are three basic types of observation^ 
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K p Index 


Feb. 8 




HOURS 

Fig. 3. Histogram of Kp indices for February 8-11, 1971, when CPLEE was passing through 
the magnetotail. Solid lines above the histograms denote periods of magnetic disturbance 
in the auroral zone- 

present in the data for any single magnetotail significant ion count rate just prior to entry into 

passage. These types can be identified with three the tail (e.g., Figure 4). The field strength ob- 

regions of space: the dusk magnetosheath, the served in the magnetosheath by Explorer 35 is 2-6 

Jiigh-latitude magneto tail , and the plasma sheet, y and highly variable in direction with a tendency 

The magnetosheath on the duskside is ob- toward the garden hoso angle (<p ^ 135° or 

served by CPLEE principally as a region of 315°). Large values of 0 are common in the 

intense lowrenergy electron fluxes, e.g., the first sheath region, especially near the tail boundary. 

2 hours of February 8 (Figure 4). These mag- Tho high-latitude magnctotail is observed by 
. netosheath electrons observed by CPLEE are CPLEE as a region with no detectable charged 

consistent with n, = 1-10 cm -5 and kT e s= 10- particle flux except for a very steady flux of 

20 ev [Burke et ql., 1973], Because of the large electrons with energies between 40 and 200 cv 

angle between the ion flow direction and the look thaT have been identified as the high-energy 

direction of the CPLEE analyzers (^>40°), in- tail of the surface photoelectron distribution 

tense fluxes of ions are not seen in the dusk mag- [ Reasoner and Burke, 1972], The magnetic 

netosheath, although analyzer A often observes a field strength has a strength generally in the 





February 8, 1971 



February 8, 13 7 f 



Fig. 4, Data for February 8, 1971: CPLEE data in counts/1.2 see averaged over six or 15 cycles (~2 or min); all data shown are from analyzer A; 
Explorer 3 5/ ARC magnetometer data averaged over 13 £ cycles (81.8 sec) displayed in solar equatorial polar coordinates. 
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Fig. 6. C 



February II, 1971 



Fig. 7. Continuation in time of Figure 6 for February 11, 1971, All CPLEE data are for analyzer B except for the 2-kev ion channel, which shows the 
effect of solar UV photons entering analyzer A near local lunar noon, which was at 2240 hours UT, February 11, 1971. 
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range of 9-13 y; during geomagnetically dis- 
turbed times the field strength can be signifi- 
cantly greater. The high-latitude tail field is 
oriented approximately parallel or antiparallel 
to the earth-sun line (9 ~ 0° or 180°; $ ~ 
0°). An example of data taken in the high- 
latitude tail is shown in Figure 4 near 1800 
hours of February 8, 1971. 

The moon is considered to be in the plasma 
sheet region when CPLEE observes a flux of 
electrons significantly above the background or 
the photociectron count rate in the energy 
range of 100 ev to 2 kev and a flux of ions in 
the range 1-10 kev. Within the capabilities of 
the CPLEE instrument it has been determined 
that there are no anisotropies in the flux of the 
plasma sheet electrons. The ions were also 
measured to be isotropic most of the time but 
were distinctly anisotropic on a few occasions 
[c.g. ? Burke and Reasoner, 1973]. 

The plasma sheet is detected in the Explorer 
35/ARC magnetometer data by a decrease in 
the magnetic field strength and an increase in 
the variability in the field direction. This is due 
to the diamagnetic effect of the plasma sheet 
plasma and has been used by Meng and Mihdov 
[1972a, h] to locate indirectly the plasma sheet 
at lunar distance. The simultaneous data from 
Explorer 35 and CPLEE show that there is a 
definite correlation between the observations by 
CPLEE and those by Explorer 35 of the plasma 
sheet. There are some noticeable exceptions that 
arc discussed below. 

Lunar Surface Effects 

Because particles striking the lunar surface 
are absorbed, a particle shadow can develop iu 
tubes of flux connecting the earth and the moon. 
Lin [1968], Van AUen and Ness [1969], and 
Anderson and Lin [1969] have reported that the 
concept of lunar shadowing is valid for high- 
energy ( E > 20 kev) electrons. Anderson [1970] 
used the high-energy electron data to develop 
a s ingl e-p articlc-sh ad 0 wing model that could be 
used to estimate the speed of flux tubes crossing 
the lunar disk. Assuming the motion to be 
E x B drift, where E is a steady state convection 
electric field in the tail and where B ^ 10 7, 
Anderson [1970] determined that E < 0.50 
mv/m, or to be more precise that the speed of 
passage of flux tubes with respect to the moon is 
generally less than 50 km/sec. Based on a cross- 


tail potential drop of < 75 kev, the steady state 
convection electric field is estimated to be < 0.25 
mv/ m in the magnetotail. 

Let us assume that Anderson's model of par- 
ticle shadowing by the moon can be extended 
to electrons or ions of plasma sheet energy, and 
let us examine the consequences of this assump- 
tion. Depending on the direction of the drift 
motion, the distance between the CPLEE and 
tho edge of the lunar disk is 1200-2400 km. In 
order for 200-ev electrons, 2-kev electrons, or 

2- kev ions from the plasma sheet to reach the 
detector, the flux tubes must travel across the 
lunar surface at a minimum speed of 24-48, 
75-150, or 2-4 km/sec, respectively. McGuire 
[1972] has shown that significant particle 
fluxes would be observed in the shadowed region 
up to 1 gyroradius from the shadow edge cal- 
culated in the simple model. Thus in the pres- 
ence of an isotropically distributed plasma sheet 
population there is no expected shadowing of 
plasma sheet ions. However, if it is assumed 
that B ~ 10 y and E <£ 0.25 mv/m, then 
plasma sheet electrons with energy greater than 
500 ev should never reach CPLEE, owing to 
shadowing. Owing to diamagnetism, B is often 

3- 6 y, and electrons with energies up to 2 kev 
should occasionally reach CPLEE. Under these 
conditions, cutoffs in the electron energy spec- 
trum would be featured in the CPLEE data. 

Generally, the CPLEE data do not show any 
of the preceding characteristics. Most observa- 
tions of ion fluxes are accompanied by observa- 
tions of an electron distribution with energies 
from UK) ev to 2 kev. There is no clear evidence 
for a cutoff energy. Occasionally, 2-kev elec- 
trons are not observed at the same time as 
200-ev electrons, but these observations arc 
consistent with variations in electron tempera- 
ture and density and are not considered to be 
evidence of shadowing. Most observations of 
decreases in the magnetic field strength at 
Explorer 35, which are indicative of the pres- 
ence of the plasma sheet, occur at times that 
CPLEE observes fluxes of plasma sheet parti- 
cles. 

There are a few occasions when tho CPLEE 
data arc consistent with the lunar-shadowing 
model. On February 10 between 1600 and 1900 
hours (Figure 6) the magnetometer data show 
the presence of the plasma sheet, and CPLEE 
detected an anisotropic ion population but did - 
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not detect significant fluxes of plasma sheet 
electrons. Between March 10 at 2300 hours and 
March 11 at 0300 hours (Figure 8) a plasma, 
sheet with an isotropically distributed ion popu- 
lation was observed, but except for three brief 
intervals the plasma sheet electrons were not 
observed. These two sots of data are more con- 
sistent. with the particle-shadowing model than 
any other data obtained by CPLEE. However, 
the cutoff in the energy spectrum predicted by 
the model did not appear in the electron spectra 
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during the three brief recoveries of electron 
fluxes (Figure 8). 

Since lunar shadowing of particles does not 
significantly affect the flux of plasma sheet 
electrons reaching the CPLEE, it must be con- 
cluded either that flux tubes cross the lunar 
surface at speeds generally greater than 40 
km/sec or that the model of particle shadowing 
may not be extended to plasma sheet, electrons. 

The CPLEE data imply that the convective 
electric field is greater than 0.25 mv/m the 


Morch 10- 1 1 * 1971 



HOURS 


lug. 8. Data for March 10-11, 1971. Between 2300 and 0300 hours the Explorer 35/ARC 
magnetometer data and the CPLEE ion data, show the presence of the plasma sheet, yet 
the CPLEE electron data show almost no plasma sheet electrons. 
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vast majority of the time and is occasional! y 
greater than 1.0 mv/in. Rapid but randomly 
oriented motions in the tail that, when they are 
averaged over several minutes are consistent 
with E < 0,25 mv/m, would explain the 
CPLEE data in a manner consistent with the 
shadowing concept. In fact, random changes 
in the solar wind flow direction of the order of 
0.5°-l-,0°/min are commonly found in the 
available solar wind data (J. H. Binsack, pri- 
vate communication, 1972). Howe and Binsack 
[19721 have shown that on a statistical basis 
the magnetotail at lunar distance tends to align 
itself with the direction of solar wind flow. In 
a study of the geomagnetic storm of April 9, 
1971, Burke et al. [1973] showed that the posi- 
tion of the magnetopause responds quickly to 
changes in flow direction. The concept that a 
1° shift in the solar wind flow causes a 1 -R B 
shift in the location of the magnetotail at lunar 
distance implies that flux tubes cross the lunar 
disk with speeds of 50-100 km/sec, which is 
sufficient to remove particle shadows. 

Other particle drifts, such as those clue to 
curvature and gradients in the magnetic field, are 
of little significance relative to the ExB drift. 
In the plasma sheet at lunar distance (away 
from the neutral sheet), magnetic field lines show 
little curvature [ Behannon , 1968]. The gradient 
drift V g — (cw 1 /eB*)Bx'VB, where w ± = % m . 
v 2 cos 2 a and a is the particle’s pitch angle. At 
lunar distance, B ~ 10 y and is directed para llel 
or antiparallel to the eart-h-sun lino. The gradient 
in the magnetic field is of the order of B/L, 
where L is the half thickness of the plasma sheet. 
If L 3 then a particle with w ± = 1 kev 
experiences a gradient, drift of ^5 km/sec. 

If the flux tubes do not cross the lunar disk 
rapidly in a random fashion, then the shadowing 
model must be invalid for describing the motion 
of plasma sheet particles. The predicted void 
region may be filled in by diffusion of plasma 
sheet, electrons from other regions of the plasma 
sheet on time scales of <JJ50 sec. Diffusion is the 
result, of the electron temperature gradient, 
from hot electrons in the unperturbed plasma 
sheet to cold electrons across the field lines. 
The calculations pertaining to the rate of diffu- 
sion expected have not been done, but the maxi- 
mum rate is the Bohra diffusion rate, there 
being lateral movement, of 1 gyroradius each 
gyroperiod, which is many orders of magnitude 


greater than that needed to fill in the predicted 
void region [ Krall and Trivelpicce , 1.973] . Al- 
ternatively, the shadowing model may be invalid 
because the plasma sheet flux tubes at lunar 
distance may be filled from a continuous plasma 
source [Frank, 1971 ; Hill, 1973], 

Spectral Distribution of the Plasma Sheet 
Population at Lunar Distance 

The ion and electron spectra shown in Figure 
9 are representative of many CPLEE spectra 
in the magnetotail. These particular spectra arc 
presented because data were obtained during 
a long period of almost constant particle fluxes, 
simultaneous Explorer 35/ARC magnetometer 
data are available, a neutral sheet crossing did 
not occur while the data were being obtained, 
and geomagnetic activity was very low. 

Tho photoelectron electron spectrum is a very 
noticeable feature of the electron spectrum in 



Fig. 9. Typical CPLEE ion and electron spec- 
tra, in the plasma sheet. Spectra taken at 1946 
hours UT, February 9, 1971. 
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Figure 9. Basically, we treat the photoelectron 
spectrum as a background spectrum that must- 
be subtracted from the data in order to obtain 
the plasma sheet electron spectrum. 

The plasma sheet electron spectrum in Figure 
9 is approximately representative of a Maxwel- 
lian distribution with a number density of 0.10 
± 0.05 cm'* and a temperature of 200 dr 50 ev. 
Beyond 750 ev the electron spectrum is best 
fitted by a power law dstribution with an expo- 
nent of —3. This trend has been shown to con- 
tinue out to 15 kev [Chase et al., 1971 ; Ander- 
son et cd., 1972], 

A reasonable fit for the ion spectrum in Fig- 
ure 9, when it is assumed that all the ions are 
protons, is a Maxwellian distribution with n t = 
n e = 0,10 ± 0.05 cm * and &T t — 2.5 ± 0.75 
kev. Owing to the low flux of ions relative to 
the CPLEE sensitivity limits, the uncertainty 
in the analysis of the ion data is large. 

Both analyzer A and analyzer B generally 
observe identical ion and electron fluxes, indi- 
cating isotropic fluxes in the plasma sheet. This 
allows a determination of the plasma pressure 
in the vicinity of the moon. For the data shown 
in Figure 9 the electron and ion pressures are 
18 d: 3 and 250 =t 50 ev/cm 3 , respectively. The 
observed magnetic field strength in the plasma 
sheet was ^5 y. When a high-latitude field 
strength of 10-12 y is assumed, a plasma’ pres- 
sure of 190-300 ev/cm n is required to maintain 
pressure balance across the high-latitude tail 
and the, plasma sheet. Thus these data are con- 
sistent with the pressure in the plasma sheet 
determined by Meng and Mihalov [19726]. 
Meng and Mihalov [19726] miscalculated the 
magnetic field pressure by a factor of 3/2, but, 
when their technique and the data used in 
their paper are taken into consideration, it can 
be determined that the particle pressure in the 
plasma sheet is typically in the range 150-350 
ev/cm 3 . 

The range of the plasma sheet electron tem- 
peratures observed by the CPLEE is generally 
175-325 ev, there being one case of a substorm 
recovery temperature of ~475 ev [Rich et al. , 
1973]. Temperatures significantly less than 175 
ev have never been observed in the plasma 
sheet. The range of ion temperatures is 1-5 
kev. Higher temperatures are not observed at. 
lunar distance. Lower ion temperatures are not 
^ observed, owing to the sensitivity limits, but are 


consistent with some data showing a plasma 
with an electron component but no detectable, 
ion component. The range of the plasma sheet 
number density is 0.05-0.20 cm" 3 . Densities as 
large as 0.5 cm" 3 have been observed but only 
during a geomagnetic storm [Burke et al., 
1973]. 

Quite obviously, our observations are in 
disagreement with results published by Nishida 
and Lyon [1972]. Their density of 1 cm" 3 and 
hT e = 0.6 kev give a total electron pressure of 
600 ev/cm 3 . If we assume a proton temperature 
of as low as 1.5 kev, the total plasma pressure 
would be sufficient to balance a high-latitude 
field of ^25 y, which is two and one half times 
tho field strength shown in their own data. 

Average Spatial Configuration of the 
Plasma Sheet 

The location and the shape of the plasma 
sheet could be determined completely with a 
few passes through the tail if the plasma sheet 
location were stationary with respect to an 
appropriate set of coordinates and if its shape 
were time independent. However, neither con- 
dition is true of the plasma sheet. The magneto- 
tail, which carries. the plasma sheet with it, has 
been shown to respond to changes in the solar 
wind flow direction [Howe, 1971; Howe and 
Binsack, 1972; Burke et aL , 1973]. At lunar 
distance, the tail shifts approximately 1 R E for 
each 1° shift in the solar wind flow direction. 
This shift is shown by comparing the neutral 
sheet encounters reported by Nishida and Lyon 
[1972] with the distribution of the N-S com- 
ponent of the solar wind flow reported by Lyon 
et at. [1968]. The plasma sheet thickness is 
known to change in response to substorms 
[Hones et al., 1973]. 

The CPLEE data used for the plasma sheet 
survey used in this paper were obtained during 
the tail passages of February, March, and May 
1971; August, September, and December 1972; 
and January 1973 plus 6 hours of data from 
July 1972. During August and September 1972, 
only electron data are available. The data pre- 
sented here consist of 262.5 hours of electron 
data and 184.5 hours of ion data. The presence 
of the plasma sheet was determined with a time 
resolution of 1 min by examining the ion and 
electron channel considered to be the most sen- 
sitive to plasma sheet fluxes. The frequency 
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Fig. 10. Probability of encountering the plasma sheet as a function of Z S6 and Z am in 
l.he interval 15 Ru < Y , m < —5 Re, based on the CPLEE data. For comparison, the 
probability in the same interval determined by Meng and Mihctlov L1972u] is also shown. 
The probability is computed once from the set of all available data and once for the set 
of data obtained during the period of Kp < Re- 


distribution of plasma sheet encounters in areal 
elements extending 1 R K in the Z direction and 
20 Ru in the V direction was computed. Smaller 
areal elements were not used because of an 
insufficiency of data. The areal elements used 
arc chosen to be approximately symmetric 
about the average central meridian of the tail, 
i.c., extending from Y am — +15 R E to Y am = 
— 5 R e . Thus this paper concentrates on the 
problem of the thickness of the plasma sheet in 
the center of the tail. An increased thickness 
near the flanks has been demonstrated clearly 
by Meng and Mihctlov [1072a], although the 
amount of increase may be uncertain. In the 
center of the tail, problems not directly related 


to the plasma sheet, such as magnetopause 
effects, should be minimal. The data have been 
organized according to the SE and SM coordi- 
nate systems. 

The results of the survey of the CPLEE data 
are presented in Figure 10. The results are 
presented from the total data set and from the 
subset of data for low Kp(< 1+ ). For compari- 
son, the results of Meng and Mihalov [1972a] 
for approximately the same areal elements, 
representing 1048 hours of their data, are also 
shown. There are two distinct differences in the 
results of the two data sets. The first difference 
is that the probabilities of encountering the 
plasma sheet at various Z sm coordinates are 
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greater in the Explorer 35/ARC magnetometer 
<iata than in the CPLEE data. (The high prob- 
ability at Z em = 6 R e and 7 R E in the CPLEE 
data set is not considered to be of great signifi- 
cance, owing to the small amount of time spent 
in those areal elements.) The second difference 
is the width of the maximum probability peak, 
which implies the plasma sheet thickness is 
smaller in the CPLEE data set than in the 
magnetometer data set. 

The systematic decrease in the plasma sheet 
encounter probability for CPLEE data com- 
pared to that for the magnetometer data could 
have several possible causes: First, there is a 
difference in the surveying technique. The tech- 
nique that Meng and Mihalov [1972a] used 
tends artificially to increase by as much as 10 
min the duration of plasma sheet encounters. 
The increase is negligible for long plasma sheet 
encounters but has some significance here, since 
the median duration for a plasma sheet encoun- 
ter is between 30 and 40 min for both the 
CPLEE and the magnetometer data sets [J. 
Mihalov, private communication, 1972], The 
sensitivity of the magnetometer using the detec- 
tion criteria of Meng and Mihalov may be great 
enough to detect regions of the plasma sheet 
where the energy density is less than that 
which can reliably be detected by the CPLEE 
(~100 cv/cm*). The difference in probabilities 
may also be due to perturbations in the high- 
latiturle tail field, perhaps induced by the solar 
wind, that cannot be distinguished with the 
magnetometer data from plasma sheet encoun- 
ters. Bulk flow of plasma in the tail, especially 
in the sunward direction, can also be recorded 
as a plasma sheet type field disturbance by the 
Explorer 35 magnetometer, whereas detectable 
particle fluxes fail to reach the CFLEE analy- 
zers. 

The results of Meng and Mihalov [1972a] 
shown in Figure 10 indicate a peak in the 
plasma sheet encounter probability ^7 li Pj wide 
and centered about the Z sni . = 0 R B plane. The 
CPLEE results, especially the results for Kp < 
1 .+, indicate two peaks in the probability func- 
tion having a total width of ^5 R Fj . It has been 
found that the lack of plasma sheet encounters 
between Z sm 0 and 2 R E cannot be attributed 
to multiple plasma sheets, a statistical anomaly, 
or seasonal effects. It is proposed that the 
apparently skewed and bimodal distribution 


derived from the CPLEE data is due to inter- 
actions of the magnetotail with the .solar wind 
that cause an aliasing problem due to a match- 
ing of the lunar synodic month and the solar 
rotation period. Solar wind flow patterns coro- 
tate with the mean solar surface rotation, and 
the solar wind flow pattern over a few days of 
a solar period may indeed be bimodal. Aliasing 
is given further support by the very clear 
correlation of the plasma sheet encounters with 
Z sv , which according to Meng and Mihalov 
[1972a] should not exist. Aliasing is expected 
to have a strong effect on the organization of 
the data in SE coordinates, owing to the close 
matching of the lunar nodical period with the 
solar rotation period. The problem of aliasing 
does not appear in the data analysis of Meng 
and Mihalov because solar wind flow patterns 
change on a time scale of 4r-8 months, and they 
used 30 contiguous months of data. By contrast, 
the CPLEE data were obtained over a 3 -month 
span in early 1971 and a 5-month span in late 
1972. 

Although aliasing gives an odd shape to the 
distribution function in Figure 10, the average 
plasma sheet configuration can still be deter- 
mined with reasonable accuracy. Based on the 
width of the peaks in the CPLEE distribution 
function, the average plasma sheet thickness in 
the center of the inagnetotail is ~5 R i{ . A 
lower bound can be set by assuming that the 
probability for areal elements outside the ele- 
ments with well-determined probabilities is 0.0 
and b^y assuming that an ideal organization of 
the data would yield a range of Z for which 
the probability is 1.0 and outside of which it is 
0.0. The lower bound is thus 3 R B . An upper 
bound can be set by assuming that the proba- 
bility for areal elements in the range 6 < 

\Z gm \ < 12 R P is a small nonzero value (c.g., 
0.2) and by assuming that an ideal organization 
of the data would yield a square distribution 
with a maximum value of 0.8 (20% of the time 
the plasma sheet is assumed to be very small, 
owing to geomagnetic activity). The upper 
bound is then 7 R E - The average plasma sheet 
thickness for the center of the inagnetotail at 
lunar distance can be staled as 5 ± 2 R K . 

The average plasma sheet thickness in the 
center of the tail determined here is consistent 
with the estimates of 7 R K given by Meng and 
Mihalov [1972a] and 6 li E implied by Nishida 
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and Lyon [1972]. The thickness value of 11 R B 
given by Meng and Mihalov, using Murayama 
[1966] neutral sheet coordinates, is more indica- 
tive of the failure of the Murayama system to 
organize the data than it is descriptive of the 
plasma sheet configuration at lunar distance. 
The instantaneous plasma sheet thickness may 
be considerably different from the average 
thickness, independent of geomagnetic activity. 
Of course, geomagnetic activity may also have 
an effect on the thickness. There is an anti- 
correlation of the encounter probabilities with 
the Kp index. This anti correlation implies that 
during geomagnetically disturbed times the 
plasma sheet on the average is thinner and/or 
has a larger domain of possible locations than 
during quiet times. 

Summary 

The CPLEE detects approximately the same 
fluxes of plasma sheet particles from its loca- 
tion on the lunar surface as it would if the moon 
were not present. This is concluded because 
the observed fluxes of particles correlate well 
with periods of magnetic field decrease that are 
indicative of .the plasma sheet, and the concept 
of pressure balance is well satisfied from the 
observations. As a con sequence of this fact, 
shadowing of particles by the lunar surface 
because flux tubes are frozen into the plasma 
and closed to the particles near the earth by 
converging field lines is not. a valid model for 
describing the particle dynamics in the plasma 
sheet unless the flux tubes cross the lunar sur- 
face with speeds generally greater than 40 km/ 
sec. If the flux tubes arc moving rapidly, the 
motion must be changing rapidly in a random 
manner in order to be consistent with the 
estimates of the steady state convection electric 
field (<£0.25 mv/m). Solar wind induced flap- 
ping motion in the magnetolail may cause such 
motion. The alternative to rapid random mo- 
tions for describing the data is to assume that 
plasma sheet particles easily diffuse across field 
lines to maintain a homogeneous population. 

Typical plasma sheet ion and electron spectra 
show that plasma sheet parameters obtained 
from the data are generally in the range 0.05- 
0.20 cm * for n, 175-325 cv for kT c , and 1-5 
kev for A 7V By comparing these data with 
plasma sheet data obtained neaT X = —20 R B 
and reported by Bame et al [1967] and Rones 


et d. [1971a, 5], the average quiet time tem- 
perature of the plasma at both locations is 
found to be approximately equivalent, and the 
X component of the plasma pressure gradient 
in the plasma sheet is found to be due mostly 
to a gradient in the number density. A recent 
report of Vela data near the center of the 
plasma sheet gives kTi zz 6 kev and kT c ~ 1 
kev (S.-I. Akasofu et al., unpublished manu- 
script, 1973), which indicate that near the 
center of the plasma sheet the X component of 
the pressure gradient may also be affected by a 
temperature gradient. The plasma sheet obser- 
vations at lunar distance may be consistent with 
the concept that the plasma at lunar distance 
is identical with plasma several earth radii from 
the center of the plasma sheet at Vela distance. 
The range of temperature variations observed 
during geomagnetically active times at lunar 
distance is smaller than that reported for the 
plasma sheet close to the earth. 

The probability of encountering the plasma 
sheet in the eont.er of the magnetotail at lunar 
distance as a function of Z 9m or Z fe has been 
developed from the CPLEE data to investigate 
the configuration of the plasma sheet. A .system- 
atic difference in the probabilities obtained 
here arid (Lose reported by Meng and Mihalov 
1 1972a] may be due partly to a difference in 
technique and instrumentation but is probably 
due partly to magnetic field disturbances in the 
high-latitude tail mistaken for the plasma sheet 
and to lunar shadowing of plasma flowing sun- 
ward from the deep tail. The odd shape of the 
distribution function obtained with CPLEE 
data has been found to be due to aliasing in- 
duced by solar wind patterns co rotating with 
the mean solar rotation period, which matches 
the lunar nodical period and the synodic month. 
This finding does not interfere with determining 
the plasma sheet thickness as 5 ± 2 R E in the 
center of the tail. The anti correlation of the 
encounter probabilities with the Kp index im- 
plies that geomagnetic disturbances on the 
average decrease the thickness and/or that 
they increase the domain of locations for the 
center of the plasma sheet. 
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MEASUREMENT OF THE LUNAR PHOTOELECTRON 
LAYER IN THE GEOMAGNETIC TAIL 


DAVID L. REASONER and WILLIAM J. BURKE 
Dept, of Space Science, Rice University, Houston, Tex. 77001, U.S. A. 


Abstract. The Charged Particle Lunar Environment Experiment (CPLEE), a part of the Apollo 14 
ALSEP, is an ion-electron spectrometer capable of measuring ions and electrons with energies between 
40 eV and 50 keV. The instrument, with apertures 26 cm above the surface, has detected a photo- 
electron gas layer above the sunlit lunar surface, with energies ranging up to 200 eV. Experimental 
data for periods while the Moon was in the Earth's magnetotail for electrons with energies 
40 eV £ < 200 eV follow a power law spectrum j (£) — /u(E/Eo) -/i with 3.5 ^ p ^ 4. In the absence 
of photoelectrons with £'>200eV, we assume that the surface potential is at least 200 V. The 
modulation of this potential in the presence of intense plasma sheet fluxes has been observed. 

Numerical solutions for the variation of electron density and potential as functions of height above 
the lunar surface were obtained. The solar photon spectrum / (/ jv ), obtained from various experi- 
mental sources, and the photoelectron yield function of the surface materials, Y (kv\ are two para- 
meters of the solution. Energy spectra at the height of the measurements for various values of Y ( hv ) 
were computed until a fit to experimental data was obtained. Using a functional form Y (hv)^ 
= \Y{){kv — W)]j{Wj2) for 6cV< hv<9t\ and Y(hv) = Yo for (Av)>9cV where IT, the lunar 
surface work function, was set at 6 eV, we calculated a value of To = 0.1 electrons photon -1 . The 
solution also showed that the photoelectron density falls by 5 orders of magnitude within 10 m of 
the surface, but the layer actually terminates several hundred meters above this height. 


1. Introduction 

The general problems of photoelectron emission by an isolated body in a vacuum and 
in a plasma have been the objects of several investigations. For example, Medved 
(1968) has treated electron sheath formations about bodies of typical satellite dimen- 
sions. Guernsey and Fu (1970) have considered the properties of an infinite, photo- 
emitting plate immersed in a dilute plasma. Grobman and Blank (1969) obtained 
expressions for the lunar surface potential due to photoelectron emission while the 
moon is in the solar wind. Walbridge (1970) developed a set of equations for obtaining 
the density of photoelectrons as well as the electrostatic potential as functions of 
height above the surface of the Moon while the Moon is in the solar wind. By assuming 
a simplified form of the solar photon emission spectrum he could provide analytic 
expressions for these quantities. 

In this paper we report on observations of stable photoelectron fluxes, with energies 
between 40 and 200 eV by the Apollo XIV Charged Particle Lunar Environment 
Experiment (CPLEE). These observations, made in the magnetotail under near 
vacuum conditions, are compared with numerically calculated photoemission spectra 
to determine the approximate potential difference between ground and CPLEE’s 
apertures (26 cm). Numerically calculated density and potential distributions, when 
compared with our measured values, help us estimate the photoelectron yield function 
of the dust layer covering the Moon. 


R. J. L. Grard ( cd .), Photon and Particle Interactions with Surfaces in Space, 369-387. Alt Rights Reserved 
Copyright (g) 1973 by D. Re ide! Publishing Company, Dovdrccht-Holiand 
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2* The Instrument 

Complete descriptions of the CPLEE instrument has been given by O’Brien and 
Reasoner (1971) and Burke and Reasoner (1972). The instrument contains two 
identical charged-particle analyzers, hereafter referred to as analyzers A and B . 
Analyzer A looks toward the local lunar vertical* and analyzer B looks 60° from 
vertical toward lunar west. 

The particle analyzers contain a set of electrostatic deflection plates to separate 
particles according to energy and charge type, and an array of 6 channel electron 
multipliers for particle detection. For a fixed voltage on the deflection plates, a five 
band measurement of the spectrum of particles of one charge sign and a single-band 
measurement of particles of the opposite charge sign are made. The deflection plate 
voltage is stepped through a sequence of 3 voltages at both polarities, plus background 
and calibration levels with zero voltage on the plates. A complete measurement of the 
spectrum of ions and electrons with energies between 40 eV and 50 keV is made every 
19.2 s. Of particular relevance to this study are the lowest electron energy passbands. 
With a deflection voltage of —35 V, the instrument measures electrons in five ranges 
centered at 40, 50, 65, 90 and 200 eV. With +35 V on the deflection plates, electrons 
in a single energy range between 50 and 150 eV are measured. At the next higher 
deflection voltages of + 350 V, the energy passbands given above are scaled upward by 
approximately a factor of 10. 

3. Observations 

In this section we present data from the February 1971 passage of the moon through 
the magnetotail. Because these are so typical, the display of data from subsequent 
months would be redundant. At approximately 0300 UT on February 8 CPLEE 
passed from the dusk side magnetosheath into the tail. The five minute averaged 
counting rates for analyzer A y channel 1, at —35 V measuring 40 eY electrons are 
plotted for this day in Figure 1 . Almost identical count rates are observed in analyzers 
A and B during this period of observation. As CPLEE moves across the magnetopause 
the counting rate drops from —200 cycle -1 to the magnetotail photoelectron back- 
ground of —35 cycle -1 (1 cycle— 1.2 s). Enhancements at —0530 h and at —0930 h 
correspond to plasma events associated with substorms on Earth (Burke and Reasoner, 
1972). There is a data gap from 1000 to 1200 h. With the exception of the short lived 
(< 1 h) enhancements the detector shows a stable counting rate when the Moon is in 
the magnetotail. 

Our contention is that these stable fluxes observed in the magnetotail during 
periods of low magnetic activity are photoelectrons generated by ultraviolet radiation 
from the Sun striking the surface of the Moon. In support of this thesis we have 
reproduced the counting rates observed in the same detector on February 10 when the 
Moon was near the center of the tail (Figure 2). First, we note that the stable count 
level is the same at the center as it was when CPLEE first entered the tail. Secondly, 
from about 0500 to 1000 h the Moon was in eclipse. During this time we observe the 
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Fig. 1. 5-ni in averaged counting rates for CPLEE, analyzer A, channel 1 at -35 V, measuring 
40 eV electrons on February 8, 1971. After 0300 UT counting rates fell from high magnetosheath 

to stable pbotoeiectron levels. 



Fig. 2. 5-min averaged deep tail counting rates of 40 eV electrons oil February 10, 1971. The lunar 
eclipse (0500-0900 UT) is marked by vanishing photoelectron counting rates. 
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counting rates go to zero. As the moon emerges from the Barth s shadow, the counting 
rates return to their pre-eclipse levels. If the stable low energy electrons were part of 
an ambient plasma, rather than photoelectrons, the counting rates would not be so 
radically altered as the moon moved across the Earth’s shadow. 

It could be argued that the observed counting rates were due to photons scattering 
within the detectors themselves and not due to external photoelectrons. This however 
is not the case. Preflight calibrations with a laboratory ultraviolet source showed 
enhanced counting rates only when the angle between the look direction of the detector 
and the source was less than 10°. Given the 60° separation between the look directions 
of analyzers A and B, it would be impossible for the Sun, essentially a point source, to 
produce identical counting rates in both analyzers simultaneously. 

A typical spectrum of photoelectrons shown in Figure 3 was observed by analyzer 



Fig. 3. Typical photoelectron spectrum observed by CPLEE at the lunar surface in the 

high latitude magnetotail. 
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A at —0400 h on February 10, shortly before the Moon entered pemimbral eclipse. 
The dark line marks the differential flux equivalent to a background count of one per 
cycle in each channel. (Channels 1-5.) For all five channels, with the deflection plates 
at — 35 V, the differential flux is well above this background level. During geomagne- 
tically quiet times no statistically significant counts are observed when the deflection 
plates are at — 350 or — 3500 V corresponding to electrons with E> 500 eV (Burke and 
Reasoner, 1972). 

With the exception of periods of ultraviolet contamination in analyzer A, we always 
observe nearly the same counting rate due to photoelectrons in analyzers A and B . 
For all purposes, we can say that the spectrum displayed in Figure 3 is just as typical 
as for analyzer B . We have found no case of anisotropy in the photoelectron fluxes. 
In all cases too, we found that the photoelectron spectra observed in both analyzers 
were close to a power law dependence on energy. If we write the differential flux in the 
form j(E)—j 0 (E/E o y fL , /* is between 3.5 and 4. In the following section the details of 
this spectrum are more carefully studied. 

Also in Figure 3 we display a schematic cross section of our instrument as it is 
deployed on the surface of the Moon. The apertures of both analyzers are elevated 
26 cm from ground. Their geometry is such that they observe only electrons with a 
component of velocity in the downward direction. Since we continually observe 
photoelectrons with energies up to -200 eV, we must assume that the lunar surface 
potential is at least 200 V during these times. This measurement will seem high to 
those familiar with the work of Wal bridge (1970) and Grobman and Blank (1969), 
who calculate a surface potential that is at least an order of magnitude lower. The 
difference is that their models deal with photoemissions from the surface of the Moon 
in the presence of the solar wind. Our measurements in the magnetotail are made 
under near-vacuum conditions. After further analysis of the problem we return to 
considerations of the surface potential. 

To summarize: During geomagnetically quiet times, when the Moon is in the 
magnetotail and not in eclipse, stable photoelectron fluxes with energies between 40 
and 200 eV are observed. These fluxes are isotropic and obey a power law, E where 
ix is between 3.5 and 4. From the fact that CPLEE is observing downward moving 
electrons w r e conclude that in the high-latitude magneto tail the lunar surface 
potential is on the order of 200 V, 

4. Numerical Analysis 

4.1. General theory 

The variations of photoelectron density and electrostatic potential above the surface 
of the Moon can be calculated numerically. We approximate the lunar surface by an 
infinite plane, with the x direction normal to the surface, and assume spatial variations 
of physical quantities only with the height. 

At a height jc above the surface the electron density is J/ (v, x) / (v, x) is the 
electron distribution function. If we assume an isotropic flux at the surface, the 
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Liouville Theorem can be used to show that the distribution function is independent 
of angles at all heights. Writing 


d 3 i> = 


2 E 


d£d*2 


m~ 


and integrating over solid angles, the density is 


n(x) — An J J 2 ^f(E,x)dE. 
0 


( 1 ) 


Since the distribution function is a constant along particle trajectories, / (£, = 

= / (E 0J x — 0), where E—£ 0 — g[(p 0 ~ </>(*)]. By changing the variable of integration 
from E to E 0 Equation (1) can be expressed 

00 

n(x) = 4x J Jim (£ 0 - q [><, - q> (x)]) / (£ 0 , x = 0) d£ 0 • (2) 

“«►(*)] 

To calculate the distribution function of photoelectrons at the surface consider the 
quantity 

j (E 0 ) dE 0 = [ J / (hv) Y (hv) e (E 0 ,hv) d Avj d£ 0 (3) 

the upward moving flux of photoelectrons emitted from the surface with energies 
between E 0 and E 0 + d£ 0< I (hv) d (hv) is the flux of photons reaching the lunar surface 
with energies between hv and hv +d (Av). T(ftv), the quantum yield function, gives the 
number of electrons emitted by the surface per incident photon with energy hv . 
q(£ 0 , hv) d E 0 is the probability that an electron emitted from the surface, due to a 
photon with energy hv , will have a kinetic energy between E Q and £ 0 -f d E 0 . q(E 0 , hv) 
is normalized so that 


DO 

j q (E 0 , hv) d E 0 = 1 . 

0 

co 

The total upward moving flux at the surface is S 7 (x=0)- j j(E 0 ) d E 0 . But 

o 


2 n ffi/2 (C 


s t (x = o) = JJJ v 0 /(£ 0 , 0 , cp , 0) dt> 0 sin 0 dtf d<? . 


ooo 


Since v 0 — y 0 [i cosO+j sintf cos<p+k sin# sin cp] and /is independent of angle. 
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Thus 


and 


/(E 0 ,x=0) = 


m 2 j(_E o) 
2nE 0 


n (x) = 2 


I 


■j(Eo) 


Jim (£„ - q Oo “ </> (*)]) ' ■ d£ y . 

^0 




(5) 

( 6 ) 


The potential as a function of height is evaluated by multiplying Poisson’s equation, 
d 2 (pjdx 2 - —4 Ttqn(x), by d(p/dx and integrating in from x=oo to get 


o 



<pix) 


n((p')d<p'. 


where we have written 



(7) 


A further integration out from the surface, gives us the potential at a point x. 


4.2. Computational methods and results 

To determine the upward moving differential flux at the surface, upon the knowledge 
of which the distribution function, number density and potential depend, we must 
first calculate the integral in Equation (3). The solar photon differential flux at I AU, 
l(hv), is taken from Friedman (1 963) for the range 2000 to 1 800 A and from Hinteregger 
(1965) from the range 1775-1 A and is plotted in Figure 4. Following the suggestion of 
Walbridge (1970) we have: 

(1) Adopted a work function W of lunar material of 6 eV. 

(2) Assumed a photoelectron yield function of the form 

0 hv < 6 eW 

. hv — 6 

Y(hv) = Y 0 —^— 6 0v<9eV (8) 

y 0 hv^ 9eV, 

w'here Y 0 is a free parameter of our calculation. 

(3) Chosen a probability function 


e(E, hv) = 


'6E(E l - E)jE v 
0 


0 < £ < Ey 

E>E l9 



where 


hv ^ W 
hv < W . 


(9) 
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In general the probability function is a complicated function depending on the 
nature of the photoemission material. However, Grobman and Blank (1969) have 
shown that for the purpose of calculating Equation (3) any broad function with zeros 
at £= 0 and E=E t and a width AE~hv will suffice. A plot of q{E , hv) is shown in 
Figure 5 for various value of E L . 



Fig. 4. Solar photon energy spectrum at 1 AU from 2000 to 1 A. 


The upward directed differential flux in electrons cm" 2 s" 1 sr" 1 eV" 1 for the values 
Y 0 = l, 0.1, 0.01 were numerically computed and have been plotted in Figure 6. We 
have also inserted the photoelectron differential flux observed by CPLEE at 26 cm. 
The Liouville theorem allows us to set a lower bound on Y 0 of 0. 1 . That is if there were 
no potential difference between the surface and 26 cm the yield function would be 0.1 
electrons photon" 1 . After estimating the potential difference between 26 cm and the 
surface we can. also determine an upper bound on Y 0 . 
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Solving the integro-differential Equation (7) for (p{x) involves an integration from 
the surface outward, with an assumed value of <p 0 . However the expression for d(p/dx 
involves an integral from infinity in to jc, or equivalently from <p= 0 to <p(x). By the 
expedient of dividing the integral into pieces in E 0 space and using an analytic 
approximation to the function ; (E 0 ) in each of these intervals, a solution was effected. 



Fig. 5. Probability function that a photon of energy hv will cause the lunar surface material to emit 
a photoclectron of energy E with different values of E\ — hv — W , 


111 this way it was only necessary to know the values of (p~cp{x) and 0 at the end 
points of the interval, and the solution would proceed. In Figure 7 we show families of 
solutions for <p(x) with several values of the parameter Y 0 . 

The value of Y 0 calculated by assuming no potential difference between the surface 
and x = 26 cm was 0.1. Figure 7 shows that for Y o = 0,l, the potential difference 
&(x=Q) — <2>(x-26 cm) is only 3 V. Obviously, we could now use an iterative pro- 
cedure, modifying our spectral measurement at 26 cm to obtain the surface spectrum 
according to the Equation /(£, x)=f ( E 0 , 0) and hence obtain a new estimate of Y 0 . 
However, the procedure is hardly justified considering the small potential difference 
( — 3 V) and the energy range of the measured photoelectrons (40-200 eV). Hence we 
conclude from our measured photoelectron fluxes and numerical analysis a lunar 
surface potential of at least 200 V and a value of the average photoelectron yield of 
Y o = 0.i electrons photon" 1 . 
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Fig, 6. Numerically computed photoclectron spectra emitted for the yield functions Yo = 1, 0.1 and 
0,01 electrons per photon. The photoelectron spectrum measured by CPLEE is found to fall close 

to the Yo =0.1 line. 


For the sake of completeness, Figure 8 shows the variation of potential and photo- 
electron density and pressure with height obtained from the numerical solution 
discussed above. For this solution, a value of y o = 0.1 and <p 0 = 200 V were chosen. 
It may be tempting to define a Debye length for the photoelectron gas according to the 
formula 


2 %n 


\l/2 

V 


fY 




= 4 cm at x = 1 cm . 



THE LUNAR PHOTOELECTRON LAYER IN THE GEOMAGNETIC TAIL 


379 


However, this characteristic length is by no means the distance over which the entire 
potential drop is developed. It is important to keep in mind that the photoelectron 
layer is not a plasma, but rather is a one-component gas and hence the concept of a 
Debye length as a potential shielding distance is not applicable. 



HEIGHT (CM) 

Fig. 7. Numerically computed potential distribution above the lunar surface for several values of 
the yield function Jo. For Yo ■■ 0. 1 the potential difference between ground and 26 cm is about 3 V. 


5. The Lunar Surface Potential cp 0 

The experimental measurements of photoelectrons at 200 eV but no significant fluxes 
in the next highest energy channel at 500 eV lead us to conclude that the lunar surface 
potential is at least 200 V. The data of Hinteregger et al. (1965) show's significant solar 
photon fluxes up to 400 eV, and presumably the lunar surface potential under vacuum 
conditions could be 400 V. However, we detected no photoelectrons with £>400 eV, 
and in fact the extrapolation of our measured spectrum (Figure 3) to 400 eV is below 
the instrument background. For this reason therefore we have adopted a conservative 
value of 200 V as the lunar surface potential for purposes of the calculations in the 
preceding and following sections. 
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The lunar surface potential can be decreased however by the presence of a hot 
ambient plasma w hich furnishes an electron return current which partially balances the 
emitted photoelectron current. In effect, the highest energy photoelectrons can escape 
from the potential well, since electrons from the ambient plasma furnish the return 



Fig. 8. Numerically computed values of electron density, potential difference and pressure as 
a function of distance from the lunar surface. 


current to balance these escaping photoelectrons. Quantitatively, if is the net 
negative flux to the lunar surface from the ambient plasma, and j(E 0 ) is the emitted 
photoelectron energy spectrum in units of electrons cm" 2 s“ l eV' 1 , then: 

<x> 

F s =jj(E 0 )dE 0 (ip) 

Ell > Q(pa 

and this equation can be solved for <p 0 , the lunar surface potential. The results of this 
calculation for 30 V<<P 0 < 200 V are shown in Figure 9. The curve was computed for 
y 0 =1.0, but can be scaled for other values of Y 0 . 

Our measurements of photoelectrons were taken during periods in the magnetotail 
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when all of the channels of the instrument except the lowest-energy electron channels 
were at background levels. Thus we can establish an upper limit to the electron flux 
from the ambient plasma for electrons with 40 eV<lj<50 keV. Figure 3 shows the 
‘background spectrum’, calculated by converting the background counting rate of 



Fig. 9. Computed curve, from the emitted photoelectron spectra shown in Figure 6, of the lunar 
surface potential versus incident electron flux for Yo= F0. The curve can be scaled directly for 

any other value of Fo. 


— J count s' 1 to equivalent flux in each of the energy channels. Integrating over this 
spectrum and converting to flux over the hemisphere gives F s ^3.4xl0 e electrons 
cm -2 s' 1 . We feel that this is a valid upper limit, as the range of measurement in 
energy includes the peak energy of the plasma sheet spectrum (~1 keV). 

We note that Vasyliunas (1968) obtained an upper limit to the electron concentra- 
tion for locations outside of the plasma sheet based on OGO-3 data. The relation 
expressed was NE^ 2 < 10~ 2 cm' 3 keV 1/2 where N is the electron density and E 0 is the 
energy at the peak of the spectrum. For an isotropic plasma where the bulk motion 
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can be neglected relative to the thermal motion, the electron flux to a probe is given 
by = Nvjln - . Applying the appropriate conversion of factors, the expression of 

Vasyliunas results in an upper limit to the electron flux of F A <5.6xl0 7 electrons 


TABLE l 


Electron flux 

(cm~ 2 s _1 ) 

Yu 

tfo(V) 

3.4 X 10 6 

1.0 

181 

3.4 X 10« 

0.1 

114 

3.4 x 10 fi 

0.01 

44 

5.6 x TO 7 

1.0 

96 

5.6 x 10 7 

0.1 

36 

5.6 x 10 7 

0.01 

8 


In Table I we show results of surface potential computations for the two electron 
flux upper limits given above and for values of Y 0 of 1.0, 0.1, and 0.01. 

The lower half-height of the channel 5 energy passband is 160 eV. Hence the surface 
potential could be as low as 160 V and still result in particle fluxes in channel 5. 
This estimate of the potential is seen to be not inconsistent with a value of T o = 0.1, 
/^<3.4x 10 6 resulting in a surface potential (Table I) of 114 V. 

One rather obvious prediction of our arguments about the surface potential is that 
when the electron flux reaching the surface of the Moon is sufficiently high, surface 
generated photoelcctrons with energies in the range of our detector should vanish. 
This happens in the solar wind and magnetosheath, but it is impossible for CPLEE 
to provide conclusive observational evidence in the presence of contamination by 
solar wind and magnetosheath electrons that the £'>40eV photo electrons are not 
returning to the surface and hence entering the detectors. Electron densities >1 cm -3 
and temperatures ~ 15 eV such as those commonly encountered in the solar wind and 
magnetosheath provide much higher fluxes in the 40-200 eV range than photoelectron 
fluxes observed in high latitude magnetotail. 

On April 9, 1971 a world-wide magnetic storm was observed by CPLEE (Burke 
et aL, 1972) in which the magnetosheath moved in to T SE =15 R E , As the magneto- 
pause moved out past the Moon, intense plasma sheet fluxes were observed. At this 
time the fluxes observed in the 40, 50 and 70 eV channels fell below photoelectron 
levels. In Figure 10 we show the count rates for the 40 eV and 500 eV electron channels 
from 1 1 :00 to 12:00 on April 9. A heavy line has been drawn at the 40 eV photoelec- 
tron level. As the 500 eV count rate rises the 40 eV count rate falls below this line. 
From 1 1 :28 to 11:35 when the 500 eV count rate dropped the 40 eV channel returned 
to the pholoelectron level. In Figure 11 we have plotted the electron spectra at 
11 : 26:40 (photoelectron) and at 11:27:19 (plasma sheet). The total incident fluxes, 
calculated by subtracting the photoelectron contribution and by assuming isotropy 
over the upper hemisphere are indicated for each of the two spectra. The data show that 
the lunar surface photoelectron yield in the range 40eV<£'<200 cV is of sufficient 



COUNTS / CYCLE 
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magnitude to maintain the potential at 40 V for an incident flux of 2.5 x ]0 7 elec- 
trons cm~ 2 s~ l . This can be used with the curve of surface potential vs incident 
electron flux (Figure 9) to compute a value of F o = 0.3. The discrepancy between this 
value and the value calculated previously (T 0 =OT) is probably due to an error 
associated with the assumption of isotropy of the incident flux. If the incident flux 
were smaller at large zenith angles, the total flux and hence the value of Y 0 computed 
by this method would be correspondingly less. 



Fig. 11. Electron spectra just prior to and after the photoelectron decrease at 11:27 shown in 
Figure 10. These spectra show in detail the modulation of the lunar surface potential by an incident 
electron flux impinging upon the lunar surface. The total fluxes, F*, were calculated by first sub- 
tracting the photoelectron contribution and then assuming that the higher energy electron flux was 

isotropic over the upper hemisphere. 


6. Discussion 

Feuerbacher et al . (1972) have measured the photoelectron yield of a lunar fine sample. 
In the photon energy range 5<£<20eV they found a yield function which reaches 
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a maximum value of ~0.08 at 15 eV then drops to 0.01 at 20 eV. B. Fitton (private 
communication, 1971) has suggested that our value of 0.1 js more a measurement of 
the CPLEE instrument case yield function than that of the lunar surface. We find it 
difficult to understand how this could be the case. 

First, were CPLEE an electrically isolated package then the requirement that the net 
current to the instrument be zero would result in a measured yield function that is 
representative of the case material. Relative to the lunar surface the case would bear 
a positive potential in order to maintain an enhanced photoelectron density in its 
immediate vicinity. 

The fact is, however, that CPLEE is not electrically isolated. CPLEE is connected 
to the central ALSEP station. Further, at any given time, only the top surfaces and 
one or two sides of ALSEP packages are illuminated by the Sun, while the remaining 
area is shadowed. These unilluminated surfaces provide receptor areas for return 
current from the photoelectron gas. Thus no large potential difference can develop 
between CPLEE and the lunar surface, and photoelectrons emitted at the lunar surface 
and at the CPLEE case are indistinguishable. Geometrical considerations of electron 
trajectories would lead us to expect that the bulk of the photoelectrons measured by 
CPLEE were emitted at the lunar surface at least several meters from CPLEE. If the 
ALSEP instrument cases had a photoelectron yield much larger or much smaller than 
the yield of the lunar surface, then one would expect a perturbation of the photo- 
electron flux in the vicinity of ALSEP. This perturbation would depend not only on 
the ratio of photoelectron yield, but more importantly on the ratio of the area of the 
ALSEP instruments to the area encompassed by the trajectories of £>40 eV photo- 
electrons from a point source. This last area is on the order of the square of the scale 
height of photoelectrons with £>40 eV (~ 10 m, Figure 11). Since the ALSEP area 
is ~3 m 2 , the ratio of areas is on the order of 2%. Thus, even if the yield of the instru- 
ment cases was a factor of 10 greater or smaller than the surface yield, the flux 
perturbation would only be 20%. 

Second, were CPLEE measuring its own photoelectrons one would expect to 
observe changes in the relative fluxes observed in the two analyzers with solar zenith 
angle as the electron cloud surrounding the instrument adjusts to changing illumina* 
tion conditions. Specifically, the ratio of the flux observed analyzer B (looking 60° west 
of vertical) to the flux observed in analyzer A (looking to the vertical) should be larger 
after than before lunar noon. Our data shows isotropic photoelectron fluxes across the 
entire magnetotail. 

1 * Summary and Conclusions 

In this paper we have reported the observation of stable, isotropic photoelectron 
fluxes 26 cm above the lunar surface. In the energy range 40 ^£^200 the flux obeys 
a power law of the form j (E)—j 0 (E/E 0 )~ fi where )x is between 3.5 and 4. Because 
these fluxes were moving down we conclude that in the near vacuum conditions of 
the high latitude magnetotail the lunar surface potential is at least 200 V. The modula- 
tion of the surface potential in the presence of intense plasma sheet fluxes has also 
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been observed. It was shown that these electrons can be explained in terms of the 
measured solar photon spectrum producing an isotropic flux of photoelectrons at the 
surface. A photoelectron yield function of Y 0 =0. 1 electron photon -1 was calculated. 
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DISCUSSION 

Feuerbacher: Did you take into account the photoemission from your box? 

Reasoner: The resolution of this question is a matter of the contact between the instrument cases 
and the local surface potential. It is certainly true that if the instrument case were truly isolated, 
then AJ= 0 would immediately imply that the measured photoelectrons were indicative of the 
instrument case yield. However, there is more area in electrical contact with the instrument that is 
not illuminated than is illuminated, and therefore these ‘cold’ areas can act as receptor areas to collect 
current from the photoelectron gas, and thus prevent large potential differences from developing 
between the instrument and the lunar surface. Granting this to be the case, then, and recalling that 
photoclectrons from both a vertical direction and a direction 60° from vertical were measured with 
equal intensities at all Sun angles, we therefore conclude that the measurement is a true representation 
of the lunar surface photoelectron flux. 

Walker: Let me get this clear: what is the basic discrepancy between Reasoner’s measurement 
and the laboratory measurements of the ESTEC group? 

Feuerbacher ; I do not think that there is a basic discrepancy between our results and those of 
Dr Reasoner. Our surface potentials have been calculated for solar wind conditions, while the CPLEE 
data refer to measurements in the high latitude magnetotail, where the incoming plasma flux is 
negligible. Our experimental data cannot possibly be extrapolated to electron energies of 200 eV, 
since the highest photon energy used was 23 eV. 

Fitton: T would like to refer to Figure 1. Could you indicate the likely effect of the cover sheet, on 
the right hand side of the figure, so far as the properties of the local photoelectron sheath measured 
by your instrument is concerned? 
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Reasoner: The photograph of Figure 1 shows the CPLEE, the ALSEP central station and a portion 
of the thermal blanket of the passive seismic experiment (PSE). This thermal blanket is to the best of 
my recollection composed of many layers of aluminized mylar and is located approximately 1.5 m 
from CPLEE. The top surface of the blanket is, however, not in electrical contact with the ALSEP 
system. Furthermore, photoelectrons emitted from this surface would produce highly anisotropic 
fluxes at the instrument, and this is not observed. 
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The purpose of this letter is to report the 
simultaneous observations of plasma flowing 
away from the earth and a southward-directed 
magnetic field across the neutral sheet at lunar 
distance during two weak magnetic bays. These 
simultaneous observations are in accord with 
the predictions of magnetic merging models for 
the generation of magnetic substorms. 

Plasma observations were made by the lunar- 
based charged particle lunar environment ex- 
periment (CPLEE) [cf, O'Brien and Reasoner , 
1971; Burke and Reasoner , 1972], The instru- 
ment contains two identical charged particle 
analyzers, one looking toward the local lunar 
vertical (analyzer A) and the other looking 
60° from the vertical toward lunar west (ana- 
lyzer B). Each analyzer measures a 15-point 
spectrum for protons and electrons with ener- 
gies between 40 ev and 20 kev every 19.2 sec. 
At 0000 on February 10, 1971, the moon was 
located at (-.V, Y , Z — ( — 62.8, 2.9, 2.5) Rb; 
approximately 5 R E above the Russell-Brody 
neutral sheet. Magnetic field data are from the 
Ames Research Center magnetometer aboard 
Explorer 35 [Mihcdov et al. } 1968], and 81.8-sec 
averages of the data are displayed in solar 
magnetospheric coordinates. 

Between 2200 on February 9 and 0200 on 
February 10, 1971, geomagnetic conditions were 
fairly quiet ( Kp — 2— f 2). Away from the 
midnight sector, auroral zone stations showed 
no signs of activity. Magnetograms from Leir- 
vogur and Kiruna, near local magnetic mid- 
night, indicate that there were two small 
substorms (Figure I). A hint of the so-called 
growth phase of the first substorm [ McPherron , 
1970] is contained in the Kinma magnetogram 
at 2245. The expansion phase began simultane- 
ously at Leirvogur and Kiruna at 2345, reached 
a maximum depression of —100 y at 2358, and 
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recovered to predisturbance conditions by 0020. 
The second event, beginning at 0116, was ob- 
served only at Leirvogur as a bay with a 100-y 
depression. 

Magnetic field data and the plasma anisot- 
ropy ratio observed at lunar distance during 
this time arc displayed in Figure 2, where B is 
the field strength in gammas, y> and 0 are the 
solar magnetospheric longitudinal and latitudi- 
nal directions of the vector, and a is the 
anisotropy, defined as the ratio of the differen- 
tial fluxes measured in the 2-kev proton chan- 
nels of analyzers A and B. Values of a. > 1 
indicate that plasma is flowing away from the 
earth, a < 1 indicates that there is a dawn to 
dusk flow across the tail, and a = 1 is the 
condition for isotropy. 

Between 2200 and 2340 the moon was in the 
southern lobe of the high-latitude tail ( B — 10 y, 
<p = ISO 0 ). The decrease in field strength and 
change in field direction at 2342 show that the 
moon moved rapidly in relation to the plasma 
sheet to the vicinity of the neutral sheet. At 
2348 the component of field across the neutral 
sheet B z changed from a northward to a south- 
ward direction. As B z changed direction, the 
plasma anisotropy ratio changed from a = 0.3 
at 2345 to a = 30 at 2350. During the period 
when the field was predominantly southward 
the anisotropy ratio was always greater than 5 
and reached values in excess of 40 on two 
occasions. By 0020 the moon moved away from 
the neutral sheet, and norma] plasma sheet fluxes 
were observed again. From 0030 to 0110 the 
moon was in the northern lobe of the high- 
latitude tail {B = 11 y, (p zr 0°). 

Correspondence between plasma and mag- 
netic field data during the second event (0110— 
0200) is obscured by a magnetic field data gap 
(0140-0150) and at least seven crossings of the 
neutral sheet. Strong anisotropies were again 
observed in the proton fluxes. However, the 
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Fig. 1. Substorms observed in northward com- 
ponent of magnetic field at Kirtma and the 
horizontal component at Leirvogur on February 
9-10, 1971. Tlie geomagnetic latitude, local time, 
and universal time are indicated. 

proton burst beginning at 0121 preceded the 
field change by about 3 min. Between 0130 and 
0135, while the moon was still in the neutral 
sheet, B z turned to the north, and a was mea- 
sured as less than 1, The final burst at 0140 
began at the same time as the field data gap. 
When magnetic field data resumed, the moon 
was away from the neutral sheet, in the plasma 
sheet, where isotropic fluxes were observed. 

Protons observed by CPLEE had a density 
of ^0.2 cm -3 , a temperature of ^1 kev, and a 
bulk speed of 250 ± 50 km/sec. The flow is 
assumed to be from the look direction of ana- 
lyzer A. This bulk speed in the presence of a 
southward-directed field of 2 y indicates a dawn 
to dusk electric field of 0,5 ± 0.1 mv/m. Qur 
measured plasma parameters are close to those 
reported by Prakash [1972]. The discrepancy 
with his calculated electric field of 0.12 mv/m 
is due to his assumption of \B Z [ = 0.6 y. On the 
other hand, our observed fluxes are between 1 
and 2 orders of magnitude less than those 
reported by Garrett et at. [1971]. Since their 
data were observed during substorms with de- 
pressions of 2000 y, it is quite possible that 
most of the flux differences only reflect the 
difference in intensity of the triggering mecha- 
nism. 

Our observations are consistent with the 
predictions of magnetic merging models proposed 
by a host of magnetospheric physicists. Plasma 


flow away from the earth while the magnetic 
field was pointing southward indicates that the 
neutral point was between the earth and the 
moon. The correlation between B z and a at 
~0130 reflects motion by the neutral point 
along the sun-earth line. 

Although the Kiruna magnetogram closely 
follows the McPherron [1970] substorm model 
(namely, the growth, expansion, and recovery 
phases), available data do not allow us to re- 



2200 2300 OOOO OIOO 0200 

HOURS {UT) 

Fig. 2. The magnetic field observed by the 
Arnes Research Center magnetometer aboard Ex- 
plorer 35 and the anisotropy ratios observed by 
CPLEE from 2300 to 0200 UT on February 9-10, 
1971. The field directions are given in solar mag- 
netospheric coordinates. 
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solye controversies regarding magnetospheric 
convection during the growth phase of the 
substorm. Between 2245 and 2345 the moon 
was in the high-latitude tail. The apparent 
anisotropy observed near 2300 is due more to 
Low counting rates than to a clearly identifiable 
flow of plasma. Flow away from the earth 
began at 2348 simultaneously with the change 
in B z from north to south. Whether these 
changes in plasma conditions reflect motion 
of an already existing neutral point from be- 
yond the lunar orbit or the creation of a 
neutral point between the earth and the moon 
remains an unsolved problem! 
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